斜向俯冲（逆冲）作用，指汇聚-俯冲或者碰撞方向与大型断裂带走向不垂直（即存在走滑分量）的逆冲作用，在板块俯冲带和褶皱冲断带普遍存在。典型的斜向俯冲带有Sumatra^[1]，Nankai^[2-3]，Chilean^[4]等；典型的斜向褶皱冲断带有Himalaya^[5-6]，Alps^[7-8]，Taiwan^[9-10]，Andes^[11-12]，Zagros^[13-14]等。这些地区既是近10a来大地震频发区^{[3, 5, 10, 15-17]}，也是研究板块构造理论的天然实验室。此外，研究斜向逆冲作用可为能源勘探^[18]、矿产资源开发^[19]、地震危险性评估^[15]、地貌演化^[19]等提供科学依据。

龙门山属于非典型性的褶皱冲断带，具有独特的地形地貌和构造特征，主要包括：①陡降地形，作为环青藏高原周缘最陡峭的地形地貌边界，部分地方可以在50km的范围内由5000m降至500m^[20]；②陡变地壳厚度，从青藏高原腹地的60~70km降至四川盆地的40~50km^[21]；③低水平缩短速率，汶川地震之前的GPS测量数据表明水平缩短速率小于3mm/a^[22]；④山前川西盆地缺乏连续完整的新生代沉积记录, 这似乎暗示该地区缺乏新生代以来的前渊凹陷和构造挠曲变形^{[20, 23-25]}。

龙门山地区年轻的高地形和低水平缩短速率、缺乏前渊凹陷之间形成强烈对比，这一特征直接启发了下地壳流动模式，Royden等^[26]认为，下地壳物质从高原内部向东流动过程中受到四川盆地坚硬克拉通的阻挡而产生侧向支撑作用，使龙门山能够在没有大规模逆冲推覆和挠曲沉降的情况下快速隆升，这种地壳增厚主要发生在下地壳，而上地壳很少或没有明显的缩短变形。该模式虽然可以很好地解释龙门山边界存在的上述4个非同寻常的特征^[26-28]，但造成了超过8万人死亡和近10亿元的经济损失的2008年汶川地震，以极端震撼的方式展示了龙门山内部强烈的变形作用和正在进行的青藏高原隆升过程^[29-31]。汶川地震所体现的上地壳缩短作用不容忽视，因此对流行的下地壳流动模式提出了严峻的挑战，同时也提醒人们需重新审视该理论和其他支撑的证据是否存在更合理的解释^[29]。汶川地震的发生，再一次引发了青藏高原东缘的隆升机制究竟是下地壳流动还是上地壳脆性缩短变形的争论^{[26, 28-30, 32-35]}。

争论的焦点在于晚新生代龙门山地壳缩短变形模式。中生代以来龙门山主要经历了晚三叠世和新生代2期重要的地壳缩短变形^{[20, 23-24]}。晚三叠世龙门山在调节SE向逆冲挤压作用下形成NE向延伸的龙门山褶皱冲断带和川西前陆盆地，表明早期造山作用为NW-SE向^{[20, 24]}。晚新生代印度-欧亚大陆碰撞的远距离效应使得龙门山褶皱冲断带被重新激活^{[20, 24]}。但是，对于这一阶段青藏高原东缘的挤压方向存在分歧：一些学者认为晚新生代造山作用方向与晚三叠世一致^[36-37]；另一种观点则指出晚新生代构造变形可能包含重要的右旋走滑分量^{[20, 30]}，后者与汶川地震同震滑移指示的动向一致^{[30-31, 38-39]}（图 1）。但这一斜向逆冲作用的活动时间、运动学和动力学机制缺乏准确的约束。

图  1  青藏高原东缘主要活动断层和大地震分布

