Ordovician island arc volcanic rocks in Halahu area of Qilian Mountain and their tectonic significance
-
摘要:
通过1:5万区域地质调查，对祁连山哈拉湖地区火山岩进行了野外地质、岩石学、锆石U-Pb同位素年龄、地球化学等研究。LA-ICP-MS锆石U-Pb同位素定年结果显示，该火山岩年龄为466.3±2.4Ma（n=9，MSWD=1.4），形成于中奥陶世。岩石地球化学研究表明，哈拉湖地区火山岩为低钾拉斑玄武系列，绝大多数样品的亲石元素Rb、Th、U、Ce、Zr、Nd和稀土元素La、Sm富集明显；高场强元素（Sr、P、Ti）强烈亏损，Zr、Hf微弱富集，Ba明显亏损。总的特征显示，少数火山岩具有过渡型洋脊玄武岩的特征，大多数火山岩样品显示岛弧火山岩的特征。此外，哈拉湖地区的岛弧火山岩与晚奥陶世的岛弧花岗岩伴生在一起。这些特征表明，与俯冲有关的奥陶纪岛弧岩石可能与拉脊山地区古大洋的闭合有关。
Abstract:Based on 1:50000 regional geological survey, the authors studied volcanic rocks of Hulahu area in Qilian Mountain in the aspects of field geology, petrology, geochronology and geochemistry. Zircon LA-ICP-MS U-Pb dating results show that the Halahu volcanic rocks were formed in Middle Ordovician (466.3 ±2.4Ma, MSWD=1.4). Geochemical data show that the Halahu volcanic rocks belong to low-K tholeiite series, and are characterized by enrichment of LILEs (Rb, Th, U, Ce, Zr and Nd) and LREEs and depletion of Sr, P, Ti and Ba. On the whole, a few samples have features similar to those of T-MORB, but most samples show affinities with island arc related to subduction. In addition, the Halahu volcanic rocks were associated with Late Ordovician arc granitoids. These features indicate that the island arc rocks were associated with Ordovician paleo-oceanic closure in Lajishan area.
-
Keywords:
- volcanic rocks /
- island arc /
- Ordovician /
- Halahu area, Qilian Mountain
-
徐淮推覆-褶皱带位于华北板块东南缘、郯庐断裂西侧（图 1），为该区重要的陆内变形带。区域内徐淮断裂-褶皱系统由北向南走向由NNE—近NS—NW变化，也称为徐淮弧。弧顶位于萧县—淮北市一带, 与华北板块的其余构造形态具有明显差别（图 1）。早在20世纪90年代，徐树桐等[3]将其称为半圆形造山带或大型褶皱逆冲推覆体；王桂梁等[4]提出区域构造表现为弧形双冲-叠瓦状逆断层特点；彭凌日等[5]与Shu等[1]认为，区域构造属于挤压背景下的薄皮构造模型，表现为区域上普遍存在一系列弧形的斜歪褶皱、逆冲推覆构造。以上均为该区构造样式成因的研究，取得较多成果，而对于区域广泛存在的滑脱层研究不足[1, 1-5]。
滑脱层的存在影响盖层与基底间界面的剪切强度，使褶皱-冲断带构造传播方式及构造样式发生变化[6]。通常滑脱层能使构造变形传播距离更远，并且不同厚度滑脱层影响后冲断层的发育数量及构造的对称性，形成隔挡隔槽式褶皱、等间距逆冲断层、反冲断层等典型构造[7-9]；滑脱层深度及强弱性质的改变将影响构造变形的传播特征及褶皱、断层构造的样式和间距变化[10-13]。
随着构造物理模拟手段的不断完善，它已成为地质学家认识构造变形过程、分析构造成因机制的重要手段[14-19]。本文选取徐淮推覆-褶皱带为研究对象，结合前人研究资料，采取物理模拟手段对区域滑脱层深度变化、滑脱层摩擦系数的控制作用进行探索，进而对区域构造样式的控制因素和形成机制进行研究。
1. 徐淮地区构造背景
1.1 区域地层
研究区沉积盖层在东西向剖面上表现为东侧盖层厚度大于西侧（图 2-b），基底为太古宇泰山群变质岩，盖层分为3个亚层，即青白口系—震旦系、寒武系—奥陶系、石炭系—二叠系，相应亚层的分布表现为下部的青白口系—震旦系多为碎屑岩-碳酸盐岩沉积，岩石组成多为叠层石灰岩与白云岩，中部的寒武系—奥陶系主要为浅海碳酸盐岩，石炭系—二叠系以滨海相灰岩及滨海湖泊-三角洲体系陆源碎屑岩沉积为主，其中山西组和上、下石盒子组是主要的含煤层系①[21-22]。结合前人野外调查，区域上滑脱层普遍发育，主要为页岩、膏岩泥灰岩及煤层，分别为青白口系泥灰岩、寒武系猴家山组与馒头组页岩、奥陶系贾汪组页岩与马家沟组膏岩、石炭系—二叠系页岩及煤层（图 3）[1, 23]。其中，青白口系泥灰岩、寒武系页岩、石炭系—二叠系页岩及煤层作为主要拆离层控制区域主要构造变形[24]。
![]()
1.2 构造特征
研究区以逆冲推覆构造为主，未卷入结晶基底，属于典型的薄皮构造（图 2-a），平面上逆冲推覆构造为一系列走向由北向南、由NNE—近NS—NW变化的逆断层及伴生斜歪褶皱[5]，如白土镇东斜歪褶皱的轴面产状为144°∠75°，高角度逆断层的产状为137°∠71°，低角度逆断层的产状为124°∠41°（图 4），表明该区受到从南东向北西的挤压应力场作用。
结合前人在该区域进行的构造应力场研究[1, 25-26]，区域受2期应力场作用，分别是早燕山期SE—NW向挤压及晚燕山期NNE—SSW向挤压力。本次研究通过61个野外共轭节理、断面擦痕求应力场，并利用赤平投影求得最大主应力轴（σ1）优势方位，2个极密点方位为200°∠6°、133°∠6°，表明该2期应力场方位与前人研究一致（图 5）。虽然晚燕山期受NNE向挤压，导致该区侏罗纪—白垩纪盆地发育，但仍以早燕山期SE—NW向挤压作用形成的构造格局为主[4, 27]。
![]()
图 5 区域最大应力场统计极密图（A与B代表研究区最大主应力轴机密点）Figure 5. Statistics of the region's maximum stress field以皇藏峪背斜为界将研究区构造划分为东部构造带与西部构造带（图 2-b），东部构造带主要为一系列斜歪褶皱，少量叠瓦状分支断层；西部构造带表现为强烈变形，发育一系列逆冲断层。总体构造具明显东西分带的特点[3, 20]。
2. 构造物理模拟实验设计
2.1 实验材料
构造物理模拟实验遵循相似原则，即几何学相似、动力学相似、运动学相似，实验材料物性、模型尺寸等需要与自然界中实际的相应参数保持一定比值[28-29]。本次实验中所用材料为石英砂、微玻璃珠及硅胶。
干燥石英砂用于模拟脆性的上覆岩层，其物理性质为脆性且遵循库伦-莫尔破坏准则，内摩擦角为31°左右，内聚力极小[30]，变形过程中铺设的石英砂颗粒将出现从应变硬化-应变软化的过渡，导致影响断层角度的内摩擦系数稍小于未变形前岩石的内摩擦系数[31-32]。石英砂颗粒直径为200~300μm，根据不同的铺沙方式（筛或直接倒入模型中），密度会有差别，总体密度在1.3~1.6g/cm3之间。人工烧结的彩色石英砂颗粒一般物性不发生改变，铺设1mm厚彩色石英砂于模型中作为观察构造变形的标志层[33]。
微玻璃珠通常用于模拟强滑脱层，表面光滑且磨圆度好，玻璃珠内聚力几乎为零，内摩擦角为25°左右[34-35]，应变软化幅度很小，变形前后内摩擦角变化不大。本次实验使用直径为300μm的微玻璃珠充当强滑脱层，铺设于模型最底部[36]。
硅胶为透明高粘度的材料，用于模拟上地壳的塑性变形，在低应变速率下具有牛顿流体性质且具有非常低的屈服强度[37]。测得室温下本实验材料粘度为1.2×104Pa.s, 密度为0.926g/cm3, 铺设2~3mm于基底之上，作为弱滑脱层，控制变形传播距离，产生各类褶皱-冲断构造[7, 38-39]。
2.2 相似条件
根据相似原理，模型与地质原型长度的相似因子L*为1×10-5（即实验室中模型1cm代表自然界中1km的长度）；地质原型中沉积岩的平均密度为2.6g/cm3, 模型中石英砂与微玻璃珠的平均密度为1.4g/cm3, 所以实验设定的密度相似因子为0.5[40]；模型与地质原型的重力加速度相同，为9.8m/s2, 重力加速度相似因子为1；关于硅胶滑脱层，粘度一般为1×104Pa.s, 而地质原型相应地层的粘度为1×1019Pa.s, 粘度相似因子为1×10-15（*代表相应模型参数与地质原型参数的比值）。
根据前人对徐淮地区变形的年代学数据及缩短率的研究[23, 25]，推算出地质原型变形速率为1~10mm/a，因而设置实验室计步器推挤速率为0.005mm/s，结合速度相似因子v*=5×104，实验中模拟地质原型的挤压速率为3.2mm/a，与推算变形速率（1~10mm/a）大小为同一数量级，符合构造模拟相似条件。
2.3 实验装置
根据控制变量法来改变滑脱层深度及滑脱层物理性质[6, 17]，设计2组实验。第一组实验设计2个模型分别为模型1与模型2，初始尺寸均为600mm×400mm×24mm，且底板水平。模型1使用4mm厚微玻璃珠作为强滑脱层，20mm厚的石英砂作为盖层，而模型2则用一层厚4mm的硅胶层为弱滑脱层，盖层材料与厚度均与模型1相同（图 6-a、b）。模型两侧为钢化玻璃，为减少边界效应的影响，在实验前用无水酒精将其擦拭干净。
第二组实验设计2个模型分别为模型3和模型4，尺寸均为600mm×400mm×24mm, 模型底板均保持水平，考虑到滑脱层深度变化的影响，模型中均设计有先存的台阶式隆起，模型3距离活动端430mm处设置有5mm高的台阶状基底隆起，而模型4则在距活动端430mm处设置为10mm高的台阶状隆起。模型3、4底部均有一层厚4mm的硅胶层作为滑脱层，盖层为石英砂，两模型材料厚度均为24mm（图 6-c、d）。
4个模型均从右侧施加挤压力，活动端推板运动速度均为0.005mm/s，根据前人对区域平衡剖面恢复所得缩短率[1, 25]，本文模型缩短率取30%，每组实验过程通过照相机定时，每隔5min照相记录。
3. 实验结果及分析
3.1 第一组实验
3.1.1 模型1
实验过程分为4个阶段。阶段一，缩短率为2.00%时，在靠近活动端出现前冲断层1及断层相关褶皱，形成构造单元Ⅰ；阶段二，缩短率为13.50%时，构造单元Ⅰ中发育1号断层的分支断层2，以及从微玻璃珠滑脱层逐渐向地表延伸的断层3、4，此时该单元具逆冲叠瓦构造特点，位于构造单元Ⅰ前方1.6cm处发育受前冲断层5控制的构造单元Ⅱ；阶段三，缩短率为20.50%时，构造单元Ⅱ形成与构造单元Ⅰ相同的构造样式，其中断层6为断层5的分支断层，调节断层运动，位于构造单元Ⅱ前方3.8cm处(相比前者间距更远)发育受前冲断层7控制的构造单元Ⅲ；阶段四，缩短率为30.00%时，构造单元Ⅲ形成构造与前者相同，其中断层9为分支断层（图 7）。
![]()
图 7 物理模拟模型1演化图a~d—缩短率分别为2.00%、13.50%、20.50%、30.00%；e—缩短率为30.00%时沿挤压方向的横切剖面图；
1~9代表相应断层发育顺序；Ⅰ~Ⅲ代表相应构造单元发育顺序Figure 7. The sequential photographs showing deformation evolution of model 1实验所切剖面见3个间距明显不同的构造单元。构造单元Ⅰ主要受1~4前冲断层控制，形成堆垛式逆冲叠瓦构造及断层相关褶皱；构造单元Ⅱ与构造单元Ⅲ构造样式与构造单元Ⅰ相同，各个构造单元内断层间距变化均表现为朝固定端（前陆）逐渐增大，与前人研究结果一致[14, 41]（图 7-e、图 8-f）。
![]()
图 8 徐淮推覆-褶皱带物理模拟中模型1(a)、模型2(b)的断层角度/缩短率关系图(1~9代表相应断层编号)、模型1与2楔长/缩短率（c）、楔高/缩短率(d)、楔角/缩短率(e)对比图及最终剖面上断层间距/断层顺序（f）对比图（断层顺序1~8代表活动端向固定端测量断层间距的顺序）Figure 8. Model 1(a), model 2(b) fault angle/shortening ratio diagram, model 1 and 2 wedge long/short ratio(c), wedge high/ short ratio(d), wedge angle/short ratio(e) contrast diagram and the final section on fault distance/order(f) contrast diagram随着缩短率逐渐增加，楔长变化表现为增长至逐渐稳定的特点，楔长增加代表构造变形向固定端（前陆）传播，直至近稳定状态，最前端的断层前方未出现新构造使楔长不变或小幅减小；楔高则随缩短率增加逐渐增加，但增高速率逐渐减小（图 8-c、d），与前人观点一致[42]。组成构造单元Ⅰ的断层1、2、3、4角度变化具有较好的相关性，表现为晚出现的断层角度变化均会对早出现的断层角度产生影响；组成构造单元Ⅱ的断层5、6中，断层6角度逐渐小幅度增加，而断层5角度先增大后稳定，随后按该变化趋势发生变化；组成构造单元Ⅲ的断层7、8、9断层角度变化速率最慢，断层角度缓慢增加直至稳定，可见断层越远离力源，角度变化越容易达到稳定状态，并受晚期新构造影响越小(图 8-a)。
地表楔形隆起主要受到基底滑脱面性质及上覆材料强度的影响[14, 43]。随缩短率增加，活动端产生的楔状隆起角度达到临界角度之后，由于前缘未出现新断层，使角度继续增大向超临界角发展，直至前缘出现新断层，楔形隆起角度将减小，并向着次临界角值发展，达到相应次临界角后将再次增大（图 8-e），可见楔体达到临界角度后仍处于非稳定状态，而不是保持该状态向固定端方向稳定传播。模型1中，缩短率为8.75%时，楔角增大至相应超临界角，随后向一定次临界角值减小，缩短率为13.75%、23.75%时，相应状态下相应次临界角值再增大，符合前人观点[44-46]。
3.1.2 模型2
按构造单元将实验过程分为5个阶段。阶段一，缩短率为2.00%时，形成由前冲断层1与断层相关褶皱组成的构造单元Ⅰ；阶段二，缩短率为13.00%时，构造单元Ⅰ中发育后冲断层，此时断层相关褶皱形态为不对称箱状，位于构造单元Ⅰ前方，构造单元Ⅱ出现，主要由一对前后冲断层及断层相关褶皱组成；阶段三，缩短率为18.25%时，构造单元Ⅰ、Ⅱ继续发育，构造单元Ⅲ前方出现由只受前冲断层控制的构造单元Ⅲ；阶段四，缩短率为23.25%时，构造单元Ⅲ未出现后冲断层，断层相关褶皱形态为尖棱斜歪，并在其前方出现由前冲断层组成的构造单元Ⅳ；阶段五，缩短率为30.00%时，构造单元Ⅳ继续发育，表现为宽阔斜歪褶皱(图 9)。
![]()
图 9 物理模拟模型2演化图a~e—缩短率分别为2.00%、13.00%、18.25%、23.25%、30.00%；
f—缩短率为30.00%时沿挤压方向的切片图；1~5代表相应断层发育顺序；
Ⅰ—Ⅴ代表相应构造单元发育顺序(1’、2’、4’、5’代表与1、2、4、5同时生成的断层)Figure 9. The sequential photographs showing deformation evolution of model 2实验所切剖面见5个明显构造单元。构造单元Ⅰ、Ⅱ、Ⅳ、Ⅴ均由前、后冲断层及断层相关褶皱组成，其中褶皱为不对称箱状；构造单元Ⅲ由前冲断层及断层相关褶皱组成，其中褶皱表现为尖棱斜歪特点。断层间距均大于模型1，表现为先减小后增大的特点（图 8-f、图 9-f）。
实验过程中楔长变化主要为增长-稳定交替变化特点，其中稳定阶段时间略小于模型1，但模型2楔长普遍大于模型1，表明模型2构造变形传播的更远；楔高逐渐增大，至缩短率13.75%达到近稳定状态，随后楔高小幅度变化；活动端地表隆起楔角在缩短率为6.25%时，达到相应超临界角后，楔角下降，随缩短率增加，楔角减小速率逐渐变小，直至楔角趋于近稳定状态，达到临界值。实验过程中模型2楔角值均小于模型1，符合Davis等关于基底摩擦系数对库伦楔角影响的认识(图 8-c、d、e)[45]。各断层角度变化趋势相似，均表现为增大-稳定（图 8-b）的变化趋势。
3.2 第二组实验
3.2.1 模型3
实验过程分为6个阶段。阶段一，缩短率为2.00%时，最先在活动端处出现构造单元Ⅰ，表现为一对前冲后冲断层，地表微微隆起；阶段二，缩短率为11.50%时，构造单元Ⅰ断层持续活动，断层相关褶皱表现为不对称箱状特点，在构造单元Ⅰ前方出现受前冲断层控制的构造单元Ⅱ；阶段三，缩短率为16.50%时，构造单元Ⅱ出现后冲断层，使断层相关褶皱表现为箱状特点，并在前方构造单元Ⅲ开始发育受前冲断层控制的断层相关褶皱，同时滑脱层深浅变化处发育构造单元Ⅲ’的断层转折褶皱；阶段四，缩短率为25.50%时，构造单元Ⅰ、Ⅱ持续活动，与Ⅰ、Ⅱ构造演化特点相同，在构造单元Ⅲ与Ⅲ’之间发育构造单元Ⅳ的受前冲断层控制的断层相关褶皱，褶皱形态宽缓；阶段五，缩短率为29.50%时，在构造单元Ⅲ’前方发育构造单元Ⅴ的前冲断层及断层相关褶皱，褶皱波长较小；阶段六，缩短率为31.70%时，构造单元Ⅳ持续活动，但未见后冲断层，构造单元Ⅴ的褶皱高度明显增加，相应断层延伸至地表(图 10)。
![]()
图 10 物理模拟模型3演化图a~f—缩短率分别为2.00%、11.50%、16.50%、25.50%、31.70%；
g—缩短率为31.70%时沿挤压方向的切片图；Ⅰ—Ⅵ—代表相应构造单元发育顺序(其中Ⅲ’代表与Ⅲ同时生成的构造单元)Figure 10. The sequential photographs showing deformation evolution of model 3实验所切剖面见7个明显的构造单元，构造单元Ⅰ、Ⅱ、Ⅲ主要受前冲断层及伴生的后冲断层控制，形成箱状褶皱；构造单元Ⅳ只受前冲断层控制，断层相关褶皱为宽缓斜歪褶皱，具明显不对称性；构造单元Ⅲ’、Ⅴ、Ⅵ受前冲断层控制，其中构造单元Ⅱ、Ⅰ’表现为断层转折褶皱，构造单元Ⅴ、Ⅵ组成逆冲叠瓦构造。深滑脱层区域各构造单元间距较大，平均间距为6.5cm，发育4个构造单元；浅滑脱层区域构造单元间距较小，从构造单元Ⅲ’至Ⅵ间距为先增大后减小；滑脱层深度变化处则对应剖面中波长最大的向斜（图 10-g、图 11-c）。
![]()
图 11 徐淮推覆-褶皱构造带物理模拟中模型3(a)、模型4(b)相应构造单元隆起高度/缩短率关系图和相邻构造单元间距/次序关系图（c）（次序1~6代表从活动端向固定端测量断层间距的顺序；Ⅰ~Ⅴ代表构造单元先后顺序）Figure 11. Model 3 (a) and model 4 (b) uplift height of the corresponding tectonic units/shortening rate diagram and adjacent structural unit spacing/sequence diagram (c)实验中随着缩短率逐渐增大，构造单元Ⅰ、Ⅱ、Ⅲ隆起高度变化为增大-稳定不断循环的特点，证明模型中导致构造活动的应力集中过程是不稳定的，使得相应构造单元的活动具有乱序性，但处于滑脱层深浅变化附近的构造单元Ⅲ’与Ⅳ的隆起高度持续增加，未表现明显的稳定期，其中晚期出现的Ⅲ’构造单元最终隆起高度甚至超过了构造单元Ⅱ、Ⅲ（图 11-a）。
3.2.2 模型4
实验过程主要划分为5个阶段。阶段一，缩短率为1.75%时，构造单元Ⅰ发育由后冲断层控制的断层相关褶皱；阶段二，缩短率为8.75%时，构造单元Ⅰ持续活动，形成明显斜歪褶皱，位于滑脱层深浅变化处见构造单元Ⅱ的断层转折褶皱；阶段三，缩短率为15.00%时，在构造单元Ⅰ与Ⅱ前方同时发育构造单元Ⅲ、Ⅲ’，两者均表现为受前冲断层控制的断层相关褶皱，构造单元Ⅰ出现2条不同深度的前冲断层，表现出叠加褶皱的特点；阶段四，缩短率为21.75%时，构造单元Ⅱ与Ⅲ’之间出现构造单元Ⅳ，表现为断层相关褶皱；阶段五，缩短率为30.00%时，在构造单元Ⅲ前方出现构造单元Ⅴ，该构造受前后冲断层控制属于冲起构造，构造单元Ⅳ进一步发育，表现为尖棱斜歪的断层相关褶皱(图 12)。
![]()
图 12 物理模拟模型4演化图a~e—缩短率分别为1.75%、8.75%、15.00%、21.75%、30.00%；f—缩短率为30.00%时沿挤压方向的切片图，
Ⅰ~Ⅴ代表相应构造单元发育顺序(其中Ⅲ’代表与Ⅲ同时生成的构造单元)Figure 12. The sequential photographs showing deformation evolution of model 4实验所切剖面见6个明显的构造单元，构造单元Ⅰ中晚期低角度前冲断层将早期一对前后冲逆断层切断，最终表现褶皱为共轴叠加褶皱；构造单元Ⅲ’与Ⅳ主要受前冲断层控制，其中构造单元Ⅲ’为尖棱斜歪褶皱，Ⅳ为箱状斜歪褶皱；构造单元Ⅱ、Ⅲ受前冲断层控制，其构造单元Ⅱ主要表现为断层转折褶皱，构造单元Ⅲ为逆冲叠瓦构造及断层相关褶皱；构造单元Ⅴ受前后冲断层控制，属于冲起构造。深滑脱层区域各构造单元间距较大，向固定端间距逐渐增加，发育3个构造单元；浅滑脱层区域构造单元间距较小，且向固定端逐渐减小；滑脱层深度变化处也对应剖面中明显波长最大的向斜(图 11-c、图 12-f)。
随着缩短率的增加，构造单元Ⅰ、Ⅱ、Ⅲ及Ⅲ’隆起高度先增大后趋于稳定，相应构造单元高度达到稳定状态后将不会改变，构造单元Ⅳ隆起高度增长最快，最终高度超过早于其出现的构造单元Ⅲ与Ⅲ’（图 11-b）。
4. 讨论
4.1 滑脱层强弱及深度对构造变形的影响
模型1基底之上为强滑脱层（微玻璃珠），随着缩短率增加，构造样式为逆冲叠瓦构造，表现为前展式传播形式，最终传播距离16.5cm，而模型2基底之上为弱滑脱层（硅胶），构造主要表现为不对称箱状褶皱及尖棱褶皱，构造传播方式与模型1相同，但变形向前缘方向传递距离明显大于模型1（图 7、图 9）。可见，滑脱层强弱差异控制着不同的构造样式。
模型3、4所设置装置底部不同高度台阶式抬升，使滑脱层深度发生变化，从而改变上覆盖层厚度与脆性层强度，使构造变形传播不再是传统的前展式，而呈乱序式，且在滑脱层深浅变化处容易应力集中而较早发生变形[17]，变形传播距离相较于模型2更远。对于模型3与4，由于基底隆起高度大小差异，模型3（低基底隆起条件）在深滑脱层区域将出现更多的构造单元；模型3在深滑脱层区域普遍出现受前后冲断层控制的箱状断层相关褶皱，模型4除构造单元Ⅰ发育共轴叠加褶皱外，均表现为受前冲断层控制的斜歪褶皱特点，在浅滑脱层区域两者构造样式基本相同，主要为断层相关褶皱与逆冲断层(图 10、图 12)。因此，滑脱层深度变化对盖层相关构造变形具有明显影响。
4.2 实验结果与实际构造变形对比
徐淮地区构造变形属于典型挤压背景下的薄皮构造，发育于早燕山期的滑脱构造，具有明显东西分带的特点，由东向西变形强度逐渐增强。以皇藏屿背斜为界，该界以西构造变形涉及盖层主要为寒武系—奥陶系与石炭系—二叠系，盖层厚度较小，发育一系列逆冲断层；该界以东构造变形涉及盖层主要为震旦系—奥陶系，盖层较厚，广泛发育轴面向SE方向倾斜的斜歪褶皱，叠瓦状分支断层较少[1]。
通过4个模型所切剖面与实际构造变形对比，认为实验模型4与区域上剖面具有较好的相似性，模型基底抬升处产生的背斜与皇藏峪背斜位置相当，虽然区域上变形含有多层滑脱层，但结合构造变形变化特点与物理模拟实验结果，认为区域上构造变形主要受2层不同深度滑脱层控制，分别是青白口系泥灰岩与寒武系页岩。以滑脱层抬升处（皇藏屿背斜）为界线，出现的褶皱间距及构造强弱的变化在模型4中有良好表现。滑脱层深度小，出现较密集的逆断层；滑脱层深度大，则构造形态以大波长的斜歪褶皱为主，断层未延伸至地表（图 13）。
![]()
图 13 物理模拟模型4（b）与徐淮地区构造剖面（a）对比图Figure 13. comparison diagram of model 4 (b) and the construction section (a) of the Xu-Hai area本文对徐淮地区构造地质原型设计了2组实验，通过改变模型的材料及边界条件，研究了滑脱层性质及滑脱层深度变化对构造变形的控制作用，取得了与地质原型相似的模拟结果。然而自然界构造变形复杂，如多期构造叠加、沉积与剥蚀作用、多层滑脱层等。本次研究从主要控制因素入手，对地质原型进行模型化与理想化。该区存在多层滑脱层，滑脱层平面展布未能考虑，因而对于盖层变形控制不全面，有待进一步探讨滑脱层平面展布等因素对区域构造的影响。
5. 结论
(1) 滑脱层控制徐淮地区薄皮构造的形成，滑脱层的存在使得变形向前传播距离更远，是形成徐淮地区构造样式的主控因素，缺少滑脱层构造变形则表现出一系列堆垛式逆冲断层特点。
(2) 徐淮地区存在多套滑脱层，模拟实验表明，控制区域上东西构造样式变化的滑脱层主要有2套，分别为青白口系泥灰岩（深滑脱层）和寒武系页岩（浅滑脱层）。以皇藏屿背斜为界线，界线以西，滑脱层深度较浅，构造变形强烈，构造间距密集，逆冲断层普遍发育；界线以东，滑脱层深度较深，构造间距较大，主要发育斜歪褶皱，断层未能延伸至地表。
致谢: 野外工作中得到哈拉湖西项目组全体成员的大力帮助，成文过程中北京大学宋述光老师给予了多次指导，审稿专家提出了宝贵的修改意见，北京大学董金龙博士给予帮助，在此表示诚挚的谢意。 -
![]()
![]()
图 3 哈拉湖地区火山岩野外特征（a）和镜下照片(b)
Pl—斜长石；Chl—绿泥石；Prx—辉石
Figure 3. Field characteristics (a)and microscopic photographs(b)of volcanic rocks in Halahu area
![]()
图 4 哈拉湖火山岩SiO2-K2O图解（a）和AFM图解（b）
Figure 4. SiO2-K2O (a) and AFM (b) diagrams for volcanic rocks in Halahu area
![]()
图 5 哈拉湖火山岩稀土元素配分图(a)和原始地幔标准化微量元素蛛网图(b)（底图据参考文献[15]）
Figure 5. Chondrite-normalized REE patterns (a) and primitive mantle-normalized trace elements spider diagrams (b) for volcanic rocks in Halahu area
![]()
图 6 哈拉湖地区火山岩（IU-PbC1014）锆石U-Pb谐和图（a）和206Pb/238U表面年龄加权平均值图（b）
Figure 6. Zircon U-Pb concordia diagram (a) and weighted average 206Pb/238U age (b) for volcanic rocks(IU-PbC1014)in Halahu area
![]()
图 7 哈拉湖火山岩构造环境判别图解
a—Ta/Yb-Th/Yb图解；b—Hf/3-Th-Ta图解。SHO—钾玄岩；CAB—钙碱性岩系；TH—拉斑岩系；TR—过渡的；VAB—火山弧玄武岩；IAT—岛弧拉斑玄武岩；ALK—碱性玄武岩；WPB—板内玄武岩；MORB—洋中脊玄武岩；D—火山弧玄武岩；C—碱性板内玄武岩；B—E型MORB/板内拉斑玄武岩；A—N型MORB
Figure 7. Tectonic discrimination diagrams for volcanic rocks in Halahu area
表 1 哈拉湖地区火山岩地球化学分析结果
Table 1 Geochemical compositions of volcanic rocks in Halahu area
样品号 ⅠC1014-2 ⅠC1014-3 IC1014 ⅠPm103/2-1 ⅠPm103/16-1 ⅠPm103/32-1 ⅠPm103/32-2 ⅠPm103/32-5 岩性 玄武岩 玄武岩 玄武岩 玄武安山岩 玄武安山岩 安山岩 安山岩 安山岩 SiO2 46.24 45.71 46.61 50.44 50.45 47.70 48.00 48.89 TiO2 0.38 0.40 0.58 0.83 0.90 0.66 0.78 0.53 Al2O3 7.57 7.97 12.39 12.19 12.27 10.51 11.70 10.51 Fe2O3 1.54 1.70 1.15 2.23 2.14 0.80 1.08 1.16 FeO 6.87 6.83 8.23 7.65 7.67 7.43 7.13 7.37 MnO 0.17 0.18 0.21 0.17 0.18 0.30 0.18 0.15 MgO 20.45 20.59 10.44 9.07 8.99 9.07 10.14 13.66 CaO 9.64 9.22 7.34 10.08 10.11 9.99 9.89 6.28 Na2O 0.49 0.49 0.17 2.10 2.34 0.89 2.07 0.19 K2O 0.19 0.18 0.42 1.59 1.43 0.04 0.33 0.03 P2O5 0.14 0.14 0.17 0.17 0.23 0.20 0.16 0.18 烧失量 5.11 5.37 11.11 2.40 2.27 11.42 7.65 10.25 总量 99.66 99.67 99.68 99.76 99.74 99.77 99.74 99.75 Mg# 82.95 83.00 68.16 65.41 65.25 67.59 70.49 75.65 A/NK 7.48 7.96 16.88 2.36 2.27 6.97 3.11 30.46 A/CNK 0.73 0.81 1.56 0.89 0.88 0.96 0.95 1.62 La 7.46 5.60 10.57 16.1 18.1 14.80 11.03 10.44 Ce 14.69 13.24 21.44 32.1 36.5 29.16 22.79 21.15 Pr 2.10 2.06 2.75 4.13 5.14 3.77 3.07 2.76 Nd 9.44 9.15 10.91 16.8 21.6 15.74 13.25 11.74 Sm 2.11 2.19 2.49 3.62 4.71 3.39 3.07 2.61 Eu 0.59 0.61 0.49 0.963 1.21 0.89 0.96 0.78 Gd 2.09 2.19 2.15 3.62 4.42 3.16 3.35 2.69 Tb 0.33 0.34 0.40 0.586 0.707 0.52 0.57 0.43 Dy 1.96 2.05 2.44 3.27 3.73 3.12 3.49 2.68 Ho 0.39 0.40 0.52 0.680 0.762 0.59 0.73 0.53 Er 1.14 1.13 1.46 1.85 2.09 1.70 2.10 1.50 Tm 0.19 0.18 0.24 0.297 0.315 0.26 0.32 0.24 Yb 1.18 1.13 1.51 2.00 2.04 1.47 2.02 1.43 Lu 0.18 0.17 0.23 0.319 0.300 0.21 0.27 0.20 Sr 59.17 57.33 145.4 428 449 241.2 272.6 156.3 Rb 3.73 1.48 14.50 35.9 28.9 3.18 9.63 1.58 Ba 37.09 26.84 111.4 380 440 23.70 104.2 36.40 Th 3.85 3.18 4.08 4.96 5.68 4.10 3.79 3.76 Ta 0.35 0.23 0.23 1.28 0.82 0.41 0.44 0.33 Nb 3.33 2.10 3.58 5.53 6.12 4.58 5.18 3.59 Zr 33.03 32.48 67.9 67.7 75.0 68.91 71.67 51.98 Hf 0.98 1.13 2.9 2.04 2.36 3.03 2.61 2.15 Cr 1177 1281 331.8 275 279 711.8 831.6 708.2 V 133.6 135.7 170.6 273 252 170.3 197.3 187.0 U 0.83 0.78 1.07 1.03 1.24 1.09 1.11 0.81 Y 10.27 10.37 13.33 17.40 19.09 14.92 17.93 13.55 Sc 24.5 24.89 27.91 38.2 41.2 26.25 34.28 39.33 注：主量元素含量单位为%，微量、稀土元素含量单位为10-6 
下载: 导出CSV
表 2 哈拉湖地区火山岩（IU-PbC1014）锆石U-Th-Pb同位素测定结果
Table 2 Zircon U-Th-Pb isotope analyses for volcanic rocks (IU-PbC1014) in Halahu area
测试号 Pb U 206Pb/238U 1σ 207Pb/235U 1σ 207Pb/206U 1σ 208Pb/232U 1σ 232Pb/238U 1σ 年龄/Ma 10-6 206Pb/238U 1σ 207Pb/235U 1σ 1 12 162 0.07545 0.00041 0.58197 0.00920 0.05594 0.00087 0.03014 0.00023 0.25267 0.00072 469 3 466 7 2 50 242 0.18414 0.00094 1.92580 0.01367 0.07585 0.00052 0.06177 0.00014 0.61816 0.00325 1090 6 1090 8 3 22 257 0.07385 0.00049 0.57869 0.01587 0.05683 0.00130 0.03671 0.00126 0.59360 0.00911 459 3 464 13 4 27 304 0.02979 0.00021 0.41911 0.00597 0.10205 0.00134 0.01808 0.00014 0.37823 0.00066 498 3 652 9 5 12 150 0.08039 0.00042 0.90074 0.01293 0.08127 0.00111 0.03683 0.00040 0.37754 0.00088 480 3 541 9 6 30 348 0.07976 0.00042 0.90553 0.01919 0.08234 0.00165 0.03635 0.00067 0.36581 0.00102 495 3 655 14 7 29 338 0.07888 0.00041 0.80625 0.01234 0.07413 0.00108 0.03137 0.00036 0.40672 0.00074 489 3 600 9 8 50 651 0.07551 0.00042 0.58558 0.00430 0.05625 0.00044 0.03269 0.00025 0.27742 0.00099 469 3 468 3 9 28 333 0.07534 0.00040 0.58593 0.01031 0.05640 0.00095 0.02570 0.00011 0.65253 0.00045 468 2 468 8 10 51 583 0.07905 0.00041 0.84468 0.00868 0.07750 0.00080 0.03047 0.00018 0.51798 0.00059 490 3 622 6 11 39 99 0.37093 0.00208 6.41381 0.03874 0.12541 0.00073 0.11708 0.00039 0.28248 0.00033 2034 11 2034 12 12 15 158 0.07948 0.00048 0.94601 0.01960 0.08633 0.00188 0.03708 0.00025 0.57647 0.00213 493 3 676 14 13 10 70 0.14135 0.00078 1.62255 0.05211 0.08325 0.00264 0.07423 0.00183 0.20732 0.00080 852 5 979 31 14 75 733 0.07966 0.00041 1.49718 0.00952 0.13632 0.00082 0.06244 0.00012 0.43148 0.00037 494 3 929 6 15 18 218 0.07444 0.00039 0.58308 0.00954 0.05681 0.00087 0.03029 0.00042 0.57697 0.00333 463 2 466 8 16 54 620 0.07514 0.00041 0.58008 0.01090 0.05599 0.00101 0.02416 0.00028 0.85623 0.00264 467 3 465 9 17 90 418 0.19838 0.00144 2.83062 0.02120 0.10349 0.00083 0.08891 0.00107 0.34191 0.00144 1167 8 1364 10 18 20 252 0.07482 0.00040 0.58176 0.00953 0.05639 0.00092 0.02340 0.00011 0.50771 0.00041 465 2 466 8 19 95 477 0.18679 0.00099 2.65184 0.02096 0.10296 0.00089 0.07092 0.00022 0.32939 0.00432 1104 6 1315 10 20 34 148 0.20761 0.00107 2.34458 0.01698 0.08191 0.00057 0.05787 0.00022 0.62972 0.00335 1216 6 1226 9 21 66 840 0.07541 0.00040 0.58401 0.00386 0.05617 0.00039 0.03284 0.00013 0.34566 0.00080 469 2 467 3 22 9 261 0.02893 0.00016 0.50801 0.00872 0.12737 0.00214 0.02838 0.00039 0.21948 0.00081 184 1 417 7 23 43 556 0.07486 0.00041 0.58123 0.00449 0.05631 0.00041 0.02048 0.00004 0.48987 0.00153 465 3 465 4 24 70 455 0.06704 0.00044 5.14905 0.04494 0.55703 0.00333 0.27698 0.00185 0.25872 0.00082 418 3 1844 16
下载: 导出CSV
-
黄汲清, 姜春发.从多旋回构造运动观点初步探讨地壳发展规律[J].地质学报, 1962, 42(2):105-151. http://www.cnki.com.cn/Article/CJFD1979-DZXE196202000.htm 李春昱, 刘仰文, 朱宝清, 等.秦岭及祁连山构造发展史[C]//国际交流地质学术论文集.北京: 地质出版社, 1978: 174-187. 肖序常, 陈国铭, 朱志直.祁连山古蛇绿岩的地质构造意义[J].地质学报, 1978, 54(4):287-295. http://www.cnki.com.cn/Article/CJFDTotal-DZXE197804002.htm 肖序常, 陈国铭, 朱志直.祁连山古板块构造的一些认识[J].地质科技, 1974, (3):73-78. 夏林圻, 夏祖春, 任有祥, 等.北祁连山石灰沟奥陶纪岛弧火山岩系岩浆性质的确定[J].岩石矿物学杂志, 1991, (1):1-10. http://www.cqvip.com/Main/Detail.aspx?id=642461 夏林圻, 夏祖春, 徐学义.北祁连山海相火山岩岩石成因[M].北京:地质出版社, 1996. Xia L Q, Xia Z C, Xu X Y. Magmageneisis in the Ordovician in back basins of the northern Qilian Mountains, China[J]. Geological Society of America Bulletin, 2003, 115:1510-1522. doi: 10.1130/B25269.1
冯益民, 何世平.祁连山大地构造与造山作用[M].北京:地质出版社, 1996. 宋述光.北祁连俯冲杂岩带的构造演化[J].地球科学进展, 1997, 12(4):351-365. doi: 10.3321/j.issn:1001-8166.1997.04.006 杜远生, 朱杰, 韩欣, 等.从弧后盆地到前陆盆地-北祁连造山带奥陶纪-泥盆纪的沉积盆地与构造演化[J].地质通报, 2004, 23(9/10):911-917. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=200409163&flag=1 张建新, 许志琴.北祁连中段加里东俯冲-增生杂岩/火山弧带及其变形特征[J].地球学报, 1995, (2):36-46. http://www.cnki.com.cn/Article/CJFDTotal-DQXB502.003.htm 左国朝, 吴汉泉.北祁连中段早古生代双向俯冲-碰撞增生模式剖析[J].地球科学进展, 1997, 12(4):315-323. doi: 10.3321/j.issn:1001-8166.1997.04.002 夏林圻, 夏祖春, 任有祥, 等.祁连秦岭山系海相火山岩[M].武汉:中国地质大学出版社, 1991. SengÖr A M C. The Palaeo-Tethyan suture:A line of demarcation between two fundam entally differen tarch itectural styles in the stru cture of Asia[J]. The Island Arc, 1992, 1:78-91. doi: 10.1111/j.1440-1738.1992.tb00060.x
SengÖr A M C, Natalin B A, Burtman U S. Evolution of the Altaid tectonic collage and Paleozoic crustal growth in Eurasia[J]. Nature, 1993, 364:209-304. http://www.researchgate.net/publication/31960368_Evolution_of_the_Altaid_tectonic_collage_and_Palaeozoic_crustal_growth_in_Eurasia
Windley B F, Alexeiev D, Xiao W, et al. Tectonic models for accretion of the Central Asian Orogenic belt[J]. Journal of Geological Society, 2007, l164:31-47. http://www.researchgate.net/publication/27246576_Tectonic_models_for_accretion_of_the_Central_Asian_Orogenic_Belt
Sengör A M C. East Asia tectonic collage[J]. Nature, 1987, 318:16-17. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0210762630/
Yin A, Harrison T M. Geological evolution of the Him alayanTibetan orogen[J]. Annual Reviews of Earth and Planetary Sciences, 2000, 28:211-280. doi: 10.1146/annurev.earth.28.1.211
Sengör A M C, Natal in B A. Paleotecton ics of Asia: Fragm en tsofa synthesis[C]//Yin A, Harrison M. The tectonic evolution of Asia. Cambridge: Cambridge University Press, 1996: 486-640.
Gehrels G E, Yin A, Wang X F. Magmatic history of the northeasttern Tibetan Plateau[J]. Journal of Geophysical Research, 2003, 108, B9:2423. http://www.researchgate.net/publication/239724885_Magmatic_history_of_the_northeastern_Tibetan_Plateau
Gehrels G E, Yin A, Wang X F. Detrital zircon geoch ronology of the north eastern Tibetan Plateau[J]. Geological Society of America Bulletin, 2003, 115:881-896. doi: 10.1130/0016-7606(2003)115<0881:DGOTNT>2.0.CO;2
董国安, 杨怀仁, 杨宏仪, 等.祁连地块前寒武纪基底锆石SHRIMP U-Pb年代学及其地质意义[J].科学通报, 2007, 52:1572-1585. doi: 10.3321/j.issn:0023-074X.2007.13.015 吴才来, 郜源红, 吴锁平, 等.柴北缘西段花岗岩锆石SHRIMP U-Pb定年及其岩石地球化学特征[J].中国科学:地球科学, 2008, 38(8):930-949. doi: 10.4037-ccn2010235/ 吴才来, 姚尚志, 杨经绥, 等.北祁连洋早古生代双向俯冲的花岗岩证据[J].中国地质, 2006, 33(6):1196-1208. http://d.old.wanfangdata.com.cn/Periodical/zgdizhi200606002 吴才来, 徐学义, 高前明, 等.北祁连早古生代花岗质岩浆作用及构造演化[J].岩石学报, 2010, 26(4):1027-1044. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201004003 Song S G, Niu Y L, Su L, et al. Tectonics of the north Qilian orogen, NW China[J]. Gondwana Research, 2013, 23(4):1378-1401. doi: 10.1016/j.gr.2012.02.004
Song S G, Niu Y L, Zhang L R, et al. Tectonic evolution of early Paleozoic HP metamorphic rocks in the north Qilian mountains, NW China:New perspectives[J]. Journal of Asian Earth Sciences, 2009, 35:334-353. doi: 10.1016/j.jseaes.2008.11.005
张建新, 孟繁聪.北祁连和北阿尔金硬柱石榴辉岩:冷洋壳俯冲作用的证据[J].科学通报, 2006, 51(14):1683-1688. doi: 10.3321/j.issn:0023-074X.2006.14.011 Le Martre R W, Bateman P, Dudek A, et al. A classification of igneous rocks and glossary of terms[C]//Recommendation of the IUGS subcommission on the systematics of igneous rocks. Oxford: Blackwell, 1989.
Winchester J A, Floyd P A. Geochemical discrimination of different magma series and their differentiation products using immobile elements[J].Chemical Geology, 1977, 20:325-343. doi: 10.1016/0009-2541(77)90057-2
Pearce J A, Harris N, Tindle A G. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks[J]. Journal of Petrology, 1984, 25(4):956-983. doi: 10.1093/petrology/25.4.956
张建新, 孟繁聪, Mattinson C G.南阿尔金-柴北缘HP/UHP变质带研究进展、问题及挑战[J].高校地质学报, 2007, 13(3):526-545. doi: 10.3969/j.issn.1006-7493.2007.03.021 董顺利, 李忠, 高剑, 等.阿尔金-祁连-昆仑造山带早古生代构造格架及结晶岩年代研究进展[J].地质论评, 2013, 59(4):731-747. doi: 10.3969/j.issn.0371-5736.2013.04.012 张旗, 孙晓猛, 周德进.祁连蛇绿岩的特征、形成环境及构造意义[J].地球科学进展, 1997, 12(4):366-393. doi: 10.3321/j.issn:1001-8166.1997.04.007 宋述光.北祁连山古大洋俯冲带高压变质岩研究评述[J].地质通报, 2009, 28(12):1769-1778. doi: 10.3969/j.issn.1671-2552.2009.12.010 王永和, 焦养泉, 李建星, 等.中祁连北缘奥陶纪岩浆弧地层[J].现代地质, 2008, 22(5):724-732. doi: 10.3969/j.issn.1000-8527.2008.05.006 黄增保, 郑建平.中祁连西段晚寒武世埃达克岩的发现及地质意义[J].岩石矿物学杂志, 2014, 33(6):1008-1018. doi: 10.3969/j.issn.1000-6524.2014.06.002 冯益民, 何世平.祁连山大地构造与造山作用[M].北京:地质出版社, 1996. Feng Y M, He S P. Orogenic process of the Qilian Mountains[J]. Acta Geosci Sinica, 1996, Special Issue:1-5. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=987984de2313aeb27adaa7fc484111ff
Wu H Q, Feng Y M, Song S G. Metamorphic deformation of blueschist belts and their tactonic implications in North Qilian Mountains[J]. China Metamorphic Geol., 1993, 11:523-536. doi: 10.1111/j.1525-1314.1993.tb00169.x
冯益民.祁连造山带研究概况——历史、现状及展望[J].地球科学进展, 1997, 12(4):307-314. doi: 10.3321/j.issn:1001-8166.1997.04.001
计量
- 文章访问数: 4299
- HTML全文浏览量: 365
- PDF下载量: 2384
