华南陆块东北部构造特征、成矿作用与动力学机制
——来自大地电磁测深的认识

2022-06-07 03:27董基恩叶高峰魏文博金胜李玉龙董小军魏有宁
地质论评 2022年3期
关键词:电性扬子华南

董基恩,叶高峰,魏文博,金胜,李玉龙,董小军,魏有宁

1)中国地质大学(北京)地球物理与信息技术学院,北京,100083;2)青海省地质调查局,西宁,810000

内容提要:通常认为,华南陆块在新元古代由扬子克拉通和华夏地块沿江南造山带碰撞拼合而成,之后经历了陆内造山、洋壳俯冲等多期岩浆—构造活动。但因华南陆块所处的特殊地质构造环境,就目前华南陆块各块体之间的接触关系、江南造山带深部构造特征及区域动力学意义等诸多地质问题争议颇多。本文依托地质调查项目和“深部探测技术与实验研究(SinoProbe)”项目完成的8d和12g两条大地电磁测深剖面,经过精细的数据处理,使用非线性共轭梯度(NLCG)算法对TE+TM模式数据联合反演得到华南陆块东北部岩石圈尺度的二维电性剖面,并用ModEM三维反演代码对全阻抗张量数据反演获得了三维电性模型。对研究区内扬子地块东部、江南造山带以及华夏地块进行电性结构研究,发现研究区内的江南造山带西南段存在扬子地块和华夏地块碰撞镶嵌的构造表现,扬子地块已越过江南造山带,在江山—绍兴断裂位置与华夏地块挤压,形成江绍断裂等逆冲型深大断裂,从电性结构推测其现今仍然为活动断裂,但东北段块体之间的接触关系被上侵的地幔物质破坏,江绍、赣东北等断裂的深部结构已被剧烈改造,推测这种深部成矿热物质上涌是形成赣东北以金银矿种为主的岩浆热液型矿床的深部动力原因;研究区东部华夏地块电性特征为高阻的上地壳以及被岩浆底侵而破坏的中下地壳,发育其中的屯溪—鹰潭—安远和上虞—大浦—政和断裂切割深度超过了50 km,为深大断裂构造。结合前人地表侵入岩填图结果,认为由于中生代以来古太平洋板块向欧亚大陆的俯冲,华南陆块东部上地壳被严重破坏,从电性特征推断可能发生过大面积的板片重融,华夏地块东北部地壳相较于西南部厚度明显更小、后期改造严重,受到的太平洋板块俯冲导致的热扰动更剧烈。

华南陆块是由扬子地块和华夏地块两个一级构造单元组成,中间沿江南造山带拼接在一起(舒良树,2012)。自显生宙加里东运动以来,在各类地球动力学背景下的不同地质作用和构造演化造就了华南当前复杂的地质环境属性和特殊的地貌风格以及资源分布情况,所以近半个世纪以来一直备受地质工作者的关注,是全球性的研究热点。古生代晚期至中生代早期,华南陆块与北侧的华北地块碰撞挤压,形成秦岭—大别造山带;中生代以来,华南陆块受到东侧太平洋板块向西的俯冲挤压,爆发了多期次的岩浆活动(毛景文等,2004);而其西南,中、新生代受到印支陆块的东北向的挤压,通过龙门山构造带与青藏地块相连。虽经过几代人的努力,其地表地质调查以及实际踏勘等工作比较充分,但由于深部资料的缺乏,使得目前对华南陆块的研究在华南陆内造山机制、华南岩石圈壳幔结构以及地球深部动力学等方面争论众多。所以在华南地区开展系统、密集的岩石圈尺度的地球物理勘探也是适应当前科研要求的重要手段。

大地电磁测深是目前主要的地球深部探测方法之一,是从导电性的角度研究地壳和上地幔结构不可缺少的方法(魏文博,2002;魏文博等,2003)。所以针对华南大陆复杂的地质环境和深部构造,在华南地区开展大地电磁测深工作,获得华南地壳的电性模型,这将对华南地区壳/幔结构及其动力学背景研究提供更具有说服力的资料,并且对于探讨地质规律、解析矿产成因、评估地质环境和地质灾害形成机制,以及研究大陆演化的动力学过程等都具有重要的意义(刘国兴等,2012)。

本文利用两条穿越华南陆块东北部地区的大地电磁(MT)剖面,获得了华南陆块东北部地区岩石圈尺度的二维及三维电阻率模型,对研究区内扬子地块、华夏地块及江南造山带的电性特征进行了地球物理解释,并根据电性特征就研究区内扬子克拉通和华夏地块接触关系、赣东北成矿带深部成矿机制及地球动力学意义等问题进行了讨论。

1 地质背景

华南陆块在西太平洋板块西向俯冲与印度洋板块北向差异运动的夹持等长期复杂的板块构造作用下,形成了目前独特的大陆构造。现今华南整体陆壳基本由扬子与华夏两个地块和其上的显生宙盖层构成,两地块之间的中间融合带就是所谓“江南造山带”。扬子地块是指华南大陆中具有前南华纪基底的扬子古微板块和华夏古微板块的西部的统称,地球化学测年研究发现扬子地块深部可能普遍存在一个新太古代结晶基底,并有古元古代岩浆活动遗迹,而且扬子地块中发现的非常古老的地质体崆岭杂岩是形成于2.90~2.95 Ga (Qiu Yumin et al.,2000;Zhang Shaobing and Zheng Yongfei,2013),说明扬子广泛残存的是早前寒武纪结晶基底,现已基本完成相对稳定克拉通陆块的演化。重磁研究发现扬子地块主要为正磁异常以及负的布格重力异常(张季生和洪大卫,2002)。结合重磁电震地球物理方法研究成果,扬子陆块总体呈现高速高密度高阻中高温的地球物理特性。

