北秦岭华阳川地区重磁三维反演及岩浆岩特征研究

图1 陕西华阳川地区地质背景及物探剖面位置示意

Fig.1 Geological background map and geophysical profile location diagram of Huayangchuan Area, Shaanxi Province

2 重磁数据处理及物性统计

2.1 数据来源及重磁数据处理

2.1.1 重力数据

2010~2011年陕西省地矿局第二综合物探大队在本区开展了1∶10万重力测量,2017~2018年中国地质调查局西安地质调查中心完成了工区部分1∶5万地面和航空重力测量,以上数据构成了本次建模的重力数据基础。

以地面重力数据为基准,在同一位置反映相同地质体的航空重力和地面重力数据中,选取均匀分布的控制点,通过多元回归分析法,建立航空重力和地面重力的转换关系;以地面重力数据为基准,将航空重力数据换算到地面,并利用地面1∶10万数据进行补空和扩边,得到本文所需的重力数据。根据本次4 km深度建模的目的,将数据网格化,网格间距为100 m,采用滑动平均滤波、趋势分析法、垂向二阶导数法和解析延拓法对华阳川地区重力数据进行了区域场与剩余场分离,经与地表地质对比,认为解析延拓方法(上延2 km)求取的剩余重力异常满足后续反演计算的要求(图2a)。

图2

图2 研究区重(a)、磁(b)异常

Fig.2 Gravity(a) and magnetic(b) anomaly map in the study area

2.1.2 磁测数据

本次采用的磁测数据为2016年中国自然资源航空物探遥感中心在华阳川矿床及其外围完成的 1∶1万高精度航空磁测数据。利用Geosoft Oasismontaj地球物理软件进行中高纬度的化磁极处理,其中输入的地磁倾角为52.46°,地磁偏角为-3.37°,并将化极异常作为三维反演的输入数据(图2b)。

2.1.3 电法数据

2018~2019年中国地质调查局西安地质调查中心在本区实施完成了大地电磁、广域大地电磁、AMT、短偏移距瞬变电磁、频谱激电5种电法手段,共16条,其剖面位置见图1。电磁法剖面同已知的地质图、钻孔、实测地质剖面信息作为初始模型,加入到三维反演模型中,对三维反演计算进行约束和控制。

2.2 物性特征

物性是连接地球物理与地质的纽带和桥梁,笔者对工区1137块物性标本的密度、磁化率进行了统计分析,为后续2.5维和三维反演计算提供依据。物性统计结果见表1和图3。

表1 地层、岩体密度及磁化率

Table 1 Statistical table of formation, rock mass density and susceptibility

地质代号	岩性	块数/块	密度/ (g·cm^-3)	磁化率/ (10^-5 SI)	地质代号	岩性	块数/块	密度/ (g·cm^-3)	磁化率/ (10^-5 SI)
ηγKD	中粗粒含斑黑云二长花岗岩	37	2.63	1430	βμ	含黄铁矿辉绿岩脉	24	2.99	4558
ηγKH	中粒含斑黑云二长花岗岩	32	2.57	590	Chd	黄铁矿、辉钼矿矿石	22	2.92	113
ηγKX	中细粒含斑黑云二长花岗岩	31	2.56	565	Chd	含辉钼矿的石英岩	41	2.66	8
ηγKS	细粒含斑黑云二长花岗岩	18	2.65	1780	Chl¹	细碧岩、细碧玢岩	70	2.85	4110
ηγKY	中粒黑云二长花岗岩	38	2.58	980	Chb¹	灰白—紫灰色石英砂岩	31	2.63	8
ηγKC	细粒黑云二长花岗岩	31	2.62	795	Chb²	灰色粘土板岩	25	2.66	10
Ar₂Dgg	花岗片麻岩	93	2.70	155	Chb³	紫—紫灰色石英砂岩	33	2.62	4
Ar₂Mgn	黑云角闪斜长片麻岩	68	2.68	338	Chb⁴	石英砂岩夹粘土板岩	32	2.59	4
Ar₂Hgn	角闪花岗片麻岩	75	2.65	395	Chb⁵	灰白色石英砂岩	32	2.62	7
Ar₂Wgn	黑云斜长片麻岩	106	2.62	1425	Chb⁶	灰绿色—青灰色板岩	31	2.58	9
ηγPt₂G	细粒二长花岗岩	47	2.61	395	Che¹	肉红色石英砂岩	32	2.64	5
Jxl	浅灰色白云岩	32	2.83	2	Chch	紫色泥质板岩	31	2.74	18

