弧形构造带形成演化特征研究

马海勇　罗安湘　王朝阳　欧阳征健　冯娟萍

摘要：為研究弧形构造带形成演化特征，以六盘山弧形构造带早古生代晚期形成的弧形构造为例，采用野外面地质调查、数值模拟等研究手段，对弧形构造带形成规律进行区域块体应力定量化的数值模拟分析。结果表明：六盘山弧形构造带受到青藏地块向北东向的推挤，北祁连地块也向北东向位移并推挤北祁连弧后盆地，北祁连弧后盆地被相对坚硬的阿拉善地块与鄂尔多斯地块阻挡。早古生代晚期六盘山区域北东向挤压的构造应力是早期弧形构造形成的主要原因，利用数值模拟方法构建的古生代晚期的数值模型，进行构造应力与地质块体形变规律数值模拟计算，使地质模型处于平衡状态后获得x方向、y方向位移，在x方向的位移，计算结果显示六盘山区域范围内古生代地质块体向东最大位移约为向北最大位移的25%，位移向量的位移方向总体表现为向北东弧形突起的构造形态。研究结果可为区域弧形构造演化研究提供定量化分析的参考和理论依据。

关键词：弧形构造;六盘山;构造演化;数值模拟;位移向量

中图分类号：TD 121文献标志码：A

文章编号：1672-9315（2022）04-0752-08

DOI：10.13800/j.cnki.xakjdxxb.2022.0415

Formation and evolution characteristics of

arcuate tectonic belt MA Haiyong LUO Anxiang WANG Chaoyang OUYANG Zhengjian FENG Juanping

（1.College of Geology and Environment，Xian University of Science and Technology，Xian 710054，China;

2.Research Institute of Petroleum Exploration and Development，Changqing Oilfield Company，Xian 710018，China;

3.Bailie School of Petroleum Engineering，Lanzhou City University，Lanzhou 730070，China）Abstract：In order to study the formation and evolution characteristics of arcuate tectonic belt，the arcuate structure formed in the Late Early Paleozoic in liupanshan arcuate tectonic belt as an example.By means of field profile geological survey and numerical simulation，an analysis was made of the formation rule of regional block stress in the initial stage of arc tectonic belt formation.The results show that the tectonic belt was pushed NE by the Qinghai-Tibet block，and the North Qilian blocks also moved NE and pushed the North Qilian back-arc basin，which was blocked by the relatively hard Alxa block and Ordos block.The tectonic stress action of the NE compression in the study area in the late Early Paleozoic is the main reason for the formation of the early arc-shaped structure.The numerical model of the tectonic belt in the late Paleozoic was constructed via the numerical simulation method to examine the tectonic stress and the deformation law of the geological block.After numerical simulation，the model was in equilibrium with the displacement in x direction and y direction  obtained.The displacement calculation results in x direction show that the maximum eastward displacement of the Paleozoic geological block in Liupanshan region is about 25% of the maximum northward displacement，and the displacement direction of the displacement vector is generally presented as the structural shape of arc-shaped protrusion in northeast direction.The results can provide a  theoretical guidance for quantitative analysis of regional arc tectonic evolution.

Key words：arcuate structure;Liupanshan;tectonic evolution;numerical simulation;displacement component

0引言

六盘山弧形构造带位于青藏高原东北缘，处于阿拉善地块、青藏地块和鄂尔多斯地块的交接部位（图1），特殊的地理位置使其经历了多期次构造变形和复杂的构造演化[1-4]。张培震认为六盘山弧形构造带受青藏地块向北东向的推挤导致中、北祁连地块也向北东向位移并推挤北祁连弧后盆地，北祁连弧后盆地被相对坚硬的阿拉善地块与鄂尔多斯地块阻挡的构造转换带[5]。王东等通过几何学、运动学特征研究，认为六盘山弧形构造带目前有7条近南北向的断层，均呈现为北东向“弓”形态[6-8]。徐化超等通过动力学研究发现六盘山弧形构造带仅新生代就经历了上新世晚期到第四纪早期、中更新世和晚更新世3个期次的构造变形和隆升，构造带相对基底垂向抬升1 500 m[9-10]。近年来许多学者利用几何学、地壳电性、磁性地层年代、GPS测量等方法对六盘山弧形构造带弧形的形成时间和机制进行了大量的研究，并在中、新生代的弧形断层走向、地壳结构及应力形变等方面取得了诸多成果[11-14]，但对于弧形构造形成初期古生代的演化过程研究相对较少[15-16]。向宏发等通过古地震方法分析六盘山区断裂特征[17]，李天斌等通过运动学、动力学研究弧形推覆构造带形成机理[18-19]，张进等通过宁南弧形构造动力分析认为早古生代晚期北东向挤压是早期弧形构造形成的主要原因[20-21]，但是都缺乏地质模型的定量数值分析的证据。张东宁等利用大陆岩石块体动力学数值模型分析应力场、位移场变化[22-25]，因此，在对六盘山古生代弧形构造演化特征系统分析的基础上，利用数值模拟方法构建了古生代晚期构造带的数值模型，进行构造应力与地质块体形变规律研究，为弧形构造带的大陆动力学定量解释提供理论和实践依据。