自从Grabau(1924)提出“华夏古陆”的认识以来,华夏地块是否真的存在,争议巨大。但是后来在浙南—闽北地区开展的早元古代同位素年代学研究使得“华夏古陆”又被重新提出(于津海等,2006)。所以目前对于华夏地块的分布以及规模还是存在不同认识(李献华,1998;邓平等,2002;舒良树,2006)。目前大家所公认的华夏地块位于中国的东南部,是华南大陆的组成部分,呈北东向分布。华夏地块主要由各个时代变质岩及中生代花岗岩为主的岩浆岩组成,其火成岩主要是在加里东以后形成的,但是在闽北、浙南等地区有古元古代和新元古代变质火成岩分布 (甘晓春等,1993;甘晓春等,1995;李献华,1998)。于津海等(2006)认为华夏地块目前的岩石主要是在新元古代形成,其中包含古老的“华夏古陆”,但被后期的构造热事件强烈改造。航磁研究发现华夏板块以负磁异常为主,间有局部正异常,这与酸性岩浆活动有关,呈现相对年轻的地壳属性(童迎世和童琼,2012)。饶家荣等(2012)结合各类地球物理资料认为华夏陆块为造山带岩石圈,总体呈低速低密度低阻高温经强烈改造的相对活动型(热)陆块。

江南造山带的地质研究工作历史悠久。李四光、黄汲清等老一批地质学家在上世纪初期先后到南岭地区进行过地质矿产研究工作,初步确定了研究区内分布的矿产种类和特点以及附近主要的构造轮廓。20世纪80年代初,水涛(1987)对江山—绍兴断裂带及其两侧开展了长期的野外地质调查,认为扬子古陆与华夏古陆之间曾为浩海分隔并发生过碰撞,由此提出两个块体碰撞对接的构想。随后杨明桂和梅文勇(1997)在该区域进行了比较系统的地质构造及成矿作用研究,确认了钦州湾至杭州湾为扬子与华夏这两个古板块的缝合带,称为钦杭结合带(江南造山带)。总的来说结合带南西起自广西钦洲湾,经湘东和赣中,往北东延伸至浙江杭州湾,总体呈NE 向反S 状弧形展布,全长约2000 km,宽100~300 km(徐德明等,2013),考虑到该结合带也是最重要的构造岩浆成矿带,并正式命之为“钦杭成矿带”。钦杭成矿带也具有分段性,北段的赣东北成矿带的优势矿种为铜、铁和贵金属,西南段的南岭成矿带是重要的钨锡多金属产地(毛景文等,2004)。

针对华南陆块特殊的地质背景及构造环境,一大批地球物理工作者就江南造山带边界问题、华南陆块东部岩石圈减薄等地质焦点问题,开展了一系列地球物理工作,并取得一些重要的成果。在地震研究方面,Zhang Zhongjie等(2013)通过清远连县—惠东港口地震剖面,建立了华南地壳上地幔P波速度结构模型,发现P波速度在上地壳横向变化大,下地壳横向变化小,并认为吴川—四会断裂是扬子与华夏的地壳分界线。江国明等(2014)通过长江中下游地区速度异常分析,得到上地幔“两高夹一低”的速度异常体,认为高速体代表拆沉的岩石圈,低速为上涌的软流圈物质,说明该地区可能发生拆沉。Dong Shuwen等 (2020)通过华南深反射地震剖面认为华夏东部发生了明显的岩浆底侵作用以及软流圈上涌,底侵为地壳增加了大量的幔源物质,弥补了地壳减薄的影响。在大地电磁测深研究方面,刘国兴等(2012)通过华南沿海地区MT剖面认为该地区深部地壳总体上呈东薄西厚的特征,西段岩石圈增厚可能由于陆内造山作用,中段岩石圈减薄可能由于软流圈物质上涌,东段岩石圈的增厚可能是太平洋俯冲导致洋壳残留物叠置形成双层岩石圈,而太平洋俯冲导致的深部热扰动作用可能引起大规模软流圈上涌。韩松(2017)通过进荣—柘贤的MT剖面研究,认为华夏沿海地区并非为岩石圈的大范围减薄,而是以局部软流圈上涌为途径来影响和改造上覆岩石圈结构。Hu Haoyuan等(2020)通过大地电磁测深建立了华南东部岩石圈3-D电性结构模型,认为赣东北成矿带从地幔延伸至地壳的低电阻率异常为成矿物质的岩浆通道,而南岭成矿带地幔和地壳的低电阻率异常通过较小的低电阻率通道连接,可能是软流圈物质上升的通道,并且上升的软流圈物质被高电阻率的下地壳阻塞,地幔上升流的热量和成矿元素可能触发了上地壳的部分熔融。在综合地球物理探测方面,严加永等(2019)使用大尺度卫星重磁资料,结合区域电性和地震资料的综合研究,认为钦杭结合带是江南造山带的南部边界,结合带南界为宁波—金华—上饶南—赣州北—郴州—北海东,北界为上海—湖州—鹰潭—临川—萍乡—衡阳—永州—桂林—梧州—钦州。

2 大地电磁数据采集及处理分析

2.1 数据采集与处理