新窗口打开| 下载CSV

图3

图3 研究区岩体、地层密度及磁化率散点分布

Fig.3 Dispersion distribution of rock mass, formation density and susceptibility in the study area

华阳区地区岩(矿)石物性具有以下特征:

1) 华山岩体和老牛山岩体具有低密度、较高磁性的特征;

2) 长城系老地层为低密度无磁或弱磁性,另熊耳群洛源组的细碧岩、细碧玢岩为高密度(2.85 g/cm³)高磁特征;

3) 太华岩群片麻岩岩套为中等密度、中等磁性特征,其中翁岔铺片麻岩套的密度低于大月坪片麻岩,磁性比大月坪片麻岩略强;

4) 燕山期侵入岩体(ηγKX、ηγKH、ηγKY、ηγKC、ηγKD、ηγKS)密度和磁化率呈正相关性(图3黑色虚线),分析其岩性可知,岩性均为含斑黑云二长花岗岩,但ηγKS为细粒含斑黑云二长花岗岩,而其他为中粒或中—细粒含斑黑云二长花岗岩,说明同一岩性的岩体,其粒度小且致密,其密度相对较大,磁化率相对较高。

5) 太古宇片麻岩套(Ar₂Dgg、Ar₂Mgn、Ar₂Hgn、Ar₂Wgn)的密度和磁化率呈负相关性。大月坪片麻岩的岩性为花岗片麻岩,其岩性单一,主要矿物为弱磁特性,随之岩石经过不同程度的变质、分化剥蚀,变为黑云斜长片麻岩、角闪斜长片麻岩和黑云斜长角闪斜长片麻岩等,岩石中的磁性矿物(黑云母、绿泥石、绿帘石等)含量逐渐增加,致使其岩石磁化率逐渐增加,而岩石的紧密性也因多种矿物的充填而产生了微小的空隙,使的密度值变小。

研究区岩体、地层及矿石具有明显的物性差异,利用重磁异常可较好地划分各岩性单元的空间展布,为下一步三维地球物理—地质模型提供了地球物理前提。

3 重磁三维反演技术流程

集成已有地质、钻孔、岩性和地球物理资料,尤其是充分利用多种地球物理数据,包括航空、地面、测井等资料,通过数据融合、多方法组合等,建立工区三维地质—地球物理模型,实现对工作区三维地质结构的刻画。具体反演流程如图4所示。

图4

图4 三维地质建模技术路线

Fig.4 3D geological modeling technology route

3.1 初始模型的建立

主要包括对地表岩性单元或地质单元进行简化,钻孔数据、岩石物性测量,岩性与物性对应关系分析,利用回归分析法对不同尺度、空间的地球物理数据融合处理,对电磁剖面的推断解释,并将多元数据信息作为三维反演的初始模型,其中钻孔信息和电法剖面提供深部主要地层单位的边界深度,一般在重磁反演中作为重要的约束,保持不变。

3.2 2.5D地质—地球物理模型构建

本次剖面反演利用了剩余布格重力异常和化极磁异常的网格数据,对网格数据进行点距100 m,线距1 km的剖面数据提取,研究区地质构造和重磁异常多为NWW、NE向或近EW向,所以二维剖面方向取为SN向,共建立了34条人机交互反演剖面,从西往东线号依次为L3853、L3903、L3913、L3923、……、L4233、L4283,其中L3853和L4283为外扩剖面。

使用澳大利亚Tensor Research公司开发的Model vision软件对研究区34条剖面进行反演拟合计算,模型EW向延伸500 m,背景密度设置为2.67 g/cm³,背景磁化率为0 SI,反演深度为地下4 km,为了消除模型的边界影响,将模型体向东西南北向各延长了10 km,获得了该区2.5维的地质反演结果。图5展示了通过华阳川韧性剪切带L4083剖面的曲线拟合结果及地质解释。

图5

图5 L4083剖面2.5D反演解释结果

Fig.5 2.5D inversion interpretation result of the section L4083