1区域地质特征

六盘山弧形构造带位于鄂尔多斯盆地西缘、青藏高原东北缘，夹持在鄂尔多斯地块、阿拉善地块、秦祁昆碰撞带之间。六盘山弧形构造带、其南西部位的秦祁昆碰撞带、其北东部位的贺兰山褶皱与其东部的鄂尔多斯盆地西缘冲断带都受到青藏地块隆升的影响，地质现象，构造复杂。六盘山弧形构造带发育的断层从北向南看，走向依次由北西逐渐变为南北，总体呈现一个镜像的“S”形状，北部的断层形态展布较为舒展，南部的断层形态展布变得收束，弧形展布是该区的显著特点。研究区弧形构造的基底是形成于中、晚元古代的秦—祁—贺三叉形裂谷的发育时期到志留纪、泥盆纪前陆盆地发育时期，盖层经历了石炭纪-三叠纪、侏罗纪-白垩纪、新近纪等不同时代的发育，形成现今的构造形态。

六盘山弧形构造带在早古生代晚期受到青藏地块推挤而产生北祁连造山运动的力学作用。由于六盘山区域被阿拉善地块、青藏地块和鄂尔多斯地块包围，在研究区东北部形成银川地堑;北祁连弧后盆地在中、北祁连地块的推挤下逐渐闭合。东部过渡带在逆冲推覆作用下形成许多逆冲推覆体，推覆体主要由下古生界的地层组成。随着中、北祁连地块继续向北推进，六盘山弧形构造带继续向北东方向推移，并且其一部分缓慢汇入银川地堑。对于六盘山弧形构造带的地质边界条件是北部坚硬的阿拉善地块和东部坚硬的鄂尔多斯地块。西部同心等地区的下古生界地层厚度很大，东部中卫、中宁、固原等地区下古生界地层相对较薄，区域整体东西方向的地层厚度存在较大变化。六盘山弧形构造带的边界形状逐渐发生了变动，诸多原因导致整个六盘山弧形构造带地层从上往下弯曲。这些现象正是早古生代晚期处于北祁连弧后盆地东部的六盘山弧形构造带发育的几个重要条件。六盘山弧形构造带古生代形成了初期的弧形构造之后，经历了晚三叠世印支运动、早侏罗世燕山运动Ⅰ、中侏罗世燕山运动Ⅱ、晚侏罗世燕山运动Ⅲ、早白垩世燕山运动Ⅳ、晚白垩世燕山运动Ⅴ、古近纪早喜山运动与新近纪晚喜山运动等8次构造运动，构造带主体向北推移[22-25]。

六盘山弧形构造带主要发育3个断裂带，分别为海原-六盘山断裂、中衛-同心断裂带、牛首山-罗山断裂带。海原-六盘山断裂带在靖远、固原和南华山向西南方向倾斜，倾角为45°～80°。断层面上大量的近乎水平的条纹和台阶，这表明了断层的反向正弦走滑运动。根据海原-六盘山断裂带发育的一系列拉脱盆地，可以确定该断裂带在中新世晚期表现为左旋走滑运动。中卫-同心断裂带在金山地区向西南方向倾斜，倾角在63°～70°之间。中卫-同心断裂属于左旋走滑断层。牛首山-罗山断裂带位置在青铜峡和牛首山这2个地区的西南方向，向西南倾斜，倾角在62°～85°之间。在该区域岩石形变大，根据横切关系显示出2个明显的构造变形期。牛首山-罗山断裂带是在渐新世之前形成的，为反向的左旋走滑断层。