研究区MT数据由中国地质大学(北京)和吉林大学分别采集于2016年和2010~2012年。主要由穿过江南造山带的短剖面及东南沿海地区的剖面组成的阵列。此次研究中,8d和12g剖面分别选取大地电磁测深点55个和43个(图1红色和蓝色圆点)。野外数据采集时,取北方向为x轴正方向,东方向为y轴正方向,垂直向下为z轴正方向。在所有测点使用加拿大凤凰公司生产的MTU-5A仪器进行了宽频大地电磁测深观测,在采集时间大于40 h的情况下可以采集到频率320~0.0002 Hz范围内的电磁场信号,根据趋肤深度的表达式,本区域宽频仪器的最大勘探深度通常在70 km左右。为了采集更长周期的大地电磁信号以获得研究区岩石圈尺度的电性结构特征,根据宽频数据质量及研究需求在其中16个测点进行了长周期观测(图1蓝色圆点),采集仪器为乌克兰科学院空间研究所研制的Lemi-417m,每个测点的采集时间在7 d以上。

图1 华南陆块东北部主要构造及MT测点分布简图Fig.1 Location of MT stations with geological division of the survey area① 绍兴—江山断裂;② 郴州—临武断裂;③ 九江—石台断裂;④ 吉首断裂;⑤ 屯溪—鹰潭—安远断裂;⑥ 上虞—政和—大浦断裂;⑦ 赣江断裂;⑧ 郯庐断裂;⑨赣东北断裂① Shaoxing—Jiangshan fault;② Chenzhou—Linwu fault;③ Jiujiang—Shitai fault;④ Jishou fault;⑤ Tunxi—Yingtan—Anyuan fault;⑥ Shangyu—Zhenghe—Dapu fault;⑦ Ganjiang fault;⑧ Tanlu fault;⑨ Northeastern Jiangxi fault

野外采集得到测点电磁场时间序列后,使用SSMT2000软件对宽频数据进行处理,并进行精细的功率谱密度挑选。针对华南多丘陵、多雨、干扰源丰富等特点,2016年进行宽频数据采集时,在湖南汉寿县以及湖北孝感市布设两台远参考站,后期处理中以远参考点磁道为参考道以及各个点互参考技术(Gamble et al.,1979)优化数据质量,对于2010~2012年SinoProbe专项东南沿海MT数据采用点间互参考进行重处理,使得数据质量得到了很大改善。由于采集时间长,宽频数据的最长周期达到10000 s。长周期数据处理采用Varentsov教授的代码(Varentsov et al.,2003),处理得到阻抗张量,并选取100 s至1000 s之间某个频点将宽频数据和长周期数据进行拼接,这样可以得到该测点完整的测深曲线,经拼接后数据质量好的点其最长周期可达到40000 s左右(图2中间幅)。华南经济发达,矿产资源丰富,因此人文噪音较强。在天然场的“死频带”10~0.1 Hz通常数据质量较差,某些测点即便进行了远参考处理还未能有效改善“死频带”数据质量,对这些受明显干扰的频点进行了删除,未参与二维、三维反演。

图2 华南陆块主要构造单元典型测点的视电阻率和相位曲线图Fig.2 Apparent resistivity and phase curves of three typical stations in main geological blocks of the South China Block)

在不同的地质体内其测深曲线具有不同特征,对图1中所示的0415、0804和3308 MT测点绘制了测深曲线图(图2)。扬子地块内测点(0415)的测深曲线表现为低阻—高阻—低阻—较高阻的曲线变化特征,表明扬子地块的基底可能为高阻;位于江南造山带内的测点(0804),XY和YX模式视电阻率曲线变化具有明显差异形态,为地质构造复杂情况下的大地电磁测深视电阻率曲线;华夏地块内的测点(3308)曲线类型可以看出华夏地块上地壳为稳定高阻特征,中下地壳电阻率明显降低。

2.2 维性、构造走向分析

在开展反演工作之前需进行数据的维性分析,主要是进行二维偏离度分析。本文选择了受电磁场畸变效应较小的Bahr二维偏离度对8d、12g剖面进行维性分析(Bahr,1988)。当二维偏离度值小于0.3时可以认为剖面满足二维假设,当二维偏离度值大于0.3时,说明地下介质存在三维效应。如图3所示,位于西南的 8d剖面二维性优于12g剖面,说明研究区内的江南造山带由西南向东北方向其地下构造逐渐复杂,逐渐趋于三维特性。在8d剖面中仅在东南沿海地区深部存在小范围的三维性特征,近似满足二维假设。而12g剖面二维性较差,在华夏地块内部存在比较明显的三维属性,因此在后期反演中也进行了三维反演,旨在进行二、三维模型对比研究,获得尽可能接近实际情况的电阻率模型。

图3 华南8d、12g剖面Bahr二维偏离度拟断面图Fig.3 Pseudo-sections of Bahr skew of profile 8d and 12g in the northeastern South China Block

相位张量分析具有不受局部电场畸变影响的优点(Caldwell et al.,2004),也被用于判断维性和构造走向。如图4所示,对0.1 s、10 s、100 s和1000 s的4个周期进行了相位张量分析,并绘制了各周期的相位张量椭圆。可以看出,研究区域内绝大部分测点反映浅部信息的高频率(0.1 s)的二维偏离角β都小于3°,说明MT数据在高频段的二维性较好(Booker,2014),部分椭圆接近退化为圆,反映了这些块体的浅表沉积层具有较好的一维特性,少数几个β>5°的孤立测点,可能因为数据质量较差所导致的。但在中频段(10 s和100 s),研究区域内部分MT测点数据表现出了一定的三维特性(β>3°),位于块体边界带或者深大断裂带附近大部分测点的β>5°的情况,体现出研究区深部的构造情况复杂的特征,存在三维效应,这与Bahr二维偏离度的分析结果基本一致。