a-青藏高原东缘周缘GPS位移场（据参考文献[22]修改）；b-青藏高原东缘龙门山褶皱冲断带和四川盆地西部区域构造位置，2008年汶川地震震中位置和震源机制解，同震地表破裂和未破裂的主断层^{[30, 38-39]}, 1970年大邑地震、1976年松潘地震和2013年芦山地震震源机制解^{[34, 42-43]}和区域地震剖面位置均在图中标明（最大水平主应力数据修改据参考文献[39, 44]）；c-小鱼洞断层地表破裂和同震滑移场（数据修改据参考文献[45-46]）。WMF-汶川-茂汶断裂；YBF-映秀-北川断裂；PGF-彭县-灌县断裂；QCF-青川断裂；MJF-岷江断裂；HYF-虎牙断裂；LRBF-龙日坝断裂；WLF-武龙断裂；BXF-宝兴断裂；DYF-大邑断层；QXF-邛西断层；PLF-平罗坝断层；LQF-龙泉断层；XPF-熊坡断层

Figure  1.  Main active faults and earthquake events on the eastern margin of the Tibetan Plateau

下载: 全尺寸图片 幻灯片

本次研究基于龙门山南段山前大邑背斜区三维地震数据，通过野外地质实测和地震解释结果，分析地下和地表构造变形走向上的差异性；结合地震相干分析建立三维构造模型，明确2期构造的三维特征及其切割关系；结合背斜区已发表的浅层地震剖面揭示的生长地层记录和大邑砾岩沉积年龄限定，确定出东西斜向逆冲作用的准确时间，并进一步探讨区域构造应力场转换时序及其运动学和动力学机制。

1.   地质背景

龙门山褶皱冲断带位于青藏高原东缘与四川盆地之间，南起雅安北至广元，长约500km，宽30~50km，呈NE向延伸，北东端与秦岭断裂带斜交，南西端被鲜水河断裂斜截。该褶皱带由4个叠瓦状逆冲推覆系统构成，自西向东分别为后山断裂（汶川-茂汶断裂）、中央断裂（映秀-北川断裂）、山前断裂（彭县-灌县断裂）及山前隐伏断裂^{[23-24, 40-41]}（图 1）。

龙门山褶皱冲断带在地形地貌、构造变形、构造演化、前陆盆地发育、构造活动性等方面沿走向的差异一直被众多学者所关注。Jia等^{[24, 47]}根据前陆盆地地层发育和构造形态特征，以安县为界将龙门山褶皱冲断带分为南、北2段，北段类似于阿尔卑斯和喀尔巴契的逆冲推覆构造带^[20]，其中晚三叠世变形强于新生代，南段以基底卷入的叠瓦状冲断为特点，新生代变形尤为突出；Li等^[48]通过对比龙门山褶皱冲断带沿走向变形、汶川地震同震破裂和同震断裂的三维构造建模等方面的差异，同样以安县为界将龙门山褶皱冲断带分为南、北2段，北段为以右旋斜向逆冲为主的单条断层且不能向下延伸至主滑脱面，而南段则是以逆冲为主的2条断层，大约在10km深处汇成一条断层向下延伸至主滑脱面。Jin等^[49]结合地形地貌、野外调查和地震剖面解释结果等差异性，以安县和都江堰为界将龙门山褶皱冲断带分为分为南、中和北3段（图 1）。本次研究集中于都江堰以南的龙门山南段。

龙门山南段山前不仅发育一系列突破至近地表的NE向构造，如大邑、平罗、熊坡和龙泉山构造^{[29, 47-48, 50-51]}，还发育NS向隐伏构造，如邛西和大塘构造^[52]（图 1）。其中, 大邑构造位于雾中山推覆构造带前缘，地表出露白垩系红层和晚新生代砾石层，呈宽缓的NE向背斜形态（图 2），Densmore等^[53]结合野外调查和遥感影像解译证实，该构造第四纪以来存在活动，且1970年大邑地震（Ms 6.2）发生于此构造上^[54]（图 2）。

图  2  大邑背斜区大邑和崇州构造及其邻区地质图