六盘山西部有许多断层和褶皱，大多数是逆走滑构造。褶皱的轴面向西倾斜，说明它们受到了东西向的压应力作用。六盘山地区构造应力整体为北东东向压应力。研究区出露地层为白垩纪和新近地层，晚白垩世地层逆冲于新近地层之上，相对年龄确定主要构造活动出现在新近纪以后（图2）。对比分析和应力分析发现，六盘山西部构造活动强烈，东部构造变形远弱于西部，地层平缓。可以推断，青藏高原隆升的远程效应空间仅限于六盘山西部。

2弧形构造形成机理

六盘山弧形构造带的野外勘查发现早古生代地层已有大范围的弧形断层与褶皱发育。下古生界地层存在不整合，在二维平面上也能观察到地层形态有一定的弧形展布特征。六盘山及其邻区早古生代晚期的弧型构造，西部可宽达200 km左右，而东部可窄至40 km左右，弧型构造的西部是散布的，而东部向南则逐渐收紧聚集，弧形构造的规模比较大，且不属于对称的弧形构造。六盘山古生代弧型构造带的构造应力主要从南东指向北西方向，即来自规模和强度都较大的北祁连加里东造山带。

从平面分布上看，六盘山及邻区弧形构造西翼较宽，南起香山、天井山地区，北至香根-达赖地区北部，南北方向约200 km。弧形构造东翼较窄，牛首山-达赖地区，小罗山-青龙山南北向，东西方向约40～50 km。

六盘山弧形构造带古弧构造的形成是早古生代晚期中、北祁连地块北东向挤压对鄂尔多斯地块西缘和阿拉善地块南部的影响。鄂尔多斯地块与阿拉善地块的边界分布非常特殊，在北东方向上形成了贺兰拉张区;北祁连弧后盆地在中、北祁连地块的挤压下逐渐关闭，形成了一系列以下古生界地层为主的逆冲推覆构造。中、北祁连地块的逆冲楔继续向北东方向挤压并进入贺兰拉张区，由于在逆冲推覆构造带的走向方向上推进时遇到的地层厚度不同，导致逆冲推覆构造带地体变形，在平面上形成了弧形的构造形态。

中生带在古弧形构造的基础上形成了现今六盘山弧形构造带的主体部分。古太平洋板块俯冲作用下，使得六盘山地区总体处于NW-NE向左旋挤压的构造应力场环境，三叠纪的印支运动使该地区隆起，侏罗纪受到有限的拉张，早白垩世早期拉张断陷，早白垩世晚期挤压反转，陆内变形巨大，多个时期多次发生逆冲推覆现象。早白垩世时期弧形构造处于早期走滑拉张沉降阶段，构造带在强烈的挤压作用下大幅抬升，处于短暂的应力松弛阶段，在近东西方向的微弱的应力松弛产生的张力作用下，导致侏罗纪己经形成的正断层再次活动，具有走滑的性质。东西向拉张作用继续加大，六盘山弧形构造带内的正断层开始变的活跃起来，六盘山盆地西侧的边界断层由早期的左行走滑转为右行走滑，六盘山盆地东侧存在断裂带，在六盘山盆地南部整体遭受北东向压应力后，六盘山弧形构造带随六盘山盆地发生了一定程度的逆时针旋转。早白垩世燕山运动晚期弧形构造带处于晚期湖盆萎缩消亡阶段，六盘山盆地继续受到北东向的挤压应力的作用，六盘山盆地地层开始整体发生缓慢抬升。

3弧形构造数值模拟

3.1平面模型的建立

构建六盘山弧形构造带古生代弧形构造形成初期的地质模型。地质模型忽略了弧形构造垂向隆升的变化及其对弧形构造形成的次要影响因素（图4）。地质模型使用二维有限差分程序FLAC软件建立平面数值模型，将数值模拟研究的重点主要集中在确认早古生代晚期形成的六盘山弧形构造带弧形的形态形成的力学机制。