图4 华南陆块东北部四个周期的相位张量分析图Fig.4 The phase tensor analysis results at four periods in the northeastern South China Block

在维性分析结果确认可以进行二维反演后,需分析剖面所在区域的构造走向,只有将MT数据旋转至构造走向方向上,才能将MT数据分解为两组相互独立的极化模式:TE和TM模式。从相位张量分析获得的椭圆结合区域地质资料可知,研究区构造走向主要为北东向。为了确定具体的构造走向角,对MT数据进行了阻抗张量分解。G—B分解是目前MT张量畸变校正应用最为广泛的方法之一,该方法将张量阻抗分解为一个区域二维张量阻抗和一个局部畸变矩阵,获得两个相互正交的电性主轴角,据此所确定的构造走向角有90°的不确定性(Groom and Bailey,1989),需要结合磁感应矢量分析或区域构造图才能确定。图5为单点全频段张量阻抗分解结果和多点分频段(0.01~0.1 s,0.1~1 s,10~100 s和100~1000 s)张量阻抗分解结果图,结果显示研究区有较明显的电性主轴方向,结合区域构造走向,可以确定剖面8d的优势主轴方向为NE35°,12g的优势主轴方向为NE40°,因此将8d和12g剖面中的所有数据分别旋转至NE35°和NE40°坐标系下进行二维反演研究。

图5 华南陆块东北部 8d和12g剖面测点构造走向分析统计玫瑰图Fig.5 The analysis of statistical roses of profiles 8d and 12g in the northeastern South China Block

3 反演

3.1 二维反演

此次二维反演采用多模式,多参数相结合的思路进行多次二维反演试算,为了让更多的数据参与反演并约束地电模型,此次反演使用以非线性共轭梯度(NLCG)为计算内核的WinGLink软件(Rodi and Mackie,2001),采用TM+TE联合反演模式,初始模型为100 Ω·m的均匀半空间,TM模式相位误差基数5%,TM模式视电阻率误差基数10%,考虑到TE模式受三维畸变效应的影响比较严重,而且在后期的反演中发现TE模式参与反演的模型容易出现大面积的高阻异常,并且数据难以拟合,因此将TE模式视电阻率误差基数增大至50%,相位误差基数设置为10%,采用不同的正则化因子τ值进行多次反演,找到模型光滑度和均方根(Root Mean Square,简称RMS)误差的最佳平衡。最终选择τ=10的正则化因子,经过200次迭代,得到8d剖面RMS误差为2.38,12g剖面RMS误差为2.37。图6为反演剖面TE和TM模式实测数据视电阻率与相位以及最终二维模型响应数据视电阻率与相位拟断面图,对比可以看出TM模式的实测数据和响应数据拟合良好,TE 模式视电阻率稍有差别,是视电阻率误差基数设置较大导致的,但是总体来说拟合情况也较好,证明最终反演模型是可靠的。

图6 华南陆块东北部TM和TE模式视电阻率与相位实测数据和模型响应拟断面图Fig.6 The pseudo-sections of observed and calculated apparent resistivity and phases of TM and TE modes data in the northeastern South China Block

3.2 三维反演

由于当下计算机技术的发展已经能够满足大地电磁三维反演的硬件需求,大地电磁测深三维反演技术已经日渐成熟,当下主流的三维反演计算代码为基于数据空间的OCCAM算法的WS3DINVMT反演程序(Siripunvaraporn et al.,2005)以及基于非线性共轭梯度算法的 ModEM(Kelbert et al.,2014)。

图7 华南东部大地电磁测深三维反演各测点拟合误差分布图Fig.7 Distribution of RMS misfits of three-dimensional MT inversion of eastern South China Block

4 研究区电性结构特征