地震剖面位置在图中均标明。WMF-汶川-茂汶断裂；WLF-武龙断裂；YBF-映秀-北川断裂；PGF-彭县-灌县断裂；BXF-宝兴断裂；PZF-彭州断层

Figure  2.  Geological maps of the DYFS and adjacent regions

下载: 全尺寸图片 幻灯片

2.   大邑背斜构造解释

2.1   野外调查和卫星影像解译

川西盆地广泛缺失晚新生代的沉积地层, 仅有零星的大邑砾岩、雅安砾岩和上更新统局部出露。在四川盆地内部，大量的沉积岩被大面积剥蚀，形成角度不整合面^{[20, 55]}。在这些不整合面之上，几十至几百米厚的以碎屑支撑的砾岩保留在龙门山山前地区^[53]。李勇等^[25]等利用电子自旋共振对大邑砾岩和雅安砾岩进行了定年，将大邑砾岩和雅安砾石的年龄分别限定为3.6~0.8Ma（上新世-早更新世）和0.6~0.2Ma（中更新世）^{[25, 53, 56-57]}。Kong等^[58]利用宇宙成因核素¹⁰Be和²⁶Al得到大邑砾岩的埋藏年龄为2.2±0.7Ma、2.1±1.4Ma和2.7±1.2Ma，这一沉积年龄与李勇等^[25]的结果一致。大邑砾岩的典型地层剖面位于岷江西岸的大邑氮肥厂附近，该剖面也是大邑砾岩的命名地。大邑构造研究区内广泛出露大邑砾岩，这为大邑构造晚新生代构造变形提供了充分的依据（图 2）。

野外调查表明，大邑县城北剖面中大邑砾岩（产状135°∠51°）与下伏白垩系灌口组红层（146° ∠55°）呈微角度不整合接触（图 3-a~d）；崇州河家山剖面中大邑砾岩（155°∠36°）与下伏灌口组（153° ∠42°）也呈微角度不整合接触（图 3-e）。其他大邑砾岩剖面，如灌县玉堂场汤家沟、崇庆街子场、白塔湖、大邑陈家场白岩沟、团结乡、邛崃水口场及名山万古乡庙坡剖面，均显示与下伏地层呈角度不整合接触，少量呈平行不整合接触。大邑背斜区大邑砾岩也发生褶皱变形，结合野外地质填图、3D扫描仪、全站仪测量等技术手段，发现这一区域的大邑砾岩及其下伏地层的褶皱变形时间应为上新世-早更新世，而最新的活动时间应为中更新世-全新世^{[25, 53, 59-60]}。

图  3  大邑构造区大邑砾岩与下伏地层不整合接触关系

a~d大邑县城北剖面角度不整合；e-大邑县何家山剖面角度不整合

Figure  3.  Structural cross sections across the DYFS, showing that the Dayi conglomerate rests disconformably on the Late Cretaceous red beds

下载: 全尺寸图片 幻灯片

2.2   大邑背斜地震解释

本次选取了6条过大邑背斜区的石油地震剖面A~F，分别经叠后偏移处理和深度转换，揭示其地下构造几何学特征。其中，地层主要依据钻井资料约束的地层标定，而近地表地层由地质图约束（图 2）。地震剖面中断层和褶皱以断面反射、断层切割、褶皱翼或膝折端终止和倾角域特征来识别（图 4-图 9）。

图  4  过崇州构造北翼的地震剖面A解释方案（剖面位置见图 2）

Pre-P-前二叠系；P-二叠系；T₁j-下三叠统嘉陵江组；T₂l-中三叠统雷口坡组；T₃m-上三叠统马鞍塘组；T₃x-上三叠统须家河组；J₁-下侏罗统；J₂-上侏罗统；K-白垩系；E-古近系；N₂-Q₁-大邑砾岩；Q-第四系

Figure  4.  Interpreted seismic profile A across the Chongzhou structure

下载: 全尺寸图片 幻灯片

图  5  过崇州构造南翼的地震剖面B解释方案（剖面位置见图 2，地层代号同图 4）