本步骤实质为将地球物理异常转化为地质信息的过程,其推断和解释的地质剖面是否合理至关重要,因为该步的结果将作为后面三维地质建模的初始模型,将影响整个三维地质模型的结果。

3.3 重磁三维反演

把拟合的34条2.5维地质—地球物理剖面作为初始模型,利用Intrepid Geophysics公司的GeoModeller,采用“基于规则建模”的原则,遵循地层层序关系,进行综合三维地质—地球物理模型建模计算,并结合已知的地质图、钻孔及地质柱状图等信息对模型进行反复修改,再进行三维反演计算,直至得出与地质规律相符合的三维模型为止。在整个模拟过程中,物性与岩性的对应关系保持不变。

通过本次区域三维反演,可以提取区域内地层、构造、岩体的三维空间分布和深部地层、岩体之间接触关系,同时可提供深部成矿信息,结合成矿模型开展深部成矿预测。针对矿区三维建模,可以全面分析控矿地层、矿体和岩体的空间关系等,尤其是华阳川铀矿区的含矿地层和周边岩体的接触关系。

4 三维地质—地球物理模型解释

三维地质模型反映了丰富的地质信息,通过对三维地质模型的分析,可为区域构造、岩体分布、控矿构造华阳川剪切带的空间分布提供了丰富的信息。

4.1 区域地层、岩体的空间分布特征

通过三维地质模型(图6)可知研究区三维地层、岩体分布情况:在研究区的东南部地层从寒武系、震旦系、蓟县系、长城系到太古宇的太华岩群均出露完整,且形成向斜构造,向斜构造的轴部位于研究区的西南角,走向为EW向;研究区的西南部出露古元古界铁洞沟组。研究区的西北和东北部基本为岩体和太华群的片麻岩套覆盖。以华阳川韧性剪切带为界,界线的西南地区广泛分布老牛山复式岩体(包括侏罗系和白垩系岩体),岩体规模巨大,埋深较深,在本次反演的地下4 000 m处为见岩体根部。界线以北,分布华山岩体(包括葱峪岩体、大夫峪岩体、翁峪岩体、西峰岩体和杜峰岩体),根据三维反演结果显示,华山岩体在研究区内的规模较大,但埋深较浅。在麻地沟一带出露小河岩体(ηγPt₂G),规模和埋深均较小。研究区岩体与地层之间主要为断层接触和侵入接触,复式岩体内不同岩性单元之间呈涌动型侵入接触关系。

图6

图6 研究区三维地质模型示意(图例同图1)

Fig.6 Diagram of three-dimensional geological model in the study area(the legend is the same as Fig.1)

4.2 华山岩体、老牛山岩体与华阳川韧性剪切带三维形态与分布特征

在整个工区建立的三维地质—地球物理模型的基础上,为进一步精细刻画华阳川铀矿矿区及外围的三维地质结构,利用剩余重力异常进行三维密度反演(图7)。反演结果直观地展示出华山岩体、老牛山岩体和华阳川韧性剪切带的三维形态和分布。为验证本次三维反演的结果,并分析华山岩体、老牛山岩体和华阳川韧性剪切带的深部接触关系,截取了横跨华阳川铀矿区、老牛山岩体和华山岩体的Z2剖面,通过对比反演密度剖面和已有电阻率剖面,发现两者具有很好的对应关系(图8),说明本次三维反演的结果可靠且符合地质规律。通过图7、图8可知:

图7

图7 华山岩体、老牛山岩体和华阳川韧性剪切带的空间分布

Fig.7 Spatial distribution of Huashan rock mass, Laoniushan rock mass and Huayangchuan ductile shear zone

图8

图8 华山岩体、老牛山岩体与太华岩群深部接触关系

Fig.8 Deep contact relationship between Huashan rock mass, Laoniushan rock mass and Taihua rock group

1) 老牛山岩体和华山岩体在深部空间上并不联通,结合同位素年龄^{[15,16,17,18,19]},华山岩体侵入早与老牛山岩体,两者可能并不同源。华山岩体下覆地层为太华岩群,认为华山岩体源区为太华岩群。