对二维及三维反演获得的电阻率模型的对比发现,12g剖面二维和三维反演的结果较为一致;虽然8g剖面两者之间差别较大,但大的结构是一致的。通常而言,三维反演会获得较为光滑的模型,对边界的分辨能力不如二维反演,因此对边界的解释以二维反演为主。结合前人研究成果(舒良树等,2004;舒良树,2006;张岳桥等,2009;张国伟等,2013;Zhang Letian et al.,2015;He Chuansong and Santosh,2016),最终得到如图8及图9所示华南陆块东北部地质—地球物理解释模型。其中8d剖面横跨扬子地块、江南造山带及华夏地块构造单元,并经过了九江—石台断裂、江山—绍兴断裂、赣江断裂带、屯溪—鹰潭—安远断裂以及上虞—政和—大浦断裂等大型断裂构造,剖面总长约800 km。综述8d剖面电性特征可概括为高阻的地壳被具有明显的电性梯度带特征的分界断裂所切割,为明显的“横向分块”的特征,而被切割的上地壳内部为较稳定高阻响应,电阻率最高达到10000 Ω·m,为较稳定的微陆块。由于剖面位于扬子地块的测点较少,剖面西段主要反映了扬子地块与江南造山带接触部位的电性特征,表现为中上地壳的高阻块体明显被挤压破坏(如图8中R1所示),可以发现向下延伸的裂隙式低阻异常切割了高阻的上地壳,该区域为江南造山带西边界断裂九江—石台断裂的电性特征,该断裂为高角度东倾逆冲断裂,切割深度达到了50 km左右,断裂与两侧电性差异较大。在九江—石台断裂与江山—绍兴断裂之间为江南造山带,带内电性特征较复杂,主要表现为被挤压破坏的高阻上地壳,而中下地壳中存在一高阻体(如图8中R2所示),深度大于100 km,推测该高阻体应为扬子地块结晶基底,该基底已经越过江南造山带,在江绍断裂附近与华夏地块碰撞挤压,形成了江南造山带东边界断裂江山—绍兴断裂,模型反映该处下地壳电阻率明显减小,这样的电性转换特征也说明了该处地壳及深部结晶基底已经被破坏,为典型的造山型地壳类型。江山—绍兴断裂表现为深大断裂构造特征,该断裂角度较缓,向东南倾斜并延伸至钦杭构造带深部,切割深度达到了近100 km,为一逆冲断裂(余心起等,2005;Zhang Letian et al.,2015),野外地质调查也认为江绍断裂带应是一条岩石圈断裂(杨明桂等,2009)。在断裂浅部存在深度约为几千米的层状低阻异常,实地地质调查研究认为可能是古生代沉积地层(余心起等,2005),并且以角度不整合覆盖于结晶基底之上(舒良树等,2004)。电性模型中华夏地块为普遍的高阻异常(如图8中R3、R4、R5所示),但是被8d剖面主要经过的赣江断裂带、屯溪—鹰潭—安远断裂及上虞—大浦—政和等断裂切割,各个被切割开的高阻块体内部并没有遭受到破坏,为稳定的高阻特性,电阻率值达到了10000 Ω·m。剖面中赣江断裂带规模较小,只发育在上地壳,切割深度50 km左右,前人综合研究结果已揭示该断裂为东倾正断层性质(梁兴和吴根耀,2006);屯溪—鹰潭—安远断裂发育在华夏地块内部,该断裂电性特征较复杂,在浅层断裂角度较缓,而深部倾斜角度较大,并变为西倾断裂,该断裂深部有热物质活动并上涌;上虞—大浦—政和断裂的电性特征显示为上地壳断裂,断裂仅仅切穿上地壳,这个与Zhang Letian等(2015)研究相符,同时前人的电阻率模型推测上虞—大浦—政和断裂可能在深部存在一个西倾的隐伏断裂 (韩松,2017),这与本次研究结果相符合。本研究区莫霍面埋深根据前人的研究(熊小松等,2009;Zhou Longquan et al.,2012)在电性剖面中标出,模型指示莫霍面深度与高阻上地壳底界面基本相符。

图8 华南陆块东北部8d剖面电性结构模型及构造解释图Fig.8 Resistivity models from 2-D and 3-D inversions of the 8d profile in the northeastern South China Block and its geological interpretation

图9 华南陆块东北部12g剖面电性结构模型及构造解释图Fig.9 Resistivity models from 2-D and 3-D inversions of the 12g profile in the northeastern South China Block and its geological interpretation

12g剖面位于研究区最北边,横跨扬子地块、江南造山带及华夏地块三大地质构造单元,经过了九江—石台断裂、江山—绍兴断裂、赣东北断裂系及上虞—政和—大浦断裂等大型构造断裂,剖面总长度约580 km。通过12g与8d剖面的对比研究发现江南造山带往东北,在较短距离内电性结构发生了剧烈变化,我们认为这种短距离的电性结构突变指示了华南陆块已经历多期次构造改造,深部构造复杂多变的特征。

在12g剖面中已没有如8d剖面中那样高阻的扬子地块(如图8中R1所示),而江南造山带内仍存在规模很小的高阻块体R2(如图9中R2所示),相比于8d剖面R2(如图8中R2所示),其带内发育的高阻体R2规模明显减小,深部延伸深度约50 km,而且高阻体为东南倾斜的特征,与此同时,该处发育的江山—绍兴断裂切割深度达50 km,为逆冲断裂构造。江南造山带及赣东北断裂带中下地壳主要为低阻异常,其低阻物质甚至已经上侵至上地壳,我们推测后期地幔物质上涌破坏了上地壳原有的高阻属性,使得上地壳发生较大规模的热物质底侵改造作用,而且越靠近江南造山带内底侵改造作用越明显,这说明挤压碰撞而导致的深部热扰动主要集中在结合带内。华夏地块整体表现为稳定的高阻块体特征,高阻体延伸深度达到了100 km,局部发育了上涌的低阻异常(如图9中C2所示),在三维模型中具有更明显的地幔物质上涌特征,但是总体来看东南部的华夏地块为较稳定的块体。上虞—政和—大浦断裂依旧是两条倾向不一致的断裂组成,并没有切穿上地壳,只是发育在上地壳的小型断裂。通过两条剖面对比发现屯溪—鹰潭—安远断裂只发育在北纬29°以南,在12g剖面中没有发现该断裂。

5 研究区构造涵义讨论

5.1 扬子、华夏地块接触关系