Figure  5.  Interpreted (b) seismic profile B across the Chongzhou structure

下载: 全尺寸图片 幻灯片

图  6  过大邑构造北翼的地震剖面C解释方案（剖面位置见图 2，地层代号同图 4）

Figure  6.  Interpreted seismic profile C across the Dayi structure

下载: 全尺寸图片 幻灯片

图  7  过大邑构造的地震剖面D解释方案（剖面位置见图 2，地层代号同图 4）

Figure  7.  Interpreted seismic profile D across the Dayi structure

下载: 全尺寸图片 幻灯片

图  8  过大邑构造南翼的地震剖面E解释方案

Figure  8.  Interpreted seismic profile E across the Dayi structure

下载: 全尺寸图片 幻灯片

图  9  沿走向地震剖面F，展示滑脱层厚度在走向上的变化和F5断层的侧断坡（剖面位置见图 2）

Figure  9.  Along-strike NE-trending seismic profile F across the Dayiand Chongzhou structures, showing the along-strike changes of the thickness and lateral ramp of F5 fault

下载: 全尺寸图片 幻灯片

根据地震反射剖面A和B、石油钻井数据及地表地质特征，发现大邑背斜北部的崇州构造主要受东倾F4逆冲断层控制。由F4断层控制的上伏背斜主要呈典型的断层相关褶皱特征：背斜前翼短而陡、后翼长而缓；背斜后翼地层与下伏断坡不平行；后翼地层倾角较断坡角（F4）小。根据这些典型的特征，推测崇州构造存在纯剪变形作用^[61-62]。考虑断层倾角（θ≈59°~64°）和后翼地层倾角（δ_b≈9°~22°），解释这一构造为典型的纯剪断层转折褶皱^[61-62]，大量滑移量自浅部滑脱层传递至断层上，塑性滑脱层系（下三叠统嘉陵江组）吸收和调节少量的构造缩短（图 4、图 5）。实际上，这一构造比纯剪断层相关褶皱端元模型复杂，深部构造对这一构造进行了改造。地震反射剖面A和B中，崇州构造东翼被一东倾的前冲断层切割，这一断层的形成跟剪切作用的调节和褶皱生长中层面偏移的后翼弯曲有关。从地震剖面A可以看出，F5断层切割了崇州构造后翼的大邑砾岩，暗示这一断层形成于大邑砾岩（N₂-Q₁）沉积之后，而崇州构造的形成时间可能对应大邑砾岩的沉积时间（图 4、图 5）。

结合地震反射剖面C、D和E，以及石油钻井数据和地表地质特征对比发现，位于大邑背斜区南部的大邑构造呈复杂的构造几何学特征。大邑构造主要受深部F1断层控制，浅部F2和F3断层对这一构造进行改造，深浅构造由浅部的下三叠统嘉陵江组塑性膏岩层相连；断层倾角（θ≈26°~48°）大于后翼地层倾角（δ_b≈10°~24°），由于深部断层F1浅部缓的断坡使得断层滑移量沿断层向浅部减少；大邑构造被下三叠统嘉陵江组发育的西倾的F2断层和东倾的F3断层切割，这2条断层的形成与剪切作用的调节和褶皱生长过程中层面偏移的后翼弯曲有关，同时这些断层将大量的位移量自滑脱层传播至近地表（图 6-图 9）。

2.3   大邑背斜变形样式和机制走向差异性