在平面分布中明确构造带弧形的弯折的力学条件是中、北祁连地块向北东挤压。利用FLAC2D构建地质模型，考虑到数值模拟软件基本计算单元网格划分，按照地质模型将区域地块划分为阿拉善地块、鄂尔多斯地块、中北祁连地块、北祁连弧后盆地等地质块体，得到了六盘山弧形构造带及其邻区的平面模型图。平面数值模型按1∶5 000的比例建立模型，最终模型分别代表东西向1 000 km，南北向760 km的研究范围。弧形构造的形成作为该模型的模拟分析重点。数值模拟的主要目的是解释其弯折的力学条件。研究区位于数值模型中心偏北东部，与模型4个边界的距离都较大，主要是为了抵消边界效应对研究区的影响。此外，贺兰拉张区的砥柱作用没有阿拉善地块和鄂尔多斯地块强烈，故将其并入北祁连弧后前陆盆地的一部分进行模拟（图5）。

3.2数值模型初始条件和边界

构建的数值模型北部为阿拉善地块，东部为鄂尔多斯地块，南部为中、北祁连地块，中部为北祁连弧后前陆盆地，由于六盘山弧形构造带形成于北祁连弧后盆地区域内部，因此可以利用代表北祁连弧后前陆盆地的网格形变来反映六盘山弧形构造的形成过程。根据区域内地块的形状特征，利用FLAC2D建立可计算的数值模型并生成计算单元网格。为了在构建的模型的北东向压力作用下网格能产生充分形变且划分尽量少的区块以减少计算负荷，地质块体划分左右边界以南北向为主，整个模型划分成子区块，通过子区块建立完整网格。模型中阿拉善地块、鄂尔多斯地块、北祁连弧后盆地、中北祁连地块分别由不同的网格组成。

由于六盘山弧形构造带北部受阿拉善地块的阻挡，东部受鄂尔多斯地块的阻挡，将模型边界设置为南部边界施加了均等应力;模型的东部和西部为滚动边界，即在东西边界处物质不发生向东和向西运动。模型北部边界为固定边界，南部边界的物质运移是自由的。南部边界的压力被设置为2.0×107 Pa。由于各向同性弹性介質的特点是受到应力作用时应力与应变呈现正相关的线性关系，应力均匀增大时应变也按一定的比率线性增大。本次数值模型将地质块体设定为各向同性弹性介质来建立模型（表1），明确六盘山弧形构造带形成的主要动力学过程。

3.3数值模型计算结果

六盘山弧形构造带的数值模型进行模拟计算后模型处于平衡状态，并得到了位移图。导出x方向、y方向位移等值线图，位移均为正值。对于y方向位移等值线图，数值模型均为正值，即相对于初始阶段的模型各网格节点向北位移，并且南部位移大于北部位移（图6）。数值模型在x方向的位移计算结果表明六盘山区域范围内古生代地质块体位移的主要趋势是向东，向东最大位移约为向北最大位移的25%（图7）。将数值模拟计算结果中x方向位移与y方向位移合并后得到位移场图，从位移向量图指示的位移方向总体表现为向北东弧形突出（图8）。

4结论

1）六盘山弧形构造带经历了长期复杂的构造演化过程，其主体形成于印度-欧亚板块碰撞所导致的青藏高原隆升对其东北缘产生的北东向的挤压应力。六盘山及其邻区早古生代晚期北东向挤压变形的相对柔软的地块遇到坚硬的阿拉善地块、鄂尔多斯地块阻挡的共同作用下经历了早期的弧形构造。数值模拟研究古生代弧形构造形成初期的地球动力学条件对弧形构造的作用机制。

2）依据前人对六盘山及其邻区早古生代晚期弧形构造研究成果构建了数值模型。六盘山弧形构造带南部边界不同时期受力大小与方向是存在变化的，其压力大小虽然变化较大，但压力方向为小幅度变化。数值模型将地质块体做为各向同性弹性介质，并对不同地质块体赋予岩石力学参数，数值模型南部边界施加的北东向应力为均等定值定向压力。

3）数值模型模拟计算的位移结果与弧形的总趋势相一致。数值模拟过程有效验证了北东向挤压作用是六盘山弧形构造带古生代弧形构造形成主要因素。

参考文献（References）：

[1]万天丰.论六盘山——贺兰山碰撞带[J].大地构造与成矿学，2020，44（5）：845-851.WAN Tianfeng.Discussion on the Liupanshan-Helanshan collision zone[J].Geotectonica et Metallogenia，2020，44（5）：845-851.