扬子和华夏地块接触关系的研究是目前华南大陆地质研究的主要争论所在。截止目前,前人已经从地质、地球物理、地球化学等各个方面开展相关研究,并取得一些重要认识。地质研究认为现今组成华南大陆的扬子地块和华夏地块在中、新元古代并不是统一的块体,而是由多个微陆块体拼合而成,主要为前寒武纪的陆壳块体,而且分布范围较大(王鸿祯和莫宣学,1996;舒良树,2006;于津海等,2006),指示板块拼接的赣东北蛇绿混杂岩带是华南陆块中唯一能确定的蛇绿混杂岩带,主要年代为距今900 Ma左右(张旗等,2001;Wang Yuejun et al.,2013b),并在武夷山地区普遍存在距今约900 Ma年的碎屑锆石(Wang Pengcheng et al.,2012;Wang Yuejun et al.,2012,2013a),综上研究认为当时依然独立的扬子和华夏古板块在距今约900 Ma年时完成内部拼合,形成各自统一的地块(张国伟等,2013)。后期的研究发现存在于华南大陆中板溪群组的同位素化学测年及分布于雪峰山、九岭等地区的新生代超基性花岗岩锆石U-Pb定年距今约850~820 Ma(张国伟等,2013),并结合区域角度不整合等构造表现认为在850~820 Ma(晋宁期)之前形成完整的扬子古陆和华夏古陆碰撞拼合,并形成新元古代的江南造山带,标志着华南大陆统一基底的形成。但需要指出的是目前分布于赣东北蛇绿混杂岩带是扬子和华夏两个微陆块在新元古代的缝合线。再后来,在Rodinia超大陆裂解的构造条件影响下,华南大陆从南华纪开始进入拉张裂谷环境,形成了政和—大浦断裂及岩浆活动(舒良树,2012),最终成为陆内海盆所分割的陆内两个地块(张国伟等,2013)。根据沉积古地理环境研究以及江南造山带周缘的碎屑锆石研究(Wang Yuejun et al.,2013b),认为在早古生代华夏地块与扬子地块之间并没有洋壳存在,对于“古华南洋”还待商榷,而且此时发生的面状分布的岩浆活动主要是由变质岩陆壳深融而成,缺少新生的幔源物质(王德滋和沈渭洲,2003;舒良树,2012;Wang Yuejun et al.,2013a),所以华南在早古生代经历的造山作用是没有洋壳参与的陆内造山作用。根据板溪群岩石属性及岩浆侵入的研究(李献华,2000;Xu Yigang et al.,2001;Shu Liangshu,2006),印支期华南大陆也发生了板块内相互作用的陆内造山作用。从中生代侏罗纪晚期开始,华南大陆开始受到西太平洋板块的俯冲,从而引发强烈的火山—岩浆侵入活动,目前华南大陆东部主要分布的花岗岩大部分形成于新生代,由于挤压而在扬子地块中形成走滑兼逆冲推覆构造(舒良树,2012;张国伟等,2013)。综上所述,目前的研究认为扬子地块和华夏地块的接触关系经历了碰撞拼合—拉张伸展—多期陆内碰撞造山的三大构造过程。

本次通过8d、12g二维地电模型对扬子地块、江南造山带以及华夏地块这三大地质构造单元接触关系的研究中发现,各地质块体具有明显“横向分块”特征,在江南造山带北段靠近长沙附近,扬子地块与江南造山带都表现为高阻特征,特别是扬子地块从上地壳至上地幔都为稳定的高阻异常(如图10 8d剖面中 R1所示),而且江南造山带的中下地壳也为高阻特征(如图10 8d剖面中R2所示),我们推测高阻体R1及R2应为扬子前寒武纪结晶基底,高阻的扬子结晶基底已经越过江南造山带西界到达江山—绍兴断裂附近,与华夏地块挤压碰撞形成了江山—绍兴深大断裂,地质调查认为江山—绍兴断裂为逆冲推覆构造(张国伟等,2013;舒良树,2012),在电性模型中也得到了同样的结论。但江南造山带西边界断裂九江—石台断裂为小型断裂特征,并没有明显切穿地壳,据此推断扬子地块与华夏地块碰撞融合主要发生在江山—绍兴断裂附近。电性模型显示8d剖面经过的江南造山带中没有发现明显的岩浆上涌的通道,上地壳未被后期的岩浆活动大范围的改造,这也与Zhao Guochun(2015)对花岗岩的时空分布研究结果对应较一致。

图10 华南陆块东北部扬子地块与华夏地块接触关系的电性特征解释图Fig.10 Interpretation of electrical characteristics of the contact relationship between the Yangtze Block and the Cathaysia Block in the northeastern South China BlockJSF:江绍断裂;GDBF:赣东北断裂;JJ—STF:九江—石台断裂;TLF:郯庐断裂;JiSF:吉首断裂;ZC—LWF:郴州—临武断裂;TX—YT—AYF:屯溪—鹰潭—安远断裂;SY—ZH—DPF:上虞—政和—大浦断裂JSF:Jiangshan—Shaoxing fault;GDBF:Northeast Jiangxi fault;JJ—STF:Jiujiang—Shitai fault;TLF:Tancheng—Lujiang fault;JiSF:Jishou fault;ZC—LWF:Chenzhou—Linwu fault;TX—YT—AYF:Tunxi—Yingtan—Anyuan fault;SY—ZH—DPF:Shangyu—Zhenghe—Dapu fault