褶皱冲断带变形样式和变形机制沿走向存在差异性^[63]，大邑背斜同样存在这些特征。大邑背斜深部，奥陶系和志留系主要由页岩、砂岩和碳酸盐岩组成，而300~900m厚的页岩可能由于其埋藏较深，经历了高温高压而成为这一区域的中部滑脱层系；浅部早三叠统嘉陵江组由碳酸盐岩、砂岩和蒸发岩组成，400~1100m厚的蒸发岩成为这一区域的浅部滑脱层系。沿走向地震剖面F（图 9）和垂直走向地震剖面A~E（图 4-图 8）均揭示浅部滑脱层由南向北逐渐加厚。为了讨论大邑背斜变形样式和变形机制沿走向的差异性，本文基于地震反射剖面建立了3个以断层相关褶皱理论为基础的构造几何学和运动学模型，包括一个单层滑脱层模型和2个底部滑脱层厚度有差异的双层滑脱层模型，以此讨论滑脱层的存在与否及其厚度变化对构造变形样式的影响（图 10）。对比单层滑脱层和双层滑脱层模型，发现浅部滑脱层存在与否对构造几何学和构造运动学特征有重要影响（图 10），而浅部滑脱层对构造变形作用的影响在物理模拟^[63-64]和数值模拟^[65]中均有讨论，在此不做深入探讨。

图  10  大邑背斜构造运动学模型

a-上部滑脱层不存在的构造运动学模型，以多膝折断层转折褶皱变形为主；b-大邑背斜北部上部和中部滑脱层均存在的构造运动学模型，以塑性变形的滑脱褶皱变形为主；c-大邑背斜南部上部和中部滑脱层均存在的构造运动学模型，以脆性变形的断层转折褶皱变形为主

Figure  10.  Conceptual kinematic models of folding in the DYFS

下载: 全尺寸图片 幻灯片

结合典型的构造特征，本文建立了以双滑脱层为基础的概念模型，以塑性变形为主的代表大邑背斜北部的崇州构造（图 10-b），以脆性变形为主的代表大邑背斜南部的大邑构造（图 10-c）。图 10-b中1展示崇州构造北部变形前的几何学特征，在变形的第一阶段，构造缩短被中间滑脱层的滑脱褶皱变形样式的加厚形式所吸收（图 10-b中2）；连续的变形导致浅部滑脱层系内的塑性流体向着背斜核部流动而使之上的沉积岩产生不对称的弯曲（图 10-b中3）；最后，背斜前后翼发育断层切割能干层，向后的逆冲断层发育在后翼，而向前的逆冲断层发育在前翼，这些构造均连接至浅部滑脱层且改造了深部的剪切断层转折褶皱（图 10-b中4）。对比大邑构造地震剖面，深部断层F1由2个典型的膝折带组成，由此建立由向斜型膝折带和背斜型膝折带组成的深部断层（图 10-c）。图 10-c中1展示变形前的构造几何学特征；在变形的第一阶段，构造缩短被以断层转折褶皱变形样式的断层上盘褶皱形成所吸收，致使构造隆升，由断层作用形成的轴面控制上覆与其平行的地层倾角^{[62, 66]}（图 10-c中2）；连续的断层滑移使得背斜膝折带进一步加宽，致使构造进一步隆升（图 10-c中3）；最后，背斜前后翼均发育断层切割能干层，这些构造均连接至浅部滑脱层且改造了深部的剪切断层转折褶皱（图 10-b中4）。

3.   大邑背斜三维构造特征

3.1   基于三维地震相干刻画断层几何学特征

Bahorich等^[67-68]最早通过三维地震反射数据提取的地震相干切片属性来刻画断层构造几何学特征和地震道之间的不连续现象，如河道沉积、不整合等。地震相干技术是1995年第65届SEG年会上推出的一种解释性处理技术，可以分辨三维地震数据体中细微断层或其他（地震）道-道之间的不连续^[67]，其出现被视为三维地震领域里的革命。通过计算在主测线和联络测线中局部波形的相似性，可以获得三维地震相干体^[67]，从而可在相干切片上很清楚地识别出断层、不同的岩性体系特征。本研究利用Landmark软件中SeisWorks三维模块来制作大邑背斜三维地震数据相干体，消除人为因素定义的大邑背斜地下断层几何学特征，包含0~1数值范围的数据体，其中0代表相邻地震道之间低相干性和低相似性，一般代表断层位置；而1代表相邻地震道之间的高相干性，两侧无差别。为了探讨大邑背斜浅部断层间的相互关系，制作了0~6km深的三维地震相干数据体，考虑到浅部断层出露深度情况，选取了3个不同深度的相干切片，分别为2.8km、3.3km、4.0km深（图 11）。根据地震相干切片中相干性的大小，提取出大邑背斜区主要的断层F2~F6，这些低相干性、低相似性的位置（断层）（图 11）和地震反射剖面中解释的断层位置一致（图 4-图 9），很好地验证了前文中的构造解释方案。地震相干切片清晰地展示了存在2组不同走向的断层，NNE向和NE走向，东倾的逆冲断层F2和F6为NE走向（N47°E-N70°E），西倾的逆冲断层F3和F4为NNE走向（N15°E-N25°E），东倾的逆冲断层F5为NNE走向（N9° E-N13° E）（图 11）。NNE走向断层F3切割了NE走向断层F6和北西侧的背斜，暗示2组走向断层形成时间上存在先后关系，NNE走向断层F3~F5发育时间晚于NE走向断层F2和F6（图 11）。

