2. 自然资源部应用地球物理重点实验室, 长春 130026;
3. 中国地质科学院地质研究所, 北京 100037
2. Key Laboratory of Applied Geophysics, Ministry of Natural Resources, Changchun 130026, China;
3. Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China
双尖子山矿集区位于内蒙古巴林左旗北部,中心地理位置为44°30′N,119°07′E,在大地构造上位于大兴安岭中南段,处于兴蒙造山带北造山带和南造山带的夹持部位(徐备等,2013),成矿潜力巨大,是亚洲最大、全球第六大银矿床(张大权等,2015).根据矿体分布情况,矿集区分为双尖子山矿段和兴隆山矿段(匡永生等,2014).矿集区位于古亚洲洋构造域与古太平洋构造域交叉叠合的部位(刘建明等,2004),在古生代期间经历古亚洲洋的发生、发展以及闭合,由华北克拉通北缘和西伯利亚板块碰撞拼接控制(刘建明等,2004;Windley et al., 2007);到中生代,受古太平洋板块大规模俯冲以及后撤的作用,该区域发生了由挤压构造体系逐渐向伸展构造体系的转变(Zorin,1999;邵济安等,1999),引起区域强烈的多期次断裂、岩浆活动,使得该区域地质情况复杂,并具有优越的成矿条件.
近年来,前人依据双尖子山矿区地表地质、矿体出露特征及钻井资料,对矿区成矿要素进行详尽的研究:(1)矿区主要赋矿地层为二叠纪大石寨组(陈良等,2009);(2)区域断裂构造发育,呈北东向、北西向及近南北向,为矿集区提供良好的赋矿空间(匡永生等,2014;王丰翔等,2016);(3)与成矿相关的岩体主要为矿区深部晚侏罗世和早白垩世斑状花岗岩及正长花岗岩(吴冠斌等,2014;欧阳荷根等,2016).同时前人依据矿石及与成矿相关花岗岩的地球化学结果,得到矿区的成矿机制,认为中晚侏罗世蒙古—鄂霍次克洋的闭合作用以及早白垩世古太平洋板块俯冲及构造体制转换作用导致软流圈物质上涌成矿(Fan et al., 2003;Wu et al., 2007;王丰翔等,2016;顾玉超等,2017;王祥东,2017).但以上认识主要依据地球化学及钻探结果,研究深度局限于1000 m以浅,对于矿区成矿要素的深部展布特征以及相互间的接触关系并不明确.因此,在前人对于双尖子山矿区深部地质特征研究不足的背景下,针对深部主要成矿要素的探测以及刻画对于推动深部找矿具有重要意义.
地球物理勘探具有勘探深度大、分辨率高的特点,可提供矿区深部丰富的地质信息(Thomas et al., 2016;Zhao et al., 2018),因此,在矿区内展开地球物理勘探工作对深部找矿具有重要作用.其中可控源音频大地电磁测深法(CSAMT)是一种基于大地电磁法和音频大地电磁法发展起来的人工源频率域测深方法,具有抗干扰能力强、探测深度大、工作效率高等特点(何继善,1990;石昆法,1999).近些年来,CSAMT方法在金属矿区确定地层结构、隐伏矿体、有利成矿位置等方面取得了卓越的成效(于昌明等,1998;Hu et al., 2013;Guo et al., 2019;Zhang et al., 2020).目前,对CSAMT实测资料的反演解释主要采用不含源的一维和二维反演得到地下电性结构,但实际地下地质情况复杂,一维、二维反演方法很难得到可靠的地电结构(Ledo,2005;林昌洪等,2011).三维反演则可以克服一维和二维反演中存在的问题,是准确描述地下探测目标全方位空间展布特征的最有效手段(Siripunvaraporn et al., 2005;Commer and Newman, 2008;刘云鹤,2011),可用于揭示矿区成矿要素深部特征.
因此,为了得到更精细的成矿要素深部展布特征,建立双尖子山矿集区成矿模式,本次研究在双尖子山矿集区布设了6条5 km长的CSAMT剖面(图 1),采用CSAMT含源三维反演方法对数据进行处理,得到了地下2000 m以浅电性结构,并结合地质、物性和重磁数据对三维电性结构中的异常体进行约束和解释,得到矿集区深部成矿要素展布特征及成矿模式,为矿集区深部找矿提供方向.
双尖子山矿集区及外围主要由呈北东向展布的二叠系火山喷发-碎屑岩和侏罗纪火山-沉积岩系组成(吴冠斌等,2014;王丰翔等,2016;江彪等,2019).区域地层单元方面,矿集区主要出露地层为白垩系玛尼吐组(K1mn)、白垩系满克头鄂博组(K1mk)、侏罗系新民组(J2x)、二叠系大石寨组(P1d)以及第四系(Q)(匡永生等,2014);赋矿地层方面,矿集区80%的金属矿产于二叠系地层,其中广泛出露的大石寨组是该区域主要的银铅锌赋矿层位,大石寨组总体呈北西向单斜构造产出(陈良等,2009),大石寨组物性采集标本多由砂岩组成;断裂构造方面,矿区内断裂构造发育,主要发育为北西向(300°~310°)、北东向(30°~50°)和近南北向(0°~10°),表现为多期构造活动特征,呈交叉网格状,是矿集区重要的控矿构造;成矿物质来源,与矿化相关的两期岩浆事件分别为早白垩世以及晚侏罗世岩浆活动(吴冠斌等,2014;欧阳荷根等,2016),表现为中酸性岩浆侵入,与成矿有关的岩体主要包括花岗闪长岩、黑云母花岗岩、二长花岗岩等(张梅等,2011;杜青松等,2017).
1.2 岩石物性特征岩石物性可作为地球物理结果解释的重要依据.本文对内蒙古东南部巴林左旗地区1043块岩石物性样本进行密度、磁化率以及电阻率测试,获得了主要地层及侵入岩的物性特征,为地球物理解释工作提供依据,如表 1所示.
研究区内沉积地层呈现为中阻、较高密度和低磁化率的物性特征;侵入岩电阻率均高于地层电阻率,密度与地层密度相似,磁化率高于地层磁化率.侵入岩和沉积地层在电阻率和磁化率方面存在明显差异,这种差异可为研究区的地层和岩体识别提供重要的依据.
2 数据采集与处理 2.1 数据采集CSAMT测线分布如图 1,6条测线平行布设,每条测线长度为5000 m,测点点距为50 m,线距为1000 m,测线覆盖双尖子山矿床.考虑测区地形及采矿、工业噪声等干扰源的影响,对不同测线采用三个不同位置的发射源(源A1B1对应测线L1的发射源,源距为1.44 km,收发距为7.08 km;源A2B2对应测线L2、3的发射源,源距为1.52 km,收发距分别为8.31 km、7.31 km;源A3B3对应测线L4、5、6发射源,源距为1.45 km,收发距分别为17.08 km、16.08 km、15.08 km).CSAMT数据采集利用加拿大凤凰公司生产的V8多功能电法测量系统,配备AMTC-30磁探头和不极化电极,采用标量观测,观测频率范围为9600~1 Hz,采集平均时长为45 min,得到每个测点的视电阻率和相位.其中超过90%的CSAMT数据质量很高(图 2a和2b).对于数据质量较差的数据(图 2c)采取剔除的方式,不参与三维反演计算.
本次三维反演采用非结构四面体网格对初始地电模型离散化处理,并基于L-BFGS对数据进行反演计算(Broyden,1970;Fletcher,1970;Goldfarb,1970;Shanno,1970).由于相位数据整体干扰较大,本文仅对视电阻率进行反演计算,估计电磁噪声为观测数据振幅的10%.采集数据中大部分测点小于10 Hz的频点突跳严重,质量较差,同时考虑到区域电阻率值较高,利用9600~10 Hz范围内的数据可以揭示2000 m以浅电性结构(表 2),因此反演频率范围为9600~10 Hz.
由于实际地形起伏对于反演结果准确性影响较大,反演初始模型以非结构四面体网格进行剖分,可准确地拟合起伏地形的变化(George et al., 1990;刘长生,2009;吴小平等,2015).反演初始模型总体大小为100000 m×100000 m×100000 m.研究区为精细剖分区域,大小为7500 m×7500 m×5000 m,其中研究区深度500 m以上空间按四面体单元大小为100~200 m,增长率为1.4(增长率设置范围:1.3~2;1.4可满足精细剖分网格)进行网格剖分,深度500~5000 m空间按四面体单元大小为150~300 m,增长率为1.4进行剖分(图 3a和图 3b).发射源处按四面体单元大小为8~15 m,增长率为1.3进行加密剖分(图 3c),减少发射源附近电磁场急剧变化引起的场源奇异性对数值模拟精度的影响(Plessix et al., 2007).反演模型共包含989309个四面体单元.研究区及扩边地层电导率设置为500 Ωm,研究区外围空气电阻率设置为1×108Ωm.最终反演迭代179次,数据拟合差为3.59(图 4).图 5为三维反演后测点RMS分布情况,大部分测点的RMS在4以下,数据整体拟合情况较好.
图 6为最终反演后得到的三维电性模型,模型大致呈双层结构,即纵向具有成层性.上部为厚度在100~1000 m的低阻-中阻层,电阻率小于2000 Ωm;下部为电阻率在2000~8000 Ωm的中高阻层以及大于8000 Ωm的高阻层.
根据地表地质与三维反演结果(图 6)发现,大石寨组、玛尼吐组以及满克头鄂博组下方为中等电阻率特征,与岩石物性测试结果一致.对比双尖子山矿集区内井ZK12-37的钻井结果(Zhai et al., 2020)(图 7),发现第四纪砾岩(-197~-147 m)和二叠纪大石寨组板岩(-147~170 m)对应浅部中低阻异常,蚀变板岩(170~319 m)对应中高阻异常,花岗斑岩(319 m以深)对应高阻异常,表明依据电阻率差异可识别与成矿有关的地质要素.
依据井ZK12-37及地表地质的约束,对L2剖面进行地质解译(见图 7),其中800 m以浅的中阻区推断为沉积地层分布(主要为大石寨组,浅层存在满克头鄂博组、新民组),依据地表出露的地层分布特征以及地层倾向,推测出L2测线第四纪地层、新民组、满克头鄂博组及大石寨组分布(图 7);0~800 m范围内的中高阻异常推断为侵入岩与板岩发生绢云母化和硅化蚀变作用的反映;0 m以下的大面积高阻区推测为侵入岩(花岗斑岩)的分布特征.结合断裂所呈现的低阻-中低阻条带特征,识别出F2、F3、F4三条地表已知断裂和一条隐伏断裂F1(图 7).
根据此剖面的电性结果,结合井资料、地表地质资料与岩石物性测试结果,可以总结出研究区域CSAMT数据解译原则为:浅部中阻区(800 m以浅)为沉积地层;中部中高阻区(0~800 m)为侵入岩与二叠纪大石寨组发生蚀变作用的响应;下部中高-高阻区推断为侵入岩;剖面中低阻条带或低阻条带可能为断裂分布.
3.2 成矿要素展布 3.2.1 地层500 m以浅切片(图 8c—h)主要分为低阻、中阻和中高阻区.其中,大石寨组作为研究区广泛出露的沉积地层,根据钻井及岩石物性结果的约束,推断低阻、中阻区主要为大石寨组的响应.根据测线切片结果可明显发现,大石寨组在L1—L3测线下方呈现为由南东向北西厚度逐渐减薄的分布特征;而在L4—L6测线下方,大石寨组在测点0—9以及测点54—72范围内,厚度较薄,为200 m左右,由南东到北西呈现为厚-薄-厚-薄的分布特征.
结合L1—L6纵向剖面结果(图 6),研究区大石寨组主要分布在地下0~1000 m,其中在兴隆山矿段南部、双尖子山矿段东部,地层较厚,厚度范围在400~1000 m,可为含矿热液提供良好容矿空间.
3.2.2 断裂依据电性切片结果(图 8c—j)和断裂的电性识别特征,推断出9条北西向断裂、2条北东向断裂以及3条近南北向断裂,与矿集区内主要为北西向、北东向和近南北向的断裂分布特征一致(吴冠斌等,2014;顾玉超等,2017),其中除了北东向的F1深断裂外,其余断裂均为存在于500 m以浅的浅层断裂,同时所有断裂均与沉积地层下方的高阻侵入岩接触,可为含矿热液提供上涌通道及成岩成矿空间.
在兴隆山矿段,断裂表现为北东向、北西向以近南北向的浅部断裂(F2—F6).图 8b为兴隆山矿段地质特征,兴隆山东矿段矿脉由北西向韧性剪切带、北东向断裂以及南北向断裂控制,这三类已查明断裂位置与北西向F6、北东向F4以及近南北向F2、F3位置基本一致;北西向断裂F6在500 m以浅的结果中倾向为南西向,该特征与王丰翔等(2016)在图 8b中E-F位置利用钻井勘查得到该位置矿体受北西向韧性剪切带控制呈南西向展布的特征一致.
在双尖子山矿段,断裂表现为由北东向深断裂F1与北西向浅断裂F13、F14组成的共轭断裂.电性结果显示,F1断裂深度可延伸至地下2000 m,与深部高阻侵入岩连通,为深部成矿热液上涌提供良好通道,同时F1断裂在矿集区的剩余重力异常结果中(图 10a)显示为一条明显的北东向重力梯度带,推测其可能是双尖山子矿段内一条提供成矿物质的主要控矿断裂.
图 9为矿集区三维电性结构的高阻体(>8000 Ωm)分布特征.高阻体主要分布于地表以下,在A、B、C三个区域隆起,分别对应兴隆山矿段西北侧、西南侧以及双尖子山矿段西部(辉绿玢岩脉和流纹岩出露范围).根据侵入岩高阻、高磁的物性特征,并结合矿集区剩余磁异常结果(图 10b)表现为大面积高磁异常的特征,推断电性结果中的高阻体分布代表了侵入岩的空间展布,表现为侵入岩在A、B、C三个区域具有高度侵位的特征.
井ZK12-37钻遇到花岗斑岩的地球化学测年结果显示,侵入岩成岩年龄为早白垩世(Zhai et al., 2020),该结果约束了高阻隆起A区侵入岩的成岩年龄.同时电性结果中深部高阻区具有连通性,推测研究区下方高阻区应为早白垩世花岗侵入岩的响应,是双尖子山矿集区的重要成矿物质来源(匡永生等,2014).
双尖子山矿段内,C区对应地表出露的辉绿玢岩脉和流纹岩的范围,在0 m切片中(图 8c)呈现为中低阻特征,推测是由于第四纪覆盖引起的;在100 m以下深度切片(图 8d),具有中高阻电性分布特征,并与深部高阻侵入岩连通.结合周围断裂的展布,推测辉绿玢岩脉及流纹岩的形成主要由于深部岩浆热液通过北东向深断裂F1上涌至浅部后,受F1断裂和北西向浅断裂F13、F14控制冷却、沉淀,并在矿段西侧F1、F13和F14共轭断裂附近出露.该形成过程限制了双尖子山矿段侵入岩的空间展布,导致侵入岩在浅部仅分布于C区.
3.3 双尖子山矿集区成矿模式大兴安岭地区自古生代开始经历了多重构造体制的叠加作用(刘建明等,2004;Wu et al., 2011).随着古亚洲洋的闭合(约250 Ma),研究区从古亚洲洋构造域逐渐转化为古太平洋构造域(Robinson et al., 1999;Xiao et al., 2003);在晚侏罗至早白垩世,大兴安岭南段地区逐渐转变为陆内伸展的构造背景(Meng,2003;Zhang et al., 2010),如图 11a所示,下地壳镁铁质物质受到上涌软流圈物质的热侵蚀,发生部分熔融形成岩浆,并沿着深大断裂上升侵位,在上升过程携带了部分古老下地壳物质,形成了含矿岩浆(王丰翔等,2016;王祥东,2017).
本文研究结果揭示了含矿岩浆侵位至双尖子山矿集区2000 m以浅分布特征.在兴隆山矿段(图 11b),早白垩世侵入岩在矿段西北侧(A区)、西南侧(B区)高度侵位至二叠纪大石寨组下方,导致上覆大石寨组发生弯曲变形,在矿段中部-南部形成穹窿的构造特征,在穹窿构造形成过程中将在浅部形成断裂,成为重要的导矿或容矿空间.侵入岩在此就位后,含矿岩浆析出含矿热液流体,由于受到成矿环境改变、大气水的参与以及与地层之间的水岩反应等因素的影响,含矿流体在浅部已存在或伴生的北东向断裂F4、北西向断裂F5、北西向韧性剪切带F6以及近南北向断裂F2、F3中沉淀(王丰翔等,2016;张志强,2018;江彪等,2019),形成以受断裂控制的热液脉状矿化矿床.
在双尖子山矿段(图 11c),早白垩世侵入岩通过北东向深断裂F1在矿段西侧高度侵位,在地表形成辉绿玢岩脉.由于受F1断裂的主要控制,导致成矿热液流体主要析出于双尖子山矿段西侧,在断裂F1以及与北西向断裂F13、F14的交汇处富集沉淀,最终形成双尖子山矿段.
4 结论(1) 双尖子山矿集区控矿要素主要为:赋矿地层大石寨组、控矿构造共轭断裂及成矿物质来源早白垩世侵入体.上覆中-低阻主要为大石寨组的反映,分布于1000 m以浅,在兴隆山矿段南部以及双尖子山矿段东部,大石寨组地层厚度较厚,具有良好的赋矿空间;在兴隆山矿段,控矿断裂埋深在500 m以浅,呈北东向、北西向及近南北向分布;在双尖子矿段,断裂为北东向和北西向分布,其中北东向F1深断裂可能为双尖子山矿段的主要控矿断裂;矿集区深部存在大面积的高阻异常,为早白垩世侵入岩的反映.
(2) 基于精细电性结构的刻画建立了双尖子山矿集区2000 m以浅的成矿模式,具体表现为:早白垩世侵入岩在兴隆山矿段西北侧以及西南侧高度侵位,使得兴隆山矿段上覆大石寨组弯曲变形,这个过程活化并伴生了浅部的断裂,侵入岩就位后析出含矿热液在浅部存在的北东向断裂、北西向断裂以及断裂交汇处,受到成矿环境变化后沉淀,形成矿床;而在双尖子山矿段,通过北东向深断裂F1的主要控制,含矿流体在F1断裂以及其与北西向断裂的交汇处富集沉积,形成了矿体主要于西矿段分布的特征.
致谢 感谢审稿专家对文章提出的宝贵建议.
Broyden C G. 1970. The convergence of a class of double-rank minimization algorithms: 2. the new algorithm. IMA Journal of Applied Mathematics, 6(3): 222-231. DOI:10.1093/imamat/6.3.222 |
Chen L, Zhang D, Di Y J, et al. 2009. A study on the regional metallogenic regularity in the mid-southern segment of Da Hinggan Mountains, China. Contributions to Geology and Mineral Resources Research (in Chinese), 24(4): 267-271, 281. |
Commer M, Newman G A. 2008. New advances in three-dimensional controlled-source electromagnetic inversion. Geophysical Journal International, 172(2): 513-535. DOI:10.1111/j.1365-246X.2007.03663.x |
Du Q S, E A Q, Li Z H, et al. 2017. Basic geological features of polymetallic deposit concentration areas in Baiyinnuoer-Shuangjianzishan-Haobugao. Nonferrous Metals Science and Engineering (in Chinese), 8(2): 63-69. |
Fan W M, Guo F, Wang Y J, et al. 2003. Late Mesozoic calc-alkaline volcanism of post-orogenic extension in the northern Da Hinggan Mountains, northeastern China. Journal of Volcanology and Geothermal Research, 121(1-2): 115-135. DOI:10.1016/S0377-0273(02)00415-8 |
Fletcher R. 1970. A new approach to variable metric algorithms. The Computer Journal, 13(3): 317-322. DOI:10.1093/comjnl/13.3.317 |
George P L, Hecht F, Saltel E. 1990. Fully automatic mesh generator for 3D domains of any shape. IMPACT of Computing in Science and Engineering, 2(3): 187-218. DOI:10.1016/0899-8248(90)90012-Y |
Goldfarb D. 1970. A family of variable-metric methods derived by variational means. Mathematics of Computation, 24(109): 23-26. DOI:10.1090/S0025-5718-1970-0258249-6 |
Guo Y C, Chen R Y, Jia B, et al. 2017. Zircon U-Pb dating and geochemistry of the granite porphyry from the Shuangjianzishan silver polymetallic deposit in Inner Mongolia and tectonic implication. Geology and Exploration (in Chinese), 53(3): 495-507. |
Guo Z W, Hu L Y, Liu C M, et al. 2019. Application of the CSAMT method to Pb-Zn mineral deposits: A case study in Jianshui, China. Minerals, 9(12): 726. DOI:10.3390/min9120726 |
Hu X Y, Peng R H, Wu G J, et al. 2013. Mineral exploration using CSAMT data: Application to Longmen region metallogenic belt, Guangdong Province, China. Geophysics, 78(3): B111-B119. DOI:10.1190/geo2012-0115.1 |
Jiang B, Zhang T, Chen Y C, et al. 2019. The source of ore-forming fluid and mechanism of metal precipitation in the Shuangjianzishan Ag-polymetallic deposit in Chifeng, Inner Mongolia. Acta Geologica Sinica (in Chinese), 93(12): 3166-3182. |
Kuang Y S, Zheng G R, Lu M J, et al. 2014. Basic characteristics of Shuangjianzishan sliver polymetallic deposit in Chifeng City, Inner Mongolia. Mineral Deposits (in Chinese), 33(4): 847-856. |
Ledo J. 2005. 2-D versus 3-D magnetotelluric data interpretation. Surveys in Geophysics, 26(5): 511-543. DOI:10.1007/s10712-005-1757-8 |
Lin C H, Tan H D, Tong T. 2011. The possibility of obtaining nearby 3D resistivity structure from magnetotelluric 2D profile data using 3D inversion. Chinese Journal of Geophysics (in Chinese), 54(1): 245-256. DOI:10.3969/j.issn.0001-5733.2011.01.026 |
Liu C S. 2009. Three-dimension magnetotellurics adaptive edge finite-element numerical simulation based on unstructured meshes[Ph. D. thesis] (in Chinese). Changsha: Central South University.
|
Liu J M, Zhang R, Zhang Q Z. 2004. The regional metallogeny of Da Hinggan Ling, China. Earth Science Frontiers (in Chinese), 11(1): 269-277. |
Liu Y H. 2011. Research on 3-D controlled source electromagnetic method inversion using nonlinear conjugate gradients[Ph. D. thesis] (in Chinese). Changchun: Jilin University.
|
Liu Y J, Li W M, Feng Z Q, et al. 2017. A review of the Paleozoic tectonics in the eastern part of Central Asian Orogenic Belt. Gondwana Research, 43: 123-148. DOI:10.1016/j.gr.2016.03.013 |
Meng Q R. 2003. What drove late Mesozoic extension of the northern China-Mongolia tract?. Tectonophysics, 369(3-4): 155-174. DOI:10.1016/S0040-1951(03)00195-1 |
Ouyang H G, Li R H, Zhou Z H. 2016. The Jurassic mineralization of the Shuangjianzisha Ag-polymetallic deposit and its significance in prospecting: Evidence from geochronology. Acta Geologica Sinica (in Chinese), 90(8): 1835-1845. |
Plessix R E, Darnet M, Mulder W A. 2007. An approach for 3D multisource, multifrequency csem modeling. Geophysics, 72(5): SM177-SM184. DOI:10.1190/1.2744234 |
Robinson P T, Zhou M F, Hu X F, et al. 1999. Geochemical constraints on the origin of the Hegenshan Ophiolite, Inner Mongolia, China. Journal of Asian Earth Sciences, 17(4): 423-442. DOI:10.1016/S1367-9120(99)00016-4 |
Shanno D F. 1970. Conditioning of Quasi-Newton methods for function minimization. Mathematics of Computation, 24(111): 647-656. DOI:10.1090/S0025-5718-1970-0274029-X |
Shao J A, Zhang L Q, Mou B L. 1999. Magmatism in the Mesozoic extending orogenic process of Da Hinggan MTS. Earth Science Frontiers (in Chinese), 6(4): 339-346. |
Siripunvaraporn W, Egbert G, Uyeshima M. 2005. Interpretation of two-dimensional magnetotelluric profile data with three-dimensional inversion: synthetic examples. Geophysical Journal International, 160(3): 804-814. DOI:10.1111/j.1365-246X.2005.02527.x |
Thomas M D, Ford K L, Keating P. 2016. Review paper: Exploration geophysics for intrusion-hosted rare metals. Geophysical Prospecting, 64(5): 1275-1304. DOI:10.1111/1365-2478.12352 |
Wang F X, Sun H J, Pei R F, et al. 2016. The geologic features and genesis of Shuangjianzishan silver-polymetallic deposit, Balinzuo Qi, Inner Mongolia. Geological Review (in Chinese), 62(5): 1241-1256. |
Wang X D. 2017. Magmatism and mineralization of Ag-Pb-Zn polymetallic deposits in the Lindong district, Inner Mongolia[Ph. D. thesis] (in Chinese). Wuhan: China University of Geosciences.
|
Windley B F, Alexeiev D, Xiao W J, et al. 2007. Tectonic models for accretion of the Central Asian Orogenic Belt. Journal of the Geological Society, 164: 31-47. DOI:10.1144/0016-76492006-022 |
Wu F Y, Han R H, Yang J H, et al. 2007. Initial constraints on the timing of granitic magmatism in North Korea using U-Pb zircon geochronology. Chemical Geology, 238(3-4): 232-248. DOI:10.1016/j.chemgeo.2006.11.012 |
Wu G B, Liu J M, Zeng Q D, et al. 2014. Occurrences of silver in the Shuangjianzishan Pb-Zn-Ag deposit and its implications for mineral processing. Earth Science Frontiers (in Chinese), 21(5): 105-115. |
Wu H Y, Zhang L C, Wan B, et al. 2011. Geochronological and geochemical constraints on Aolunhua porphyry Mo-Cu deposit, northeast China, and its tectonic significance. Ore Geology Review, 43(1): 78-91. DOI:10.1016/j.oregeorev.2011.07.007 |
Wu X P, Liu Y, Wang W. 2015. 3D resistivity inversion incorporating topography based on unstructured meshes. Chinese Journal of Geophysics (in Chinese), 58(8): 2706-2717. DOI:10.6038/cjg20150808 |
Xiao W J, Windley B F, Hao J, et al. 2003. Accretion leading to collision and the Permian Solonker suture, Inner Mongolia, China: Termination of the central Asian orogenic belt. Tectonics, 22(6): 1069. |
Xu B, Zhao P, Bao Q Z, et al. 2013. Preliminary study on the pre-Mesozoic tectonic unit division of the Xing-Meng Orogenic Belt (XMOB). Acta Petrologica Sinica (in Chinese), 30(7): 1841-1857. |
Yu C M. 1998. The application of CSAMT method in looking for hidden gold mine. Acta Geophysica Sinica (in Chinese), 41(1): 133-138. |
Zhai D G, Williams-Jones A E, Liu J J, et al. 2020. The genesis of the giant Shuangjianzishan epithermal Ag-Pb-Zn Deposit, Inner Mongolia, northeastern China. Economic Geology, 115(1): 101-128. DOI:10.5382/econgeo.4695 |
Zhang D Q, Jiang B, Wang D H, et al. 2015. A summary of resources characteristic and metallogenic regularity of silver deposits in China. Acta Geologica Sinica (in Chinese), 89(6): 1008-1025. |
Zhang J H, Gao S, Ge W C, et al. 2010. Geochronology of the Mesozoic volcanic rocks in the Great Xing'an Range, northeastern China: Implications for subduction-induced delamination. Chemical Geology, 276(3-4): 144-165. DOI:10.1016/j.chemgeo.2010.05.013 |
Zhang J M, Zeng Z F, Zhao X Y, et al. 2020. Deep Mineral Exploration of the Jinchuan Cu-Ni Sulfide Deposit Based on Aeromagnetic, Gravity, and CSAMT Methods. Minerals, 10(2): 168. DOI:10.3390/min10020168 |
Zhang M, Zhai Y S, Shen C L, et al. 2011. Metallogenic system of copper polymetallic deposits in the middle-southern part of Da Hinggan Mountains, China. Geoscience (in Chinese), 25(5): 819-831. |
Zhang Z Q. 2018. Geological and geochemical characteristics and metallogenetic mechanism of the Shuangjianzishan Ag-polymetallic deposit, Inner Mongolia[Master's thesis] (in Chinese). Beijing: China University of Geosciences (Beijing).
|
Zhao Z, Zhou X P, Guo N X, et al. 2018. Superimposed W and Ag-Pb-Zn (-Cu-Au) mineralization and deep prospecting: Insight from a geophysical investigation of the Yinkeng orefield, South China. Ore Geology Reviews, 93: 404-412. DOI:10.1016/j.oregeorev.2018.01.017 |
Zorin Y A. 1999. Geodynamics of the western part of the Mongolia-Okhotsk collisional belt, Trans-Baikal region (Russia) and Mongolia. Tectonophysics, 306(1): 33-56. DOI:10.1016/S0040-1951(99)00042-6 |
陈良, 张达, 狄永军, 等. 2009. 大兴安岭中南段区域成矿规律初步研究. 地质找矿论丛, 24(4): 267-271, 281. |
杜青松, 鄂阿强, 李志华, 等. 2017. 白音诺尔-双尖子山-浩布高多金属矿集区基本地质特征. 有色金属科学与工程, 8(2): 63-69. |
顾玉超, 陈仁义, 贾斌, 等. 2017. 内蒙古双尖子山银多金属矿床花岗斑岩年代学、地球化学特征及构造意义. 地质与勘探, 53(3): 495-507. |
何继善. 1990. 可控源音频大地电磁法. 长沙: 中南工业大学出版社.
|
江彪, 张通, 陈毓川, 等. 2019. 内蒙古赤峰双尖子山银多金属矿床成矿流体来源及金属沉淀机制探讨. 地质学报, 93(12): 3166-3182. DOI:10.3969/j.issn.0001-5717.2019.12.011 |
匡永生, 郑广瑞, 卢民杰, 等. 2014. 内蒙古赤峰市双尖子山银多金属矿床的基本特征. 矿床地质, 33(4): 847-856. DOI:10.3969/j.issn.0258-7106.2014.04.014 |
林昌洪, 谭捍东, 佟拓. 2011. 利用大地电磁三维反演方法获得二维剖面附近三维电阻率结构的可行性. 地球物理学报, 54(1): 245-256. DOI:10.3969/j.issn.0001-5733.2011.01.026 |
刘长生. 2009. 基于非结构化网格的三维大地电磁自适应矢量有限元数值模拟[博士论文]. 长沙: 中南工业大学.
|
刘建明, 张锐, 张庆洲. 2004. 大兴安岭地区的区域成矿特征. 地学前缘, 11(1): 269-277. DOI:10.3321/j.issn:1005-2321.2004.01.024 |
刘云鹤. 2011. 三维可控源电磁法非线性共轭梯度反演研究[博士论文]. 长春: 吉林大学.
|
欧阳荷根, 李睿华, 周振华. 2016. 内蒙古双尖子山银多金属矿床侏罗纪成矿的年代学证据及其找矿意义. 地质学报, 90(8): 1835-1845. DOI:10.3969/j.issn.0001-5717.2016.08.013 |
邵济安, 张履桥, 牟保磊. 1999. 大兴安岭中生代伸展造山过程中的岩浆作用. 地学前缘, 6(4): 339-346. DOI:10.3321/j.issn:1005-2321.1999.04.017 |
石昆法. 1999. 可控源音频大地电磁法理论与应用. 北京: 科学出版社.
|
王丰翔, 孙红军, 裴荣富, 等. 2016. 巴林左旗双尖子山银多金属矿床基本地质特征及成因机制. 地质评论, 62(5): 1241-1256. |
王祥东. 2017. 内蒙古林东地区银铅锌多金属矿床成岩成矿作用[博士论文]. 武汉: 中国地质大学(武汉).
|
吴冠斌, 刘建明, 曾庆栋, 等. 2014. 内蒙古双尖子山铅锌银矿床银的赋存状态及其指示意义. 地学前缘, 21(5): 105-115. |
吴小平, 刘洋, 王威. 2015. 基于非结构网格的电阻率三维带地形反演. 地球物理学报, 58(8): 2706-2717. DOI:10.6038/cjg20150808 |
徐备, 赵盼, 鲍庆中, 等. 2013. 兴蒙造山带前中生代构造单元划分初探. 岩石学报, 30(7): 1841-1857. |
于昌明. 1998. CSAMT方法在寻找隐伏金矿中的应用. 地球物理学报, 41(1): 133-138. DOI:10.3321/j.issn:0001-5733.1998.01.015 |
张大权, 江彪, 王登红, 等. 2015. 中国银矿的资源特征及成矿规律概要. 地质学报, 89(6): 1008-1025. DOI:10.3969/j.issn.0001-5717.2015.06.002 |
张梅, 翟裕生, 沈存利, 等. 2011. 大兴安岭中南段铜多金属矿床成矿系统. 现代地质, 25(5): 819-831. DOI:10.3969/j.issn.1000-8527.2011.05.001 |
张志强. 2018. 内蒙古双尖子山银多金属矿床地质地球化学特征与成矿机制[硕士论文]. 北京: 中国地质大学(北京).
|