沿结合带往东北,在研究区结合带的北段部位展布的12g剖面中,扬子地块与华夏地块接触部位发现由上地幔延伸而上形成较大面积的低阻异常(如图10 12g剖面中C1、C2、C3所示),该区域的扬子地块、江南造山带以及华夏地块地壳都被地幔上侵物质严重破坏,特别是扬子地块和江南造山带破坏更是严重,已经无法确认九江—石台、赣东北等断裂深部延伸情况,只能在浅部看到断裂分布。在该区域由于深部热物质的上涌底侵而改造了扬子地块和华夏地块固有的块体接触关系,因为当下的MT测深研究是对现有的构造情况的一种电性反映,但是可以说明结合带的东北段部位受到的深部热扰动比西南段剧烈,目前研究认为造成这种差异是因为在146~130 Ma,太平洋板块沿着平行于中国大陆边缘方向向西南俯冲,首先改变了华南陆块东北部地壳厚度及结构(Li Xianhua et al.,2013;Zhang Yuzhi et al.,2013)。

综上所述,在研究区西南部江南造山带内,扬子地块与华夏地块存在较明显的块体碰撞挤压作用,且这样的挤压作用主要集中在江山—绍兴断裂附近,也是形成该断裂的主要深部原因,而往东北,特别是靠近赣东北断裂带附近存在明显的深部热物质上涌特征,主要为低阻异常,已改造扬子地块与华夏地块碰撞接触关系,我们推测太平洋板块的俯冲对研究区华南陆块东北部的影响较西南部剧烈。

5.2 赣东北成矿带深部机制讨论

江南造山带经历了多期构造造山成矿作用,在晋宁期等早期形成的大多数矿床已经被破坏改造,目前主要的矿床形成于中生代,在印支期进一步矿化富集,而在燕山期由于太平洋板块的俯冲导致大量花岗岩集中爆发侵入,形成大规模矿床(毛景文等,2004;Zhou Xinmin et al.,2006;Li Xianghua et al.,2007;毛景文等,2011)。本文涉及的赣东北成矿带现存的矿床主要分布于扬子、华夏两个古地块拼合处,该矿带形成与燕山期酸性火山侵入岩有关的铜、铅锌、金、银等矿床,主要受赣东北等NNE向断裂控制(徐磊等,2012)。

为了进一步研究赣东北成矿带深部机制关系,将穿过赣东北成矿带的12g剖面进行了深度为200 km的深部成矿机制研究(如图11所示)。陈毓川等(2014)认为赣东北断裂带与地幔沟通,部分幔源物质参与了该成矿带的成矿成岩,且成矿作用受断裂控制。虽然该矿带形成于燕山期,但成矿系统是在深部过程驱动下形成的、具有自组织的能量及物质迁移汇聚系统,系统在形成和演化过程中,在岩石圈不同尺度上留下“痕迹”,这种“痕迹”可以通过地球物理方法进行探测(吕庆田等,2020)。电性模型中发现赣东北成矿带深部存在明显的岩浆热液底侵现象,指示深部成矿物质从地幔上涌(如图11红色箭头所示),目前成矿研究认为赣东北成矿带成矿花岗岩为I型花岗岩,即幔源物质参与成矿,在175 Ma左右,由于伊泽奈奇板块向西北俯冲,使得江南造山带成为活动大陆边缘导致俯冲板块重熔,形成高钾钙碱质岩浆(毛景文等,2011),我们结合电性特征推测幔源成矿物质在上侵过程中由于赣东北等断裂的存在使得成矿岩浆中加入了较少地壳物质后形成德兴、银山、水平等铜矿,该上侵通道已经破坏或撕裂了赣东北成矿带高阻的上地壳,部分存在的上地壳厚度已经不足10 km,该通道中因超壳断裂存在,使得地壳物质参与的比重应该很小,形成的与地幔相通的高钾钙碱质岩浆为赣东北成矿带提供幔源成矿元素,并且成矿岩浆继续沿着因赣东北断裂形成的冲断层、层间滑脱及次级断裂网络等上侵,最终富集成铜、锌、金等矿床,吕庆田等(2019)通过对长江中下游成矿带成矿系统的探测同样发现成矿带普遍存在壳内巨型逆冲断裂系统等岩浆/流体从深部迁移的主要“通道”;而陆内造山期形成的断裂网络和岩体接触带等成为成矿系统“末端”矿质沉淀的“场所”。可以说发育在赣东北成矿带中下地壳的低阻物质区是导致赣东北成矿的深部原因,这也与研究区天然地震研究结果一致(He Chuansong and Santosh,2016) 。

图11 赣东北成矿带深部机制解释模型及德兴铜矿成矿模式简图(右图据毛景文等,2009修改)Fig.11 Interpretation model of deep mechanism of the Northeast Jiangxi Metallogenic Belt and the Metallogenic model of Dexing copper deposit (the right figure modified from Mao Jingwen et al.,2009&)

5.3 地球动力学意义

由华夏地块和扬子地块在元古代碰撞拼合而成了华南陆块,在古生代受到华北板块的挤压,自中生代以来受到了天平洋板块向西的俯冲影响,但是新生代构造叠加使得中生代太平洋板块俯冲的地质迹象被破坏,所以对华南中生代构造演化争议最多(毛建仁等,2014)。所以当下华南地区的构造模式是在多期构造活动下的总体反映,本文对地球动力学意义的讨论将以华南陆块东北部现今电性模型特征为基础,对华南现今构造特征影响最大、但争议最多的中生代岩浆活动做出对电性约束,讨论其地球动力学内容。