图  11  未解释（左）和解释的（右）大邑背斜区三维地震相干深度切片（深度分别为2.8km、3.3km和4.0km）

Figure  11.  Uninterpreted (left) and interpreted (right) depth coherence slices at depths of 2.8km, 3.3km, and 4.0km from 3D seismic data in the DYSF

下载: 全尺寸图片 幻灯片

3.2   三维构造建模刻画断层几何学特征

基于大邑背斜三维地震数据构造解释成果，结合相干分析和钻井数据，本文在GOCAD地质三维建模平台上构建了大邑构造三维构造模型^[69]，以便通过三角面网格识别三维断层构造几何学特征，确认断层之间、断层和背斜或向斜间的空间关系，评价断层几何学沿走向上的变化，讨论深部和浅部构造关系等地下三维构造问题。将Discovery或Land-mark中地震解释模块解释的地层数据和断层数据以离散点方式导出，在GOCAD地质三维建模平台导入这些点数据，使用DSI（离散元插值）^[69]算法分别建立断层面模型和地层模型，确立层位和断层、断层和断层的空间切割关系，实现地下三维可视化，该三维构造模型中的断层系统是现今大范围地震灾害的研究基础^[70]。首先进行地层点的拾取（图 12-a），其次对地层和断层进行拾取（图 12-b），然后基于断层离散点和断层切割关系建立断面模型（图 12-c），最后基于地层离散点，以及断层和地层的接触关系建立地层模型（图 12-d），从而建立大邑背斜三维构造模型，主要包括7条主要断层和6个地层底界面（分别为二叠系、下三叠统嘉陵江组、中三叠统雷口坡组、上三叠统须家河组、下侏罗统自流井组和上侏罗统遂宁组）（图 12）。在该三维构造模型中，未加入地震剖面中识别的走向延伸小于3km的次级小断层，例如图 4和图 5中左上侧的断层。

图  12  大邑背斜三维构造模型

Figure  12.  3D structural modeling of the DYFS

下载: 全尺寸图片 幻灯片

大邑背斜中的大邑构造和崇州构造几何学特征在三维构造模型中得以清楚展示。在图 12-d中，由下三叠统嘉陵江组和中三叠统雷口坡组限定的浅部滑脱层系被突出，而蒸发膏岩层在大邑构造较薄，在崇州构造区较厚。大邑构造和崇州构造浅部构造主要受浅部滑脱层系（4~7km深）上部的隐伏逆冲断层改造和控制，而这一滑脱层可向西延伸与龙门山褶皱冲断带深部滑脱层相连。大邑构造主要受深部隐伏的断层转折褶皱控制，其断坡连接了中部9~10km深的志留系滑脱层和浅部4~8km深的下三叠统滑脱层，中部滑脱层也可向西延伸，与龙门山褶皱冲断带深部滑脱层相连。大邑背斜三维构造模型同样展示了2组不同走向的断层F1~F6，包括NE走向的逆冲断层F1、F2和F6，NNE走向的逆冲断层F3、F4和F5（图 12）。