2) 华阳川韧性剪切带大致呈“漏斗形”,NW走向,EW向长6 km,SN宽为0~200 m不等;在深部2 km时,宽度变窄至40~50 m。剪切带的密度随着深度变化逐渐变大,到深部2.5 km处接近太华岩群密度。且2018年实施的广域电磁剖面结果,佐证了本次推断华阳川韧性剪切带的形状及规模。

3) 结合同位素年龄,华阳川韧性剪切带形成于400 Ma,而华山岩体、老牛山岩体大约在140 Ma侵入就位,演化过程中,韧性剪切带受岩体侵位时挤压而变形;在中—新生代以来,由于华山岩体、老牛山岩体发育节理及脆性断层,早期韧性剪切带由叠加脆性变形再次被活化。

5 结论

本文尝试在地形地质“双复杂”地区开展了区域和矿区两个尺度的三维反演与建模工作。集成航空、地面和测井等不同来源、不同属性、不同维度的数据资料,通过多方法数据整合处理、多元数据融合、数据验证、联合反演等工作,构建了工作区三维地质—地球物理模型,获取了华阳川地区深部矿体、地层、构造、岩体的空间结构与分布特征,为华阳川铀矿区勘探开发与深、边部找矿提供地球物理依据。

1) 从区域重磁三维反演结果可知:研究区内地层、岩体等的深部的形态和空间分布规律,有助于地质人员快速了解和掌握区内的地质背景。

2) 由华阳川矿区及外围三维密度反演结果可获取华山岩体、老牛山岩体和华阳川韧性剪切带的三维形态和分布特征。结合地质资料,初步推断解释了华阳川韧性剪切带的形成过程,及与华山岩体、老牛山岩体的接触关系。并利用野外最新的广域电磁剖面佐证了本次三维反演的结果。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

高成, 康清清, 江宏君, 等.

秦岭造山带发现新型铀多金属矿:华阳川与伟晶岩脉和碳酸岩脉有关的超大型铀—铌—铅—稀土矿床

[J]. 地球化学, 2017, 46(5):446-455.

Gao

, Kang

Q Q

, Jiang

H J

, et al.

A unique uranium polymetallic deposit discovered in the Qinling orogenic belt: The Huayangchuan super-large U-Nb-Pb-REE deposit associated with pegmatites and carbonatites

[J]. Geochimica, 2017, 46(5):446-455.

[2]

赵海杰, 毛景文, 叶会寿, 等.

陕西黄龙铺地区碱性花岗斑岩及辉绿岩的年代学与地球化学:岩石成因及其构造环境示踪

[J]. 中国地质, 2010, 37(1):12-27.

Zhao

H J

, Mao

J W

, Ye

H S

, et al.

Geochronology and geochemistry of the alkaline granite porphyry and diabase dikes in Huanglongpu area of Shaanxi Province: Petrogenesis and implications for tectonic environment

[J]. Chinese Geology, 2010, 37(1):12-27.

[3]

毛景文, 谢桂青, 张作衡, 等.

中国北方中生代大规模成矿作用的期次及其地球动力学背景

[J]. 岩石学报, 2005, 21(1):169-188.

Mao

J W

, Xie

G Q

, Zhang

Z H

, et al.

Mesozoic large-scale metallogenic pulses in North China and corresponding geodynamic settings

[J]. Acta Petrologica Sinica, 2005, 21(1):169-188.

[4]

焦建刚, 汤中立, 钱壮志, 等.

东秦岭金堆城花岗斑岩体的锆石U-Pb年龄、物质来源及成矿机制

[J]. 地球科学:中国地质大学学报, 2010, 35(6):1011-1022.

Jiao

J G

, Tang

Z L

, Qian

Z Z

, et al.

Metallogenic mechanism,magma source and zircon U-Pb age of Jinduicheng granitic porphyry, East Qinling

[J]. Earth Science-Journal of China University of Geosciences, 2010, 35(6):1011-1022.

DOI:10.3799/dqkx.2010.114 URL [本文引用: 3]

[5]

孙晓明, 刘孝善.

金堆城钼矿区两类不同花岗岩的关系及其成因的研究

[J]. 地质找矿论丛, 1987, 2(2):34-45.

Sun

X M

, Liu

X S

Study on the relationship between two different types of granites and their genesis in Jinduicheng molybdenum mining area

[J]. Contributions To Geology and Mineral Resources Research, 1987, 2(2):34-45.