[2]周鼎武，赵重远，李银德，等.鄂尔多斯西南缘地质特征及其与秦岭造山带的关系[M].北京：地质出版社，1994.

[3]张培震.六盘山地区的构造特征及地壳缩短量[M].北京：地震出版社，1990.

[4]刘凤山，吴中海，张岳桥，等.青藏高原东缘新构造与活动构造研究新进展及展望[J].地质通报，2014，33（4）：403-418.LIU Fengshan，WU Zhonghai，ZHANG Yueqiao，et al.New progress and prospects of neotectonics and active tectonics synthetically study on eastern edge Qinghai-Xizang plateau[J].Geological Bulletin of China，2014，33（4）：403-418.

[5]张培震，邓起东，张国民，等.中国大陆的强震活动与活动地块[J].中国科学D辑（地球科学），2003，33（增刊）：12-20.ZHANG Peizhen，DENG Qidong，ZHANG Guomin，et al.Active tectonic blocks and strong earth quakes in continent of China[J].Science in China（Series D），2003，33（S）：12-20.

[6]王東，张快乐，李伟.六盘山弧形构造带几何学、运动学及年代学研究[J].西部资源，2013（6）：84-86.WANG Dong，ZHANG Kuaile，LI Wei.Geometry，kinematics and chronology of Liupanshan arc-shaped structural belt[J].Western Resources，2013（6）：84-86.

[7]张韬，林畅松，杜玉良，等.六盘山弧形构造带铜金属矿床的成矿环境[J].地球学报，2002，23（6）：553-558.ZHANG Tiao，LIN Changsong，DU Yuliang，et al.Mineralization environments of the Cu polymetallic deposits in the Liupanshan arc structure belt[J].Acta Geoscientia Sinica，2002，23（6）：553-558.

[8]李强，江在森，武艳强，等.海原-六盘山断裂带现今构造变形特征[J].大地测量与地球动力学，2013，33（2）：18-22.LI Qiang，JIANG Zaisen，WU Yanqiang，et al.Present-day tectonic deformation characteristics of HaiYuan-Liupanshan fault zone[J].Journal of Geodesy and Geodynamics，2013，33（2）：18-22.

[9]王双绪，蒋锋云，郝明，等.青藏高原东缘现今三维地壳运动特征研究[J].地球物理学报，2013，56（10）：3334-3345.WANG Shuangxu，JIANG Fengyun，HAO Ming，et al.Investigation of features of present 3D crustal movement in eastern edge of Tibet plateau[J].Chinese Journal of Geophysics，2013，56（10）：3334-3345.

[10]徐化超，王辉，曹建玲.青藏高原东北缘主要断裂滑动速率及其动力学意义[J].地震，2018，38（3）：13-23.XU Huachao，WANG Hui，CAO Jianling.Slip rates of the major faults in the northeastern Tibetan plateau and their geodynamic implications[J].Earthquake，2018，38（3）：13-23.

[11]李文辉，高锐，王海燕，等.六盘山断裂带及其邻区地壳结构[J].地球物理学报，2017，60（6）：2265-2278.LI Wenhui，GAO Rui，WANG Haiyan，et al.Crustal structure beneath the Liupanshan fault zone and adjacent regions[J].Chinese Journal of Geophysics，2017，60（6）：2265-2278.

[12]张晓亮，师昭梦，蒋锋云，等.海原-六盘山弧型断裂及其附近最新构造变形演化分析[J].大地测量与地球动力学，2011，31（3）：20-24.ZHANG Xiaoliang，SHI Zhaomeng，JIANG Fengyun，et al.Research on late tectonic deformation evolvement of Huaiyuan-Liupanshan arc fault and its surrounding area[J].Journal of Geodesy and Geodynamics，2011，31（3）：20-24.

[13]LI Wei，DONG Yunpeng，GUO Anlin，et al.Chronology and tectonic significance of Cenozoic faults in the Liupanshan arcuate tectonic belt at the northeastern margin of the Qinghai-Tibet plateau[J].Journal of Asian Earth Sciences，2013，73：103-113.