图 13展示了大邑背斜区下三叠统嘉陵江组和上三叠统须家河组底界的地下三维构造模型。下三叠统嘉陵江组三维构造模型揭示，大邑背斜区的大邑构造被深部NE走向的断层F1切割，而崇州构造被另一条深部NE走向的断层切割（图 13-a）。大邑背斜区大邑构造发育西倾的逆冲断层F2和F6，以及东倾的反冲断层F3，而崇州构造发育西倾的逆冲断层F5和东倾的反冲断层，中三叠统、上三叠统、侏罗系和白垩系均被这些断层切割。另外，大邑构造西翼被NNE走向的逆冲断层F3切割，大邑构造背斜区和向斜区也被该断层切割；崇州构造被2条NNE走向的逆冲断层F4和F5限制，断层F4控制NNE走向的背斜区，而断层F5同样切割了大邑构造向斜区西翼（图 13-b）。断层和断层切割关系，以及与相应背斜、向斜的切割关系暗示，大邑背斜区存在2期构造事件。结合地震剖面解释、地震相干分析和三维构造建模，可以确定大邑背斜区发育的NNE向构造晚于NE向构造，而NNE向构造的发育与近EW向的构造缩短作用相关。

图  13  大邑背斜区下三叠统嘉陵江组（T1j）（a）和上三叠统须家河组(T3x)（b）底界三维构造模型

Figure  13.  Contour map of the bottom of the Lower Triassic Jialingjiang (T₁j) (a) and the Upper Triassic Xujiahe (T₃x) Formation (b) generated from 3D seismic interpretation in the DYFS

下载: 全尺寸图片 幻灯片

4.   晚新生代东西向地壳缩短和构造应力场转换

构造应力场转换是现今地球科学研究的热点，晚新生代龙门山褶皱冲断带右旋走滑特征和2008年汶川地震同震斜向挤压逆冲作用暗示，晚新生代存在构造应力场转换。虽然Burchfile等^[20]和Wang等^[52]已就中生代和新生代龙门山褶皱冲断带的最大主应力场方向问题进行过讨论，但构造应力场演化过程及其相应的构造特征仍然不清。自中生代以来，龙门山褶皱冲断带经历了2期构造缩短事件^{[20, 24]}，晚三叠世龙门山在调节SE向逆冲挤压作用下形成了NE东向延伸的龙门山褶皱冲断带和川西前陆盆地，表明早期造山作用为NW-SE向^{[20, 24]}（图 14-a）。晚新生代印度-欧亚大陆碰撞的远距离效应使龙门山褶皱冲断带重新被激活^{[20, 24, 40]}，但对这一阶段青藏高原东缘的挤压应力方向的认识存在分歧。Avouac等^[36]和Zhang等^[37]认为，晚新生代造山作用方向与晚三叠世一致，而另一种观点指出，晚新生代龙门山褶皱冲断带构造变形可能包含重要的右旋走滑分量^{[20, 30, 53]}，后者与汶川地震同震滑移所揭示的斜向逆冲认识一致^{[30-31, 38-39, 44, 46-48, 71]}，暗示晚新生代龙门山褶皱冲断带发生了区域应力场的转变。

图  14  青藏高原东缘龙门山褶皱冲断带及其邻区构造应力场演化

a-晚三叠世水平最大主压应力（σ_hmax）方向为NW-SE向；b-渐新世-早上新世σ_hmax为NW-SE向；c-晚上新世-全新世为EW向挤压过程；d-现今测量的和汶川地震反演的近EW向挤压应力场。WMF-汶川-茂汶断裂；YBF-映秀-北川断裂；PGF-彭县-灌县断裂；QCF-青川断裂；MJF-岷江断裂；HYF-虎牙断裂；LRBF-龙日坝断裂；WLF-武龙断裂；BXF-宝兴断裂；DYF-大邑断层；QXF-邛西断层；PLF-平罗坝断层；LQF-龙泉断层；XPF-熊坡断层