目前研究认为中生代以来的太平洋板块向欧亚大陆的俯冲导致了华夏地块从东南沿海到雪峰山近东西向的岩浆活动( Zhou Xinmin et al.,2006;Li Xianghua et al.,2007;刘凯等,2016;夏炎等,2016),并提出了多种俯冲模式假设。在研究区二维电性模型中华夏地块浅地表至10 km深度主要为高阻响应,推测这是研究区密集分布火成岩及华夏地块残存上地壳的电阻率特征响应。在实地地质调查中可能无法完全正确圈定花岗岩范围,但是从地表至地下10 km的电性特征来看华南陆块东部大部分范围都被侵入岩占据,这也说明华南陆块东部上地壳破坏严重,可能发生过大面积的板片重融。

通过8d和12g剖面对比发现华夏地块东北部地壳厚度明显小于其西南部,说明华夏地块东北部较西南部受到了更大的深部热扰动。目前的地球动力学研究认为在146~130 Ma,太平洋板块沿着平行于中国大陆边缘方向向西南俯冲,首先改变了华南陆块东北部地壳厚度及结构,俯冲洋壳熔融在长江中下游形成早期的矿化和埃达克岩侵入(Li Xianhua et al.,2013;Zhang Guowei et al.,2013),我们的模型支持这个观点,电性模型指示研究区东北部地壳相较于西南部厚度明显更小、受到的热扰动更剧烈、时间更久。在130 Ma左右继续俯冲,并于华南陆块东南部下方形成板片窗,导致软流圈广泛上涌使得大陆岩石圈逐渐变薄,地壳的大规模熔融使得华南形成早白垩世长英质钾质火山岩(Zhao Liang et al.,2016)。印支期以来华南大陆从古亚洲洋域到太平洋域的转变,使其大地构造格局由东西向为主转变为北东向为主的构造体系(Wang Yuejun et al.,2013a;毛建仁等,2013),在120 Ma左右,华南陆块前后受到太平洋板块俯冲发生80°偏转转向西北的正向俯冲体系(Sun Weidong et al.,2007;毛建仁等,2014;Zhao Liang et al.,2016)的强烈作用,而后俯冲的大洋板片在重力作用下逐步后撤,导致华南陆块东部出现拉张裂谷环境,我们的电性模型被屯溪—鹰潭—安远和上虞—大浦—政和等断裂切割的高阻微地块分散拼接特征支持这样的观点,裂谷环境使得华南陆块中出现了零星A 型花岗岩类,而且这些A型花岗岩类分布特征显示出由内陆向沿海、由南西向北东逐渐年轻化趋势(夏炎等,2016)。所以华夏地块沿海火成岩带中向洋、向东北大致体现出年轻化趋势说明了太平洋板块先俯冲后再后撤,后撤过程中弧后伸展作用控制了华南陆块东部岩浆活动内陆向沿海迁移的过程( Zhou Xinmin et al.,2006;He Zhenyu and Xu Xisheng,2012;Liu Lei et al.,2012,2014;夏炎等,2016),而且不同纬度板片后撤的时间也不同,在侏罗至白垩纪大致表现为南早北晚的非同步后撤(Liu Lei et al.,2016),这也是导致华南陆块中形成在侏罗至白垩纪形成大规模岩浆热液型矿床的主要深部动力学原因。

6 结论

本文通过华南陆块东北部大地电磁测深研究,得出了华南陆块东北部岩石圈尺度的二维及三维电阻率模型,通过地电模型的地质—地球物理研究,就研究区华南陆块各块体之间的接触关系、江南造山带深部构造特征及区域动力学意义等地质问题进行了讨论,得出以下结论:

(1)研究区西南段靠近长沙附近,在江南造山带内存在扬子地块与华夏地块碰撞挤压构造,高阻的扬子结晶基底已经越过江南造山带到达江山—绍兴断裂附近,与华夏地块挤压碰撞形成了江山—绍兴逆冲推覆深大断裂构造,而研究区东北段赣东北断裂附近,存在明显的深部物质上涌和底侵作用,已改造扬子与华夏地块原有接触关系。

(2)赣东北成矿带深部存在明显的岩浆热液上涌现象,幔源成矿物质在上侵过程中由于赣东北等断裂的存在使得成矿岩浆上侵过程中混入少量分布于赣东北断裂两侧的地壳物质,形成的与地幔相通的高钾钙碱质岩浆为赣东北成矿带提供幔源成矿元素,最终富集成铜、锌、金等矿床,中下地壳的低阻异常区是导致成矿的深部物源。

(3)华夏地块东北部地壳厚度明显小于其西南部,说明华夏地块东北部较西南部受到了更大深部热扰动,认为在146~130 Ma,太平洋板块沿着平行于中国大陆边缘方向向西南俯冲,首先改变了华南陆块东北部地壳厚度及结构,在120 Ma左右,太平洋板块俯冲发生80°偏转转向西北的正向俯冲体系后,而后俯冲的大洋板片在重力作用下逐步后撤,我们的电性模型高阻微地块的分散拼接特征支持这样的观点。后撤形成的裂谷环境使得华南陆块中出现了由内陆向沿海、由南西向北东逐渐年轻化零星A 型花岗岩类,大洋板块南早北晚的非同步后撤使得华南陆块在侏罗至白垩纪形成大规模岩浆热液型矿床的主要深部动力学原因。

猜你喜欢
电性扬子华南
江苏扬子电缆集团有限公司
河北省山区典型地层的电性特征及物探找水的方向
华南风采
2019南京扬子当代艺术博览会
记华南女院前三任校长
华南掠影
苏萌娜 初心不渝 情牵华南
用于燃烧正电性金属的合金的方法
岔路失羊
红雪