[14]戴洪宝，许继影.海原断裂东段至六盘山断裂西段GPS剖面地壳变形与应变积累分析[J].大地测量与地球动力学，2016，36（4）：343-345.DAI Hongbao，XU Jiying.Analysis of the GPS section deformation and strain accumulation of the eastern Haiyuan and Western Liupanshan faults[J].Journal of Geodesy and Geodynamics，2016，36（4）：343-345.

[15]陈虹，胡健民，渠洪杰，等.川滇南北向构造带早中生代构造变形研究[J].中国科学D辑（地球科学），2011，41（9）：1281-1294.CHEN Hong，HU Jianmin，QU Hongjie，et al.Early Mesozoic structural deformation in the Chuandian N-S tectonic belt，China[J].Science China Earth Sciences，2011，41（9）：1281-1294.

[16]张岳桥，廖昌珍，施玮，等.鄂尔多斯盆地周边断裂运动学分析与晚中生代构造应力体制转换[J].地质学报，2006，80（5）：639-647.ZHANG Yueqiao，LIAO Changzhen，SHI Wei，et al.Fault kinematic analysis and change in late Mesozoic tectonic stress regimes in the peripheral zones of the Ordos Basin，North China[J].Acta Geologica Sinica，2006，80（5）：639-647.

[17]向宏發，池田安隆，张晚霞，等.六盘山东麓的古地震研究[J].中国地震，1999，15（1）：74-81.XIANG Hongfa，YASUTAKA Ikeda，ZHANG Wanxia，et al.Study on paleoearthquakes of the eastern Liupanshan piedmont fault zone[J].Earthquake Research in China，1999，15（1）：74-81.

[18]李天斌.宁夏南部弧形推覆构造带特征及演化[J].地质力学学报，1999，5（3）：22-27.LI Tianbin.The feature and formation mechanism of arciform thrust-nappe structure zone of the southern NingXia[J].Journal of Geomechanics，1999，5（3）：22-27.

[19]李天斌.鄂尔多斯盆地西缘逆冲推覆构造特征及演化[D].北京：中国地质大学（北京），2006.LI Tianbin.The characteristics and evolution of thrust nappe structure in the west margin of Ordos basin[D].Beijing：China University of Geoseienees（Beijing），2006.

[20]张进，马宗晋，任文军，等.宁夏中南部古生代弧型构造[J].大地构造与成矿学，2004（1）：29-37.ZHAGN Jin，MA Zongjin，REN Wenjun，et al.Paleozoic arcuate thrust belt in central and southern NingXia[J].Geotectonica et Metallogenia，2004，28（1）：29-37.

[21]ZHANG Jin，ZHANG Yinping，XIAO Wenxia，et al.Linking the Alxa terrane to the eastern gondwana during the early Paleozoic：constraints from detrital zircon U-Pb ages and cambrian sedimentary records[J].Gondwana Research，2015，28（3）：1168-1182.

[22]陈连旺，杨树新，谢富仁，等.中国大陆构造应力应变场现今年变化特征的数值模拟[J].中国地震，2005，21（3）：341-349.CHEN Lianwang，YANG Shuxin，XIE Furen，et al.Numerical simulation on annual change patterns of present-day tectonic stress-strain field on the Chinese mainland[J].Earthqueake Research in China，2005，21（3）：341-349.

[23]刘鎏，魏东平.中国大陆及邻区板内应力场的数值模拟及动力机制探讨[J].地震学报，2012，34（6）：727-740.LIU Liu，WEI Dongping.Numerical simulation of the intraplate stress field in China mainland and adjacent areas and its dynamic implication[J].Acta Seismologica Sinica，2012，34（6）：727-740.

[24]张东宁，许忠淮.中国大陆岩石层动力学数值模型的边界条件[J].地震学报，1999，21（2）：21-27.ZHANG Dongning，XU Zhonghuai.Boundary conditions of the numerical model of lithosphere dynamics in mainland China[J].Acta Seismologica Sinica.1999，21（2）：21-27.

[25]马德云，高振敏，杨世瑜，等.北衙金矿区构造应力场数值模拟[J].大地构造与成矿学，2003，27（2）：160-166.MA Deyun，GAO Zhenmin，YANG Shiyu，et al.Numerical simulation for tectonic stress field of the Beiya gold deposit[J].Geotectonica et Metallogenia，2003，27（2）：160-166.