Figure  14.  The evolution of the tectonic stress field in Longmen Mountain and adjacent areas

下载: 全尺寸图片 幻灯片

本次研究表明，大邑背斜区在晚新生代的确存在2期构造事件，分别对应早期的NE向构造和晚期的近NS向构造。平行于龙门山褶皱冲断带的NE向构造存在明显的不连续和被NS向构造切割和限制的特征，而NS向构造在龙门山南段，特别是在四川盆地西南部，呈现为主要的断裂系统，可连续追踪120km^{[43, 52]}。已发表的四川盆地西南部NE向构造的热年代学数据表明，熊坡、平罗和其他NE向构造主要形成于渐新世-中新世^{[52, 55, 72-73]}，而地震相干分析和三维构造建模的结果显示，这一区域内NS向构造切割和限制了NE向构造，由此推断，区域内NS向构造形成于中新世之后。地震剖面和野外测量结果表明，NS走向的断层F5切割了大邑砾岩，暗示NS向构造形成于大邑砾岩底部沉积（约3.6Ma）之后。基于大邑背斜区NS向构造延伸特征及其倾滑的位移量，笔者认为，这一构造吸收了近EW向的构造缩短作用（图 14），而大邑背斜区和邛西背斜区NS向断层和褶皱的最新活动性则反映了龙门山南段局部或区域性水平最大主压应力方向从渐新世-早上新世的NW-SE向（如熊坡、龙泉山背斜区反映的构造应力场方向）^{[20, 52, 55, 72-73]}（图 14-b）转变为晚上新世-全新世的EW向挤压过程^[52]（图 14-c）。在龙门山褶皱冲断带和四川盆地西南部，直接测量地壳应力的方法（例如震源机制解、井眼崩落、水力压裂法、断裂滑动反演等）显示，测得的最大水平主压应力方向在85°~105°之间^{[44, 74]}，与笔者推测的EW向挤压方向结果一致（图 14-d）。位于研究区东部的龙泉山背斜区为一典型的断坡背斜，其走向由NE向转变为近NS向。这一转变可能指示在EW向挤压应力作用下，背斜的重新活动形成近NS向构造^{[30, 51]}。这一最新的近EW向挤压应力形成的NS向构造与青藏高原东缘龙门山褶皱冲断带主体NE向构造的方向存在斜向逆冲分量，而这一滑动分量与汶川地震反映的斜向逆冲同震破裂样式一致^{[30, 38, 71]}（图 14-d）。

龙门山褶皱冲断带南段延伸120km的晚上新世-全新世NS向构造（大邑背斜区和邛西背斜区）切割了渐新世-早上新世的NE向构造和晚更新世名邛冲积扇^{[20, 52, 55, 72-73]}，同时，该NS向构造距汶川地震震中只有20km。在龙门山褶皱冲断带和四川盆地西南部地区重新定位的汶川地震余震揭示，震中以北及其附近存在NS向的地震带，震中可能位于NS向地震带与NE向地震带的交汇处^{[39, 75]}。因此，NS向地震带在未来可能存在一定的地震风险。这一区域发育120km长的NS向构造和汶川地震同震变形均反映出青藏高原东缘最新的EW向地壳缩短过程。

5.   结论

（1）浅部塑性滑脱层自北向南减薄，导致大邑背斜区北部的崇州构造深部以塑性变形为主，而南部的大邑构造深部以脆性断层转折褶皱变形为主。大邑背斜前后翼均发育断层切割能干层，浅部断层活动时间为大邑砾岩沉积（约3.6Ma）之后。

（2）大邑背斜区发育近NS向和NE向2期构造，且近NS向构造切割了NE向构造，反映了龙门山南段局部或区域水平最大主应力转换过程，从渐新世-早上新世的NW-SE向转变为晚上新世-全新世的近EW向。

（3）龙门山南段山前发育的120km长的NS向构造和汶川地震同震变形均反映了青藏高原东缘最新的EW向地壳缩短过程。