基于数值模拟的中生代川东褶皱带构造演化与形成机制探究

基于数值模拟的中生代川东褶皱带构造演化与形成机制探究一、引言1.1研究背景与意义川东褶皱带作为扬子板块西部重要的板内变形带，位于四川盆地与雪峰隆起带之间，其独特的地质构造特征一直是地球科学领域研究的焦点。该褶皱带呈北北东—北东向延伸，总体向北西凸出呈弧形，构造带宽约400km，长约600km。以齐岳山断裂为界，西带为隔档式褶皱，向斜宽、背斜窄，主要由三叠系至侏罗系组成；东带为隔槽式褶皱，背斜宽、向斜窄，卷入地层为南华系至三叠系。这种独特的构造格局是多种地质作用长期相互影响的结果，对其深入研究有助于理解板块运动和地质演化过程。从板块运动角度来看，川东褶皱带处于多个板块的相互作用区域，其构造变形记录了板块之间复杂的碰撞、俯冲和挤压过程。通过研究川东褶皱带，能够为揭示扬子板块与周边板块的相互作用机制提供关键线索，从而进一步完善板块构造理论。例如，有研究认为川东褶皱带的形成与太平洋板块的平板式俯冲有关，这一观点若能得到充分证实，将对板块运动的动力学模型产生重要影响。在地质演化历史中，川东褶皱带经历了多期构造运动，从新元古代的基底形成，到古生代海相沉积盖层的发育，再到中生代的强烈褶皱变形，每一个阶段都在其地质构造上留下了深刻的印记。对这些构造运动的研究，可以帮助我们重建该地区的地质演化历史，了解地壳物质的运动和再分配过程，以及不同地质时期的古地理、古气候环境变化。然而，由于地质构造的复杂性以及地质历史时期的漫长性，仅依靠传统的地质调查和分析方法，难以全面、深入地揭示川东褶皱带的形成机制和演化过程。数值模拟技术的发展为解决这一难题提供了新的途径。数值模拟能够通过建立数学模型，对地质构造过程进行定量分析和模拟，从而弥补传统研究方法的不足。通过数值模拟，可以在计算机上重现川东褶皱带在不同地质时期的构造变形过程，分析各种地质因素对褶皱形成和演化的影响，如岩石力学性质、滑脱层的分布和性质、应力场的变化等。数值模拟技术在川东褶皱带研究中具有重要的作用和意义。一方面，它可以帮助我们验证和完善现有的地质理论和模型。通过将数值模拟结果与实际地质观测数据进行对比，可以检验理论模型的正确性和可靠性，对模型进行修正和完善，从而提高我们对地质过程的认识和理解。另一方面，数值模拟还可以预测地质构造的未来演化趋势，为资源勘探、工程建设和地质灾害防治等提供科学依据。在油气勘探领域，了解川东褶皱带的构造演化历史和现今构造格局，有助于预测油气的分布规律，提高油气勘探的成功率；在工程建设方面，数值模拟可以评估地质构造对工程稳定性的影响，为工程选址和设计提供参考；在地质灾害防治方面，通过模拟地震、滑坡等地质灾害的发生机制和演化过程，可以制定有效的防灾减灾措施，减少灾害损失。1.2研究目的本研究旨在利用数值模拟技术，深入剖析中生代川东褶皱带的构造演化过程和形成机制。通过构建合理的数值模型，模拟不同地质时期的构造应力场、岩石力学性质以及滑脱层的作用，再现川东褶皱带的形成过程，揭示其构造变形的动力学机制。具体而言，本研究的目的包括以下几个方面：揭示构造演化过程：精确确定川东褶皱带在中生代不同阶段的构造变形特征，如褶皱的形态、规模、轴向变化等，梳理其构造演化的详细历程，明确各阶段的关键构造事件和演化顺序。以往的研究虽对川东褶皱带的构造特征有一定认识，但对于其在中生代漫长地质历史时期的动态演化过程，仍缺乏系统且深入的了解。本研究将通过数值模拟，填补这一研究空白，为后续的地质分析提供全面且准确的构造演化信息。分析形成机制：深入探讨控制川东褶皱带形成的主要因素，如板块运动、岩石力学性质、滑脱层的分布和性质等，明确这些因素在褶皱形成过程中的相互作用和影响机制。川东褶皱带的形成是多种复杂因素共同作用的结果，然而，目前对于各因素的具体作用方式和贡献程度，尚未达成一致结论。本研究将运用数值模拟手段，对这些因素进行定量分析，为揭示川东褶皱带的形成机制提供科学依据。验证和完善地质模型：将数值模拟结果与实际地质观测数据进行细致对比，验证和完善现有的川东褶皱带地质模型，提高对其地质构造的认识和理解水平。现有的地质模型在解释川东褶皱带的某些构造现象时，存在一定的局限性。通过数值模拟与实际观测数据的紧密结合，可以对这些模型进行优化和改进，使其更加符合地质实际情况，从而为地质研究提供更可靠的理论支持。为资源勘探和地质灾害防治提供参考：基于对川东褶皱带构造演化和形成机制的深入研究，预测其地下地质构造的分布规律，为油气等资源勘探提供有价值的指导；同时，评估该地区的地质灾害风险，如地震、滑坡等，为地质灾害防治提供科学的决策依据。川东地区蕴含着丰富的油气资源，准确掌握其地质构造分布规律，对于提高油气勘探效率和成功率具有重要意义。此外，该地区地质灾害频发，深入了解地质构造与地质灾害之间的关系，能够为制定有效的防灾减灾措施提供科学指导，保障人民生命财产安全和社会经济的可持续发展。1.3国内外研究现状川东褶皱带作为扬子板块西部重要的板内变形带，其独特的地质构造特征一直是国内外学者关注的焦点。国内外学者在川东褶皱带的地质特征、形成时代和机制等方面展开了大量研究，取得了丰富的成果。同时，数值模拟技术在川东褶皱带研究中的应用也逐渐增多，为深入理解其构造演化提供了新的手段。在地质特征研究方面，早期的研究主要集中在褶皱带的地表形态和构造样式描述。丁道桂等学者于1991年指出，川东褶皱带以齐岳山断裂为界，西带为隔档式褶皱，向斜宽、背斜窄；东带为隔槽式褶皱，背斜宽、向斜窄。随着研究的深入，学者们对褶皱带的地层分布、岩石组合以及断裂构造等进行了详细的调查和分析。研究发现，川东褶皱带内出露的地层从新元古代板溪群到白垩系均有分布，不同地层的岩石组合和变形特征存在明显差异。在断裂构造方面，除了齐岳山断裂外，还识别出了华蓥山断裂等一系列对褶皱带构造演化具有重要影响的断裂。关于川东褶皱带的形成时代，学术界存在一定的争议。Yan等学者在2003年通过构造分析认为，川东构造带内多层逆冲—褶皱构造系形成于中侏罗世至晚白垩世之间；冯向阳等学者同年提出，川东构造带内出现的变形构造起始于印支期—早燕山期，定型于晚燕山期—喜马拉雅期；Wang等学者在2005年根据在雪峰隆起的北东向逆冲断层内获得的5个绢云母40Ar/39Ar同位素年龄数据（217～195Ma），认为雪峰构造带及其西侧川东前陆变形发生在中三叠世到早侏罗世之间；Li和Li在2007年通过对华南岩浆岩的同位素年代学研究后，提出宽约1300km的华南褶皱构造带（包括川东构造带）形成于250～150Ma间，是太平洋区板块平板式俯冲的结果。胡召齐等人于2009年通过对川东地区中生界关键地层之间接触关系的研究，认为川东“侏罗山式”褶皱带形成在晚侏罗世末至早白垩世初期间。在形成机制研究方面，学者们提出了多种观点。一些学者认为，川东褶皱带的形成与板块运动密切相关，是扬子板块与周边板块相互作用的结果。如Li和Li提出的太平洋区板块平板式俯冲导致川东褶皱带形成的观点。另一些学者则强调滑脱层在褶皱带形成中的关键作用。张小琼等人在2015年的研究中指出，中生代川东褶皱带发育着两种不同的褶皱组合型式，以NE向齐岳山断裂为界，东侧为隔槽褶皱，西侧为隔档褶皱，二者在成因上均与不同埋深的滑脱带密切相关。通过数值模拟发现，当滑脱带在隔档褶皱区和隔槽褶皱区分别处于寒武系和角度不整合面向下3-4km的基底深处时，模拟结果能很好地再现该褶皱带的总体形态特征。此外，还有学者认为岩石力学性质、应力场变化等因素也对川东褶皱带的形成起到了重要作用。张必龙等人于2009年采用有限差分法FLAC对川东隔档式褶皱、隔槽式褶皱进行数值模拟，发现互层单元的粘聚力差异与上覆压力是控制隔档式褶皱、隔槽式褶皱形成的主要因素。当互层单元粘聚力差异较小时，出现隔槽式褶皱；当互层单元粘聚力差异较大时，出现隔档式褶皱。并且随着粘聚力差异变大，褶皱形态从箱状向尖棱状过渡。当施加上覆压力后，随着上覆压力的逐渐增大，由粘聚力差异大导致的隔档式褶皱会逐渐向隔槽式褶皱转化，而粘聚力差异小导致的隔槽式褶皱则基本没有变化。随着计算机技术的发展，数值模拟在川东褶皱带研究中得到了广泛应用。数值模拟可以通过建立数学模型，对地质构造过程进行定量分析和模拟，从而弥补传统研究方法的不足。在川东褶皱带的数值模拟研究中，学者们主要采用有限元法、有限差分法等方法，模拟褶皱带的形成过程、应力应变分布以及滑脱层的作用等。张小琼等学者利用FLAC6.0软件模拟了川东褶皱带的两阶段形成过程，即隔槽褶皱区和隔档褶皱区依次形成，结果表明断坡倾角和滑脱带的埋深对褶皱带的形态和形成过程具有重要影响。解国爱等人在2013年通过物理模拟研究了川东侏罗山式褶皱构造带的形成机制，发现滑脱层的存在使得褶皱变形主要发生在滑脱层之上，并且褶皱的形态和规模与滑脱层的性质和厚度密切相关。尽管国内外学者在川东褶皱带的研究中取得了一定的成果，但仍存在一些不足之处。在形成时代的确定上，由于缺乏直接的同位素年代学证据，目前的观点主要基于地质构造分析和地层对比，存在较大的不确定性。在形成机制方面，虽然提出了多种观点，但对于各种因素在褶皱带形成过程中的相对重要性以及它们之间的相互作用机制，还缺乏深入的认识。在数值模拟研究中，目前的模型大多是基于简化的地质条件建立的，对于一些复杂的地质过程，如岩石的非线性变形、多期构造运动的叠加等，还难以进行准确的模拟。此外，数值模拟结果与实际地质观测数据的对比验证还不够充分，需要进一步加强。二、中生代川东褶皱带地质构造特征2.1区域地质背景川东褶皱带主体位于扬子克拉通西部，在扬子板块的构造格局中占据重要位置。其大地构造位置独特，西接四川盆地，东邻雪峰隆起带，北侧为秦岭造山带，向南延伸到滇黔桂交界区。这种特殊的地理位置使其成为多个构造单元相互作用的关键地带，受到了不同方向构造应力的影响，造就了其复杂的地质构造特征。从板块运动的角度来看，扬子板块在地质历史时期经历了复杂的演化过程，川东褶皱带作为其西部的重要变形带，是板块运动的直接产物。在新元古代，扬子板块经历了基底的形成和演化，川东地区靠近雪峰隆起西缘出露的新元古代板溪群，便是扬子板块浅变质基底的一部分。板溪群大部分由紫红色与灰绿色层状杂砂岩—板岩组成，局部可见千枚岩，其中褶皱较紧闭，但变质程度仅达低绿片岩相。这一基底的形成和特征，为后续盖层的沉积和构造变形奠定了基础。在古生代，扬子板块处于相对稳定的海相沉积环境，川东地区也不例外。区内大面积出露的海相盖层由南华系至中三叠统组成，岩性主要为碳酸盐岩和碎屑岩。这些海相沉积地层在长期的地质历史中，记录了当时的沉积环境和构造背景信息。例如，碳酸盐岩的发育表明当时的沉积环境可能为温暖、清澈的浅海，而碎屑岩的存在则反映了陆源物质的输入和沉积环境的变化。中生代是川东褶皱带构造演化的关键时期。在这一时期，川东地区经历了强烈的构造变形，形成了现今所见的褶皱构造格局。上三叠统至上白垩统为陆相碎屑岩，其中上三叠统至上侏罗统属于秦岭造山带南部的前陆沉积，广泛出露在川东构造带的东部。这一时期的沉积特征反映了区域构造背景的变化，从海相沉积转变为陆相沉积，暗示了板块运动和构造应力场的调整。四川盆地作为中国重要的含油气盆地之一，对川东褶皱带的形成和演化产生了深远影响。在地质历史时期，四川盆地经历了多期构造运动，其内部的构造变形和沉积演化与川东褶皱带密切相关。在印支运动时期，四川盆地整体抬升，导致川东地区海相沉积结束，转入陆相沉积阶段。这一构造运动使得川东地区的地层发生了褶皱和断裂，为川东褶皱带的形成奠定了基础。在燕山运动时期，四川盆地受到来自不同方向的构造应力作用，进一步加剧了川东地区的构造变形，使得褶皱构造更加复杂和强烈。雪峰隆起带是扬子板块东南缘的一个重要构造单元，对川东褶皱带的形成同样具有重要作用。雪峰隆起带在地质历史时期经历了多次构造运动和隆升过程，其隆升导致了川东地区的地壳变形和应力调整。在新元古代，雪峰隆起带的基底形成，与川东地区的基底具有一定的亲缘关系。在古生代和中生代，雪峰隆起带的隆升对川东地区的沉积和构造演化产生了重要影响。由于雪峰隆起带的阻挡，使得川东地区的沉积物质来源和沉积环境发生了变化，同时也导致了川东地区的构造应力场发生改变，促进了褶皱构造的形成。一些学者认为，雪峰隆起带的向西推覆作用是川东褶皱带形成的重要动力来源之一。在这种推覆作用下，川东地区的地层发生了多层次拆离推覆及递进挤压，从而形成了隔槽式与隔档式褶皱。2.2地层分布特征川东褶皱带内出露的地层从新元古代到白垩纪均有分布，不同时代的地层具有各自独特的岩性特征和沉积环境，记录了该地区漫长而复杂的地质演化历史。新元古代板溪群是区内出露的最老地层，主要分布在靠近雪峰隆起西缘的区域。其大部分由紫红色与灰绿色层状杂砂岩—板岩组成，局部可见千枚岩。板溪群中的褶皱较为紧闭，这表明其在形成后经历了强烈的构造变形作用。而变质程度仅达低绿片岩相，反映了当时的变质条件相对较低，可能是由于区域构造热事件的强度有限或作用时间较短所致。板溪群作为扬子板块的浅变质基底，为后续地层的沉积和构造演化奠定了基础。在其之上，沉积了南华系至中三叠统的海相盖层，岩性主要为碳酸盐岩和碎屑岩。这些海相沉积地层的形成与当时的海洋环境密切相关，碳酸盐岩的大量发育，暗示了温暖、清澈且水体能量相对稳定的浅海环境，有利于碳酸盐生物的生长和沉积。而碎屑岩的存在，则反映了陆源物质的输入，可能是由于周边陆地的风化剥蚀作用，将碎屑物质搬运至海洋中沉积下来。不同类型岩石的交替出现，也反映了沉积环境的周期性变化，可能与海平面的升降、气候的变迁等因素有关。中生代是川东褶皱带构造演化的关键时期，这一时期的地层特征对于理解褶皱带的形成机制具有重要意义。上三叠统至上白垩统为陆相碎屑岩，标志着沉积环境从海相转变为陆相。上三叠统至上侏罗统属于秦岭造山带南部的前陆沉积，广泛出露在川东构造带的东部。在晚三叠世，印支运动使得整个扬子板块发生了重大的构造变革，川东地区也不例外。区域构造应力场的改变，导致了海相沉积的结束，转入陆相沉积阶段。此时，川东地区处于秦岭造山带南部的拗陷型前陆盆地环境，来自造山带的碎屑物质在盆地中大量堆积，形成了上三叠统至上侏罗统的陆相碎屑岩沉积。这些地层的岩性主要为灰黄、灰绿色薄至厚层粉砂岩夹岩屑砂岩、页岩及菱铁矿结核，底部主要为泥页岩、炭质页岩夹薄煤层。粉砂岩和页岩的大量出现，表明沉积环境相对安静，水体能量较低；而岩屑砂岩的存在，则反映了陆源碎屑物质的快速堆积。菱铁矿结核和薄煤层的发育，暗示了当时的沉积环境可能为潮湿的沼泽或湖泊，有利于植物的生长和保存。白垩系红层局部出现在川东构造带东部的小型盆地内，如沅麻（沅陵—麻阳）盆地和江汉盆地西北部的当阳盆地。白垩系又分为上、下白垩统，其地层时代归属已经得到了古生物证据和地层对比的支持。在沅麻盆地，根据芷江燕子岩剖面发现的双壳类化石，将白垩系分为下统和上统；辰溪县柘溪西南恐龙足印化石的发现，限定了沅麻盆地小洞组近顶部地层应属上白垩统。在江汉盆地西部，五龙组中所含孢粉经研究后被定位于早白垩世晚期；红花套组中采获的介形类和跑马岗组介形虫类化石、轮藻化石，表明这两套地层应属晚白垩世地层，白垩系改为三分后这两套地层应属于上白垩统。白垩系红层的形成，可能与当时的气候干旱、氧化环境增强有关。在这种环境下，陆源碎屑物质在盆地中堆积，并受到氧化作用的影响，形成了红色的砂岩和泥岩。白垩系红层的分布和特征，也反映了川东地区在白垩纪时期的构造活动和沉积环境的变化，小型盆地的形成可能与区域构造应力场的局部调整有关，导致地壳发生差异性升降运动，形成了这些小型的沉积盆地。2.3褶皱构造样式2.3.1隔档式褶皱隔档式褶皱是川东褶皱带西带的典型构造样式，其形态独特，具有鲜明的特征。在地表出露的隔档式褶皱主要由三叠系至侏罗系组成，其显著特点是向斜宽、背斜窄，整体呈现出一种独特的“窄背宽向”的几何形态。向斜和背斜多为具有平底或平顶的箱状褶皱，这种箱状形态的形成与地层的力学性质以及构造应力的作用方式密切相关。在构造变形过程中，由于地层的岩石力学性质差异，以及受到的水平挤压应力和垂直上覆压力的共同作用，使得褶皱在发育过程中形成了较为平缓的底部和顶部，从而呈现出箱状的形态。隔档式褶皱的走向具有一定的规律性。自南向北，其走向由南北至北北东逐渐过渡到北东向，总体呈向北西突出的弧形。这种弧形的走向是区域构造应力场和基底构造共同作用的结果。在区域构造应力场中，受到来自不同方向应力的影响，使得褶皱在形成过程中发生了一定程度的弯曲和旋转，从而呈现出弧形的形态。而基底构造的不均匀性，也对褶皱的走向产生了重要影响，使得褶皱在不同地段的走向发生了变化。褶皱轴面的倾向在不同地段也有所不同。一般来说，轴面较陡立，但多以向西倾为主，仅在华蓥山断裂带附近轴面向东倾。轴面倾向的这种变化与断裂带的活动密切相关。华蓥山断裂带作为一条重要的构造边界，其活动对周边地层的变形产生了显著影响。在断裂带附近，由于受到断裂活动的牵引和挤压作用，使得褶皱轴面的倾向发生了改变，从而呈现出与其他地段不同的特征。以华蓥山背斜为例，它是川东隔档式褶皱中的典型代表。华蓥山背斜核部出露了古生界地层，这表明其在构造演化过程中经历了更为复杂的变形历史。与其他背斜相比，华蓥山背斜的核部地层更为古老，反映了其在早期构造运动中受到的影响更为强烈。其轴面在华蓥山断裂带附近向东倾，这一特征与华蓥山断裂带的活动密切相关。在断裂带活动过程中，对背斜的轴面产生了牵引和挤压作用，使得轴面发生了向东的倾斜。这种轴面倾向的变化，不仅影响了背斜的形态和构造特征，也对其内部的地层结构和岩石变形产生了重要影响。华蓥山背斜的两翼地层倾角相对较大，这是由于在褶皱形成过程中，受到强烈的挤压作用，使得地层发生了较大程度的弯曲和变形。较大的地层倾角也反映了该背斜在构造演化过程中所承受的构造应力强度较大，经历了较为复杂的构造变形过程。2.3.2隔槽式褶皱隔槽式褶皱是川东褶皱带东带的典型构造样式，与隔档式褶皱形成鲜明对比，具有独特的地质特征。其卷入地层为南华系至三叠系，呈现出背斜宽、向斜窄的显著特征，与隔档式褶皱的“窄背宽向”形成了鲜明的反差。这种背斜与向斜宽度的差异，反映了不同的构造变形机制和地层响应。在隔槽式褶皱的形成过程中，由于受到的构造应力和地层力学性质的影响，使得背斜部位的地层相对较为舒展，而向斜部位的地层则受到强烈的挤压和收缩，从而导致了背斜宽、向斜窄的形态特征。各个背斜与向斜同样以箱状褶皱为主，与隔档式褶皱中的箱状褶皱类似，其形成与地层的力学性质以及构造应力的作用密切相关。在构造变形过程中，地层受到水平挤压应力和垂直上覆压力的共同作用，使得褶皱在发育过程中形成了较为平缓的底部和顶部，从而呈现出箱状的形态。不同的是，由于隔槽式褶皱的地层组成和构造应力环境与隔档式褶皱有所不同，其箱状褶皱的具体形态和规模也存在一定的差异。背斜核部多出露南华系—奥陶系，而向斜核部主要为三叠系。这种地层出露特征反映了隔槽式褶皱在构造演化过程中的地层变形和隆升情况。在构造运动过程中，背斜部位由于受到的挤压作用相对较小，地层相对较为稳定，使得较老的南华系—奥陶系得以保存并出露于核部；而向斜部位则受到强烈的挤压和沉降作用，使得较新的三叠系地层沉积并保存于核部。隔槽式褶皱带的走向自南向北由北北东逐渐变化成北东东向，总体呈向西北凸出的弧形，与隔档式褶皱带的弧形走向相似，但具体的走向变化和形态特征存在一定差异。这种弧形走向同样是区域构造应力场和基底构造共同作用的结果。在区域构造应力场的作用下，地层发生了弯曲和变形，形成了弧形的褶皱带。而基底构造的不均匀性，也对褶皱带的走向和形态产生了重要影响，使得褶皱带在不同地段的走向和形态发生了变化。轴面同样陡立，但多向南东倾，这与隔档式褶皱轴面多向西倾的特征不同。轴面倾向的差异反映了隔槽式褶皱和隔档式褶皱在构造应力作用方式和变形机制上的差异。在隔槽式褶皱的形成过程中，受到的构造应力方向和作用方式使得轴面发生了向南东的倾斜。这种轴面倾向的差异，不仅影响了褶皱的形态和构造特征，也对其内部的地层结构和岩石变形产生了重要影响。以桑植-石门复向斜为例，它是隔槽式褶皱的一个典型实例。桑植-石门复向斜内褶皱具“复杂褶皱”组合样式，主要形成于中三叠世后期印支运动。在印支运动期间，区域构造应力场发生了重大变化，使得该地区的地层发生了强烈的褶皱变形。复向斜内的褶皱基底变形主要以逆冲断裂为主，这表明在褶皱形成过程中，基底地层受到了强烈的挤压和剪切作用，导致了逆冲断裂的发育。这种复杂的褶皱组合样式和基底变形特征，反映了桑植-石门复向斜在构造演化过程中的复杂性和多阶段性。其形成过程不仅受到了印支运动的影响，还可能受到了后期构造运动的叠加和改造，使得褶皱的形态和构造特征更加复杂多样。2.4断裂构造特征齐岳山断裂是川东褶皱带的重要构造边界，对褶皱带的构造演化和变形特征具有重要影响。在平面上，齐岳山断裂呈S形分布，这种独特的形态是由其在不同构造阶段受到的应力作用和地层变形所导致的。它由多条不连续断层组成，且在北部呈右阶左行雁列式，在南部左阶右行排列。这种排列方式反映了断裂在形成和演化过程中受到了不同方向应力的作用，导致断层在平面上呈现出复杂的组合形态。齐岳山断裂的走向为北北东向，与川东褶皱带的总体走向基本一致。其倾角较陡，这表明在断裂形成过程中，受到了强烈的挤压和剪切作用，使得地层发生了急剧的错动和变形。齐岳山断裂在晚侏罗世主滑脱层抬升导致隔槽式褶皱带向隔档式转变过程中形成雏形，这一时期的构造运动使得区域应力场发生了重大变化，导致主滑脱层抬升，从而引发了褶皱带的变形和断裂的形成。随后，齐岳山断裂经历了晚白垩世黔中地块的改造作用和新近纪末期的差异隆升，最终形成现今的形态。在晚白垩世，黔中地块的运动对齐岳山断裂产生了重要影响，导致其进一步变形和演化；而新近纪末期的差异隆升，则使得断裂两侧的地层发生了不同程度的升降运动，进一步塑造了断裂的现今形态。齐岳山断裂控制形成了西侧背斜和西侧斜坡两类转换带，并分别发育“基底冲断一断层相关褶皱一滑脱”和“基底逆冲一滑脱”两类构造变形模式。在西侧背斜转换带，基底冲断作用导致地层发生强烈的变形，形成了断层相关褶皱，同时伴随着滑脱作用，使得地层在不同层次之间发生了相对滑动。这种变形模式使得背斜构造更加复杂，地层的连续性受到破坏，形成了一系列的断层和褶皱组合。在西侧斜坡转换带，基底逆冲作用使得地层向上逆冲，同时伴随着滑脱作用，使得地层在斜坡上发生了滑动和变形。这种变形模式导致斜坡地区的地层结构不稳定，容易引发地质灾害。华蓥山断裂同样是川东褶皱带内的一条重要断裂，走向北北东，与齐岳山断裂的走向具有一定的相似性。这种相似性反映了它们在区域构造应力场中的共同作用和响应。华蓥山断裂倾向东，倾角较陡，其力学性质表现为压扭性。这表明在断裂形成过程中，受到了强烈的挤压和扭动作用，使得地层发生了复杂的变形。华蓥山断裂对川东褶皱带西带的隔档式褶皱具有重要控制作用。在褶皱形成过程中，华蓥山断裂的活动导致了地层的变形和错动，使得隔档式褶皱的形态和分布受到了影响。华蓥山背斜核部出露古生界地层，且轴面在华蓥山断裂带附近向东倾，这与华蓥山断裂的活动密切相关。在断裂活动过程中，对背斜的轴面产生了牵引和挤压作用，使得轴面发生了向东的倾斜，同时也导致了背斜核部地层的出露和变形。三、数值模拟原理与方法3.1数值模拟在地质学中的应用概述数值模拟作为一种重要的研究手段，在地质学领域的应用日益广泛且深入，为地质研究带来了全新的视角和方法。随着计算机技术的飞速发展，数值模拟技术能够处理日益复杂的地质问题，成为地质学研究不可或缺的工具。在地质构造研究方面，数值模拟发挥着关键作用。地质构造的形成和演化是一个漫长而复杂的过程，受到多种因素的综合影响。通过数值模拟，可以在计算机上构建地质模型，模拟不同地质时期的构造应力场、岩石力学性质以及各种地质作用的相互关系，从而再现地质构造的形成过程。在研究板块运动对地质构造的影响时，数值模拟能够模拟板块的碰撞、俯冲和分离等过程，分析这些过程中地壳的变形、应力分布以及岩浆活动等现象，为理解板块构造理论提供了重要的定量依据。例如，在研究造山带的形成机制时，数值模拟可以模拟板块碰撞过程中地壳的缩短、加厚和隆升，揭示造山带的构造演化历史。通过模拟不同的边界条件和岩石力学参数，可以探讨造山带的形态、结构以及内部构造的形成原因，为造山带的研究提供了重要的参考。在矿产资源勘探领域，数值模拟同样具有重要意义。矿产资源的形成与地质构造、岩石性质以及流体运移等因素密切相关。通过数值模拟，可以建立矿产资源形成的地质模型，模拟成矿过程中物质的迁移、富集和沉淀，预测矿产资源的分布规律。这有助于指导矿产资源的勘探工作，提高勘探效率和成功率。在石油勘探中，数值模拟可以模拟油气的生成、运移和聚集过程，分析储层的物性特征和油气分布规律，为油气勘探提供重要的依据。通过建立三维地质模型，结合地球物理和地球化学数据，数值模拟可以预测潜在的油气藏位置和规模，为勘探决策提供科学支持。地质灾害预测也是数值模拟的重要应用领域之一。地震、滑坡、泥石流等地质灾害严重威胁着人类的生命财产安全和社会的可持续发展。数值模拟可以模拟地质灾害的发生机制和演化过程，预测灾害的发生概率和影响范围，为地质灾害的防治提供科学依据。在地震研究中，数值模拟可以模拟地震波的传播、地震的破裂过程以及地震对建筑物的影响，评估地震灾害的风险。通过建立地震地质模型，结合地震监测数据，数值模拟可以预测地震的发生时间、地点和震级，为地震预警和灾害防范提供重要的支持。对于川东褶皱带的研究而言，数值模拟技术具有独特的优势和重要的意义。川东褶皱带的地质构造复杂，受到多种地质因素的影响，传统的研究方法难以全面、深入地揭示其形成机制和演化过程。数值模拟技术能够综合考虑各种地质因素，建立精确的地质模型，模拟川东褶皱带在不同地质时期的构造变形过程，分析褶皱、断裂的形成机制以及它们之间的相互关系。通过数值模拟，可以探讨板块运动、岩石力学性质、滑脱层的分布和性质等因素对川东褶皱带形成的影响，为解释其独特的构造特征提供定量依据。与传统地质研究方法相比，数值模拟具有以下优势：一是能够处理复杂的地质过程和多因素相互作用，而传统方法往往难以全面考虑这些因素；二是可以进行定量分析，提供精确的数值结果，而传统方法多以定性描述为主；三是可以在不同的假设条件下进行模拟实验，探索各种可能性，而传统方法受限于实际地质条件，难以进行这样的探索。3.2常用数值模拟方法介绍3.2.1有限差分法有限差分法（FiniteDifferenceMethod，FDM）是一种将连续问题离散化的数值方法，其核心原理是用差分来近似替代微分。在地质学研究中，地质体的物理过程通常可以用偏微分方程来描述，但这些方程往往难以直接求解。有限差分法通过将连续的求解区域划分为有限个网格点，把偏微分方程在这些网格点上离散化，从而将连续的无限自由度问题转化为离散的有限自由度问题，最终通过求解代数方程组得到近似解。以一维热传导问题为例，热传导方程为：\frac{\partialT}{\partialt}=\alpha\frac{\partial^{2}T}{\partialx^{2}}，其中T是温度，t是时间，x是空间坐标，\alpha是热扩散系数。为了使用有限差分法求解，首先将时间和空间进行离散化。将时间t划分为t_{n}（n=0,1,2,\cdots），空间x划分为x_{i}（i=0,1,2,\cdots）。然后，用差分近似导数，例如，对于\frac{\partialT}{\partialt}，可以用向前差分近似为\frac{T_{i}^{n+1}-T_{i}^{n}}{\Deltat}，对于\frac{\partial^{2}T}{\partialx^{2}}，可以用中心差分近似为\frac{T_{i+1}^{n}-2T_{i}^{n}+T_{i-1}^{n}}{\Deltax^{2}}，其中\Deltat和\Deltax分别是时间步长和空间步长。将这些差分近似代入热传导方程，就得到了离散的差分方程：\frac{T_{i}^{n+1}-T_{i}^{n}}{\Deltat}=\alpha\frac{T_{i+1}^{n}-2T_{i}^{n}+T_{i-1}^{n}}{\Deltax^{2}}。通过求解这个差分方程，就可以得到在不同时间和空间点上的温度近似值。在模拟地质构造变形时，FLAC（FastLagrangianAnalysisofContinua）软件是一款常用的基于有限差分法的软件。它能够模拟地质材料在复杂受力条件下的力学行为，考虑岩石的非线性、大变形等特性。在模拟川东褶皱带的构造变形时，研究人员可以利用FLAC软件建立地质模型，将川东地区的地层结构、岩石力学参数等信息输入模型中。通过设置边界条件和加载方式，模拟不同地质时期的构造应力作用下地层的变形过程。FLAC软件会根据有限差分法的原理，将模型划分为有限个网格单元，在每个单元上求解力学平衡方程，从而得到整个模型的应力、应变和位移分布。通过对模拟结果的分析，可以深入了解川东褶皱带的形成机制和演化过程，如褶皱的形成、发展和变形特征，以及断裂的产生和演化等。有限差分法的优点在于算法简单直观，易于理解和实现，对于一些简单的地质模型能够快速得到结果。然而，它也存在一些局限性。例如，有限差分法对网格的依赖性较强，网格的划分方式和大小会直接影响计算结果的精度和稳定性。在处理复杂的地质边界条件和不规则的地质体形状时，有限差分法的网格划分往往比较困难，可能会导致计算误差增大。此外，有限差分法在处理大规模问题时，计算量较大，计算效率相对较低。3.2.2有限元法有限元法（FiniteElementMethod，FEM）是一种广泛应用于工程和科学计算的数值方法，在地质学领域也发挥着重要作用。其基本原理是将连续的求解域离散为有限个单元的组合体，把连续的无限自由度问题离散成以未知场函数的结点值为未知量的有限自由度问题，通过求解高阶代数方程组来求出这些未知量。有限元法的实现过程主要包括以下几个关键步骤。首先是结构的离散化，将分析结构系统划分成有限个单元体，并在单元体的指定点设置节点，把相邻的单元体在节点处连接起来组成单元的集合体，以代替原来的结构。一般来说，单元划分越细，表述变形情况就越精确，越接近实际变形，但同时计算量也会越大。例如，在模拟川东褶皱带的地质构造时，可以将该区域的地层划分成三角形、四边形或四面体等不同形状的单元，这些单元的大小和分布需要根据研究区域的地质特征和计算精度要求进行合理设计。选择合适的位移模式也是有限元法的重要环节。通过设定位移模式，可以导出用节点位移表示单元内任意点位移的表达式。常用的位移模式有线性位移模式、二次位移模式等，不同的位移模式适用于不同的问题和单元类型。在模拟川东褶皱带的变形时，需要根据岩石的力学性质和变形特点选择合适的位移模式，以准确描述地层的变形行为。建立单元结点力和结点位移之间的关系是有限元法的核心步骤之一。根据材料力学和弹性力学的基本原理，通过对单元进行力学分析，建立起单元结点力与结点位移之间的刚度矩阵关系。这个刚度矩阵反映了单元的力学特性，是求解整个结构平衡方程的基础。在川东褶皱带的模拟中，需要考虑岩石的非线性力学行为，如塑性变形、蠕变等，因此刚度矩阵的建立需要采用相应的非线性本构模型。计算结点荷载，包括集中力、面力、体积力等，将这些荷载等效到节点上，以便后续建立平衡方程。集合所有单元的刚度方程，建立整个结构的平衡方程。这个平衡方程是一个高阶代数方程组，通过求解该方程组，可以得到节点的位移解。最后，由位移解进一步计算出应变和应力，从而得到整个结构的力学响应。在模拟复杂地质问题时，有限元法具有显著的优势。它能够灵活地处理复杂的几何形状和边界条件，对于川东褶皱带这种地质构造复杂、边界条件多样的区域，有限元法可以通过合理的单元划分和边界条件设置，准确地模拟地质体的形态和边界情况。有限元法还可以考虑多种物理场的耦合作用，如应力场、温度场、渗流场等。在川东褶皱带的形成过程中，岩石的变形不仅受到构造应力的作用，还可能受到温度变化和地下水渗流等因素的影响，有限元法能够综合考虑这些因素，更全面地揭示地质过程的本质。有限元法在模拟地质构造演化、矿产资源勘探、地质灾害评估等方面都有广泛的应用。在研究川东褶皱带的构造演化时，有限元法可以模拟不同地质时期的构造应力场变化，分析褶皱和断裂的形成机制，预测地质构造的未来发展趋势。在矿产资源勘探中，有限元法可以通过模拟地质体的物理性质和力学行为，预测矿产资源的分布规律，为勘探工作提供指导。在地质灾害评估方面，有限元法可以模拟地震、滑坡等灾害的发生过程，评估灾害对地质体和建筑物的影响，为灾害防治提供科学依据。3.3本次研究采用的数值模拟方法及软件在本次对中生代川东褶皱带的研究中，我们选用了基于有限差分法的FLAC（FastLagrangianAnalysisofContinua）软件作为主要的数值模拟工具。FLAC软件在处理复杂地质力学问题方面具有独特的优势，这使其成为研究川东褶皱带构造演化的理想选择。从算法原理上看，FLAC采用显式有限差分格式来求解运动方程和连续性方程。这种算法能够高效地处理大变形和非线性问题，非常适合模拟地质构造演化过程中岩石的复杂力学行为。在川东褶皱带的形成过程中，岩石经历了强烈的挤压和变形，呈现出非线性的力学响应。FLAC软件能够准确地模拟这种非线性行为，通过迭代计算逐步逼近真实的变形状态，从而为我们提供详细的应力、应变和位移信息。与其他数值模拟方法相比，有限差分法在处理大变形问题时具有较高的计算效率和稳定性，能够快速地得到可靠的模拟结果。FLAC软件提供了丰富的本构模型，这些模型可以准确地描述不同岩石类型的力学特性。在川东褶皱带中，岩石类型多样，包括砂岩、页岩、灰岩等，它们具有不同的力学性质。FLAC软件中的弹性模型、弹塑性模型、粘弹性模型等，可以根据岩石的实际情况进行选择和调整，从而精确地模拟不同岩石在构造应力作用下的变形和破坏过程。对于脆性的砂岩和灰岩，我们可以选用合适的弹塑性模型来描述其在受力过程中的屈服和破坏行为；对于具有一定粘性的页岩，则可以采用粘弹性模型来考虑其随时间变化的力学响应。通过合理选择本构模型，FLAC软件能够真实地反映川东褶皱带中岩石的力学行为，为研究构造演化提供准确的基础数据。在模拟复杂地质构造方面，FLAC软件也表现出了卓越的能力。它能够处理复杂的几何形状和边界条件，这对于模拟川东褶皱带这种具有复杂地质构造的区域至关重要。川东褶皱带内存在着众多的褶皱、断裂和滑脱层，这些地质构造的几何形状和分布特征对构造演化过程有着重要的影响。FLAC软件可以通过灵活的网格划分技术，准确地模拟这些复杂的地质构造，并且能够考虑它们之间的相互作用。在模拟褶皱和断裂的相互作用时，FLAC软件可以根据地质模型的特点，合理地设置网格和边界条件，准确地模拟断裂对褶皱形态和发展的影响，以及褶皱对断裂活动的控制作用。这种对复杂地质构造的精确模拟能力，使得FLAC软件能够为我们揭示川东褶皱带构造演化的内在机制提供有力的支持。FLAC软件还具有良好的后处理功能，能够直观地展示模拟结果。通过可视化的方式，我们可以清晰地观察到褶皱带在不同演化阶段的应力、应变和位移分布情况，以及不同地质构造要素之间的相互关系。这有助于我们深入理解川东褶皱带的构造演化过程，为进一步的分析和研究提供了便利。我们可以通过后处理功能生成应力云图、应变等值线图和位移矢量图等，直观地展示褶皱带内应力、应变和位移的分布特征，从而快速地识别出构造变形的关键区域和主要影响因素。四、川东褶皱带数值模拟模型构建4.1模型建立的依据和假设构建川东褶皱带数值模拟模型的首要依据是丰富且详实的地质资料。通过对川东地区的地质调查，我们获取了关于地层分布、褶皱构造、断裂构造等方面的第一手资料。这些资料为模型的建立提供了坚实的基础，使得模型能够尽可能真实地反映川东褶皱带的地质特征。在研究地层分布时，我们详细了解了新元古代板溪群、南华系至中三叠统的海相盖层以及上三叠统至上白垩统的陆相碎屑岩的分布范围、厚度变化等信息，这些信息对于准确设定模型中的地层参数至关重要。前人的研究成果也是模型建立的重要参考。众多学者在川东褶皱带的地质构造、形成机制等方面进行了深入研究，他们的研究成果为我们提供了宝贵的思路和方法。一些学者对川东褶皱带的褶皱构造样式、断裂构造特征进行了详细的分析，这些研究成果帮助我们更好地理解了川东褶皱带的地质构造规律，从而在模型建立过程中能够更准确地模拟这些构造特征。在岩石力学参数方面，由于川东褶皱带内岩石类型多样，不同岩石的力学性质存在显著差异。为了准确描述这些岩石的力学行为，我们参考了相关的岩石力学实验数据和经验公式。对于砂岩、页岩、灰岩等主要岩石类型，我们根据已有的研究成果，确定了它们的弹性模量、泊松比、粘聚力、内摩擦角等力学参数。在缺乏实测数据的情况下，我们采用了合理的经验值，并结合地质背景进行了适当的调整。边界条件的假设对于数值模拟结果的准确性同样至关重要。在川东褶皱带的数值模拟中，我们根据区域构造背景和地质演化历史，对模型的边界条件进行了合理假设。考虑到川东褶皱带位于扬子板块西部，受到了来自不同方向的构造应力作用。我们假设模型的边界受到水平挤压应力的作用，应力方向与川东褶皱带的总体走向相关。在模型的底部，我们假设为固定边界，以模拟岩石的深部约束条件；在模型的侧面，根据实际地质情况，设置了相应的位移边界条件或应力边界条件。为了简化模型，我们还做出了一些必要的假设。假设岩石为连续介质，忽略岩石中的微小孔隙和裂隙对整体力学行为的影响。虽然在实际地质情况中，岩石中存在着大量的孔隙和裂隙，它们对岩石的力学性质和变形行为有着重要的影响。但在本次模拟中，为了降低模型的复杂性，我们在一定程度上忽略了这些因素。同时，我们假设模型在模拟过程中不受温度、流体等因素的影响。尽管在实际地质过程中，温度和流体的作用不可忽视，但在本次研究中，为了集中研究构造应力对川东褶皱带形成的影响，我们暂时忽略了这些因素，以便更清晰地分析构造应力与褶皱形成之间的关系。4.2模型参数设置4.2.1岩石力学参数岩石力学参数的准确设定对于数值模拟结果的可靠性至关重要。在川东褶皱带的数值模拟中，不同地层岩石具有各自独特的力学性质，这些性质直接影响着褶皱和断裂的形成与演化过程。对于新元古代板溪群，其主要岩石类型为紫红色与灰绿色层状杂砂岩—板岩，局部为千枚岩。根据相关岩石力学实验研究以及类似地质背景下岩石的力学参数参考，设定其弹性模量为E_1=8\times10^4MPa，泊松比\mu_1=0.25，粘聚力C_1=20MPa，内摩擦角\varphi_1=35^{\circ}。板溪群作为扬子板块的浅变质基底，其较高的弹性模量反映了其相对较强的抵抗变形能力，而粘聚力和内摩擦角则决定了其在受力过程中的剪切强度和稳定性。南华系至中三叠统的海相盖层，岩性主要为碳酸盐岩和碎屑岩。其中，碳酸盐岩（如灰岩）的弹性模量E_2=6\times10^4MPa，泊松比\mu_2=0.23，粘聚力C_2=15MPa，内摩擦角\varphi_2=30^{\circ}；碎屑岩（如砂岩）的弹性模量E_3=4\times10^4MPa，泊松比\mu_3=0.22，粘聚力C_3=10MPa，内摩擦角\varphi_3=28^{\circ}。碳酸盐岩和碎屑岩力学参数的差异，导致它们在构造应力作用下的变形行为有所不同。碳酸盐岩相对较高的弹性模量使其在受力时更倾向于发生脆性变形，容易产生断裂；而碎屑岩较低的弹性模量则使其更容易发生塑性变形，形成褶皱。上三叠统至上白垩统的陆相碎屑岩，主要由粉砂岩、页岩等组成。粉砂岩的弹性模量E_4=3\times10^4MPa，泊松比\mu_4=0.24，粘聚力C_4=8MPa，内摩擦角\varphi_4=25^{\circ}；页岩的弹性模量E_5=2\times10^4MPa，泊松比\mu_5=0.26，粘聚力C_5=6MPa，内摩擦角\varphi_5=22^{\circ}。陆相碎屑岩的力学性质相对较弱，尤其是页岩，其较低的弹性模量、粘聚力和内摩擦角，使其在构造应力作用下容易发生塑性流动和变形，常常作为滑脱层影响上覆地层的变形。这些岩石力学参数并非固定不变，在实际模拟过程中，需要根据模拟结果与实际地质观测数据的对比进行适当调整。如果模拟得到的褶皱形态与实际地质情况差异较大，可能需要对岩石的弹性模量和粘聚力等参数进行微调，以更好地拟合实际地质构造。在模拟隔档式褶皱时，如果模拟结果中背斜和向斜的宽度比例与实际情况不符，可以尝试调整相关地层岩石的力学参数，观察模拟结果的变化，直到得到与实际地质构造较为吻合的结果。4.2.2边界条件设定边界条件的合理设定是数值模拟的关键环节之一，它直接影响着模型内部的应力分布和变形状态，进而决定了模拟结果的准确性和可靠性。在川东褶皱带的数值模拟中，边界条件的设定需要综合考虑区域构造背景、地质演化历史以及岩石力学性质等多方面因素。在模型的底部，由于受到深部岩石的约束，假设为固定边界，即底部节点在三个方向上的位移均为零。这种假设能够模拟深部岩石对上部地层的支撑作用，限制了模型底部的变形，使其符合实际地质情况中深部岩石相对稳定的特点。在实际地质构造中，深部岩石的变形相对较小，其对上部地层的约束作用显著，固定边界条件能够较好地反映这一特征。模型的侧面边界条件根据实际地质情况进行了相应设置。考虑到川东褶皱带位于扬子板块西部，受到了来自不同方向的构造应力作用，假设模型的左右两侧边界受到水平挤压应力的作用。应力方向与川东褶皱带的总体走向相关，大致为北北东—北东向。这种水平挤压应力的施加，模拟了板块运动对川东地区的构造影响，是导致川东褶皱带形成的主要动力来源之一。在模拟过程中，通过调整水平挤压应力的大小和方向，可以研究不同构造应力条件下川东褶皱带的形成和演化过程。当增大水平挤压应力的大小时，观察褶皱的发育程度和形态变化，分析应力大小对褶皱形成的影响。在顶部边界，考虑到上覆地层的压力作用，施加了相应的垂向压力。垂向压力的大小根据地层的厚度和岩石密度进行计算，以模拟上覆地层对下部地层的加载作用。上覆地层的压力会影响岩石的力学行为和变形方式，在模拟中准确考虑这一因素，能够使模拟结果更加符合实际地质情况。随着垂向压力的增加，岩石的抗压强度增大，褶皱的形成和发展可能会受到一定的抑制，通过模拟可以观察到这种变化趋势。在模型的加载方式上，采用逐步加载的方式模拟地质构造演化过程。这种加载方式能够模拟地质历史时期构造应力的逐渐积累和释放过程，使模拟结果更具真实性。在模拟中生代川东褶皱带的形成过程时，根据地质演化历史，将整个过程划分为多个阶段，每个阶段施加一定的构造应力，观察地层在不同阶段的变形响应。通过逐步加载，可以清晰地看到褶皱和断裂的形成和发展过程，以及它们在不同构造阶段的变化特征。4.3网格划分与模型初始化在构建川东褶皱带数值模拟模型时，网格划分是一项关键步骤，其质量直接影响到模拟结果的精度和计算效率。本研究采用了结构化网格划分方法，根据川东褶皱带的地质构造特征和模型的几何形状，将模型区域划分为规则的四边形网格。这种网格划分方式具有较高的计算效率和稳定性，能够准确地描述地质体的几何形状和边界条件。在网格划分过程中，遵循了以下原则：在地质构造复杂的区域，如褶皱和断裂发育的部位，采用较小的网格尺寸，以提高模拟结果的精度。这是因为在这些区域，岩石的变形和应力分布较为复杂，需要更精细的网格来准确描述其变化。在齐岳山断裂和华蓥山断裂附近，以及褶皱的核部和翼部，减小网格尺寸，使网格能够更好地捕捉到这些区域的构造特征和应力应变变化。在地质构造相对简单的区域，适当增大网格尺寸，以减少计算量，提高计算效率。这样可以在保证模拟精度的前提下，合理控制计算资源的使用，提高模拟的可行性和效率。为了确定合适的网格尺寸，进行了网格敏感性分析。通过设置不同的网格尺寸，对同一模型进行多次模拟，并对比模拟结果的差异。结果表明，当网格尺寸过小时，虽然计算精度会提高，但计算量会大幅增加，计算时间显著延长；当网格尺寸过大时，模拟结果的精度会受到影响，无法准确反映地质构造的细节特征。经过反复测试和分析，最终确定在褶皱和断裂附近采用500m×500m的网格尺寸，在其他区域采用1000m×1000m的网格尺寸。这样的网格划分方案在保证计算精度的同时，能够有效控制计算量，提高计算效率。完成网格划分后，对模型进行了初始化。初始化过程主要包括设置初始应力场和位移场。根据区域地质背景和前人研究成果，假设模型初始处于重力平衡状态，通过计算岩石的自重应力，确定初始应力场。在计算自重应力时，考虑了不同地层岩石的密度差异，根据岩石力学参数中设定的密度值，计算每个网格单元的自重应力。初始位移场设置为零，即模型在初始时刻没有发生位移。在模型初始化过程中，还对模型的边界条件进行了再次检查和确认，确保边界条件的设置符合实际地质情况和模拟要求。对模型底部的固定边界条件、侧面的水平挤压应力边界条件以及顶部的垂向压力边界条件进行了仔细检查，确保边界条件的施加准确无误。同时，对模型的材料参数进行了检查，确保岩石力学参数等材料信息的输入正确，以保证模拟结果的可靠性。五、模拟结果与分析5.1模拟结果展示5.1.1应力应变分布特征在数值模拟过程中，获取了不同演化阶段的应力应变分布云图，这些云图为深入理解川东褶皱带的构造变形机制提供了关键信息。在模拟的初始阶段，模型受到水平挤压应力的作用，应力主要集中在模型的边界部位，尤其是受到挤压的两侧边界。随着模拟的进行，应力逐渐向模型内部传递，在褶皱开始发育的区域，应力集中现象愈发明显。在齐岳山断裂和华蓥山断裂附近，应力值显著增大，这是由于断裂带的存在使得岩石的连续性受到破坏，应力在这些部位发生了集中和重新分布。在褶皱的核部和翼部，也出现了明显的应力集中现象。在背斜核部，岩石受到强烈的挤压和弯曲，导致应力集中；而在向斜翼部，由于地层的倾斜和变形，应力也相对集中。通过对不同阶段应力分布云图的对比分析，可以清晰地观察到应力集中区域的迁移和变化，以及应力值的逐渐增大过程，这反映了褶皱在形成和发展过程中，岩石所承受的应力不断积累和调整。应变分布云图同样展示出独特的特征。在模拟初期，应变主要集中在模型的薄弱部位，这些部位通常是岩石力学性质相对较弱的区域，如页岩层或存在节理、裂隙的区域。随着水平挤压应力的持续作用，应变逐渐向整个模型扩展，并且在褶皱发育的区域，应变值明显增大。在隔档式褶皱和隔槽式褶皱的形成区域，应变分布呈现出不同的特征。在隔档式褶皱区域，向斜部位的应变相对较大，这是因为向斜部位的地层受到的挤压和弯曲作用更为强烈，导致岩石发生了较大程度的变形。而在隔槽式褶皱区域，背斜部位的应变相对较大，这与隔槽式褶皱的构造特征有关，背斜宽、向斜窄，背斜部位的地层在形成过程中经历了更复杂的变形过程，从而导致应变集中。随着模拟的进一步进行，应变分布逐渐趋于稳定，但在褶皱和断裂的关键部位，仍然保持着较高的应变值，这表明这些部位的岩石变形较为强烈，是构造变形的主要区域。通过对不同阶段应力应变分布云图的详细分析，可以发现应力和应变的分布与褶皱构造的发育密切相关。在褶皱形成的初期，应力集中在边界和薄弱部位，为褶皱的起始提供了动力条件；随着褶皱的发展，应力和应变逐渐向褶皱区域集中，使得褶皱不断发育和演化。在褶皱形成的后期，应力和应变在褶皱和断裂部位达到相对稳定的状态，但仍然保持着较高的值，这反映了褶皱和断裂构造的稳定性和持续性。5.1.2褶皱形态演化过程为了直观地展示川东褶皱带褶皱形态的演化过程，制作了一系列图片和动画。在模拟的起始阶段，模型呈现为水平的地层结构，没有明显的褶皱形态。随着模拟的开始，水平挤压应力逐渐施加，地层开始发生微小的变形，在模型的局部区域出现了一些轻微的弯曲，这是褶皱形成的初始迹象。这些微小的弯曲主要出现在岩石力学性质相对较弱的部位，如页岩层或存在节理、裂隙的区域，因为这些部位更容易受到应力的作用而发生变形。随着挤压应力的持续作用，褶皱开始逐渐发育。在隔档式褶皱区域，向斜和背斜的雏形开始显现，向斜部位的地层逐渐下凹，背斜部位的地层逐渐上凸。此时，褶皱的形态较为简单，轴向也相对较为平直。在隔槽式褶皱区域，背斜和向斜的发育特征与隔档式褶皱有所不同，背斜部位的地层开始明显展宽，向斜部位的地层相对变窄，褶皱的形态呈现出背斜宽、向斜窄的特征。随着模拟的进一步进行，褶皱的形态逐渐变得复杂。在隔档式褶皱区域，向斜和背斜的规模不断增大，轴面开始发生倾斜，褶皱的形态逐渐从简单的箱状向更为复杂的形态转变。在隔槽式褶皱区域，背斜和向斜的边界更加清晰，背斜的顶部逐渐变得平缓，向斜的底部则更加狭窄，褶皱的形态更加典型。在褶皱形态演化的过程中，还可以观察到褶皱之间的相互作用。相邻的褶皱在发育过程中，会相互影响和干扰，导致褶皱的形态和轴向发生变化。一些褶皱可能会发生合并或分叉现象，使得褶皱的结构更加复杂。在模拟的后期，褶皱形态逐渐趋于稳定，形成了与实际地质情况相似的隔档式和隔槽式褶皱构造格局。隔档式褶皱呈现出向斜宽、背斜窄，轴面多向西倾的特征；隔槽式褶皱呈现出背斜宽、向斜窄，轴面多向南东倾的特征。通过对褶皱形态演化过程的详细观察和分析，可以清晰地了解川东褶皱带在中生代的构造演化历史，以及不同地质因素对褶皱形成和发展的影响。5.2结果分析与讨论5.2.1与实际地质构造的对比验证将数值模拟结果与实际地质构造特征进行对比，是验证模拟准确性和可靠性的关键步骤。在褶皱形态方面，模拟得到的隔档式褶皱和隔槽式褶皱的形态与实际地质情况具有较高的相似度。模拟的隔档式褶皱呈现出向斜宽、背斜窄的特征，轴面多向西倾，这与实际地质中川东褶皱带西带的隔档式褶皱形态一致。华蓥山背斜作为川东隔档式褶皱的典型代表，其在模拟中的形态特征，如核部地层的出露、轴面的倾斜方向以及背斜和向斜的宽度比例等，都与实际地质观测结果相符。模拟的隔槽式褶皱呈现出背斜宽、向斜窄的特征，轴面多向南东倾，与实际地质中川东褶皱带东带的隔槽式褶皱形态相吻合。桑植-石门复向斜在模拟中的形态和构造特征，也与实际地质情况较为一致，其褶皱的复杂性和组合样式在模拟结果中得到了较好的体现。在断裂构造方面，模拟结果同样与实际地质构造具有较好的一致性。模拟的齐岳山断裂在平面上呈S形分布，由多条不连续断层组成，北部呈右阶左行雁列式，南部左阶右行排列，这与实际地质中齐岳山断裂的形态和组合特征相符。模拟的华蓥山断裂走向北北东，倾向东，倾角较陡，力学性质为压扭性，这与实际地质中华蓥山断裂的特征一致。并且，模拟结果能够较好地反映齐岳山断裂和华蓥山断裂对褶皱构造的控制作用。在模拟中，齐岳山断裂控制形成了西侧背斜和西侧斜坡两类转换带，并分别发育“基底冲断一断层相关褶皱一滑脱”和“基底逆冲一滑脱”两类构造变形模式，这与实际地质中的构造变形模式相符合。华蓥山断裂对川东褶皱带西带隔档式褶皱的控制作用在模拟中也得到了体现，如华蓥山背斜核部出露古生界地层，且轴面在华蓥山断裂带附近向东倾，模拟结果与实际地质情况一致。通过对比模拟结果与实际地质构造特征，可以发现模拟结果在褶皱形态和断裂构造等方面都能够较好地再现实际地质情况。这表明本次数值模拟所采用的模型、参数设置以及模拟方法是合理的，能够较为准确地反映中生代川东褶皱带的构造演化过程，为进一步研究川东褶皱带的形成机制和演化历史提供了可靠的依据。5.2.2影响褶皱形成的因素探讨岩石力学性质是影响褶皱形成的重要因素之一。不同岩石的弹性模量、泊松比、粘聚力和内摩擦角等力学参数差异，导致它们在构造应力作用下的变形行为各不相同。弹性模量较高的岩石，如灰岩，在受到水平挤压应力时，抵抗变形的能力较强，更容易发生脆性变形，形成断裂；而弹性模量较低的岩石，如页岩，抵抗变形的能力较弱，更容易发生塑性变形，常作为滑脱层，促进上覆地层的褶皱形成。在川东褶皱带中，海相盖层中的碳酸盐岩和碎屑岩由于力学性质的差异，在褶皱形成过程中扮演了不同的角色。碳酸盐岩的脆性变形使得褶皱的形态更加尖锐，而碎屑岩的塑性变形则使得褶皱的形态更加平缓。粘聚力和内摩擦角也对褶皱的形成和稳定性产生重要影响。粘聚力较大的岩石，在受力时能够保持较好的整体性，有利于褶皱的形成和维持；而内摩擦角较大的岩石，则在剪切作用下更难发生滑动，影响褶皱的发育和演化。构造应力的大小、方向和作用方式对褶皱的形成和演化起着关键作用。在川东褶皱带的模拟中，水平挤压应力是导致褶皱形成的主要动力。随着水平挤压应力的增大，褶皱的幅度和波长逐渐增大，褶皱的形态也更加复杂。当水平挤压应力超过岩石的屈服强度时，岩石会发生塑性变形，形成褶皱和断裂。构造应力的方向也会影响褶皱的轴向和轴面的倾向。在川东褶皱带中，由于受到来自北北东—北东向的构造应力作用，褶皱的轴向总体呈北北东—北东向，轴面的倾向也与构造应力的方向相关。隔档式褶皱的轴面多向西倾，隔槽式褶皱的轴面多向南东倾，这与构造应力的作用方向和地层的力学性质有关。构造应力的作用方式，如持续加载或脉冲加载，也会对褶皱的形成产生影响。持续加载会使得褶皱逐渐发育和演化，而脉冲加载则可能导致褶皱的突然变形和调整。滑脱层在褶皱形成过程中具有重要作用，它能够降低上覆地层与下伏地层之间的摩擦力，使得上覆地层更容易发生滑动和变形，从而促进褶皱的形成。在川东褶皱带中，不同深度的滑脱层对褶皱的形态和分布产生了显著影响。一般认为，川东构造带西部的隔档式褶皱主要受志留系滑脱层控制，而东部的隔槽式褶皱则受海相盖层与变质基底之间的滑脱层控制。模拟结果表明，当滑脱层位于合适的深度时，能够有效地促进褶皱的形成和发育。在隔档式褶皱区域，志留系滑脱层的存在使得上覆地层在水平挤压应力作用下更容易发生滑动和弯曲，形成向斜宽、背斜窄的褶皱形态；在隔槽式褶皱区域，海相盖层与变质基底之间的滑脱层使得上覆地层能够更自由地变形，形成背斜宽、向斜窄的褶皱形态。滑脱层的厚度和连续性也会影响褶皱的形成。较厚且连续的滑脱层能够提供更好的滑动条件，促进褶皱的发育；而较薄或不连续的滑脱层则可能限制褶皱的形成和发展。5.2.3两阶段构造演化模型分析结合模拟结果，两阶段构造演化模型能够合理地解释川东褶皱带的形成过程。在第一阶段，隔槽褶皱区首先形成。这一阶段，由于区域构造应力场的作用，以及海相盖层与变质基底之间滑脱层的影响，使得该区域的地层发生了强烈的变形。在水平挤压应力的作用下，地层沿着滑脱层发生滑动和弯曲，形成了背斜宽、向斜窄的隔槽式褶皱。这一过程中，岩石的力学性质和构造应力的大小、方向等因素共同作用，决定了隔槽式褶皱的形态和特征。随着构造演化的进行，进入第二阶段，隔档褶皱区开始形成。此时，区域构造应力场可能发生了一定的变化，或者由于第一阶段隔槽褶皱区形成后，应力场发生了重新分布。在新的应力场作用下，以及志留系滑脱层的控制作用，使得隔档褶皱区的地层发生变形，形成了向斜宽、背斜窄的隔档式褶皱。在这一阶段，齐岳山断裂的雏形开始形成，它在晚侏罗世主滑脱层抬升导致隔槽式褶皱带向隔档式转变过程中起到了重要作用。两阶段构造演化模型与模拟结果和实际地质构造特征相符合。模拟结果清晰地展示了隔槽褶皱区和隔档褶皱区依次形成的过程，以及褶皱形态和断裂构造的演化特征。在实际地质构造中，也可以观察到隔槽式褶皱和隔档式褶皱在空间上的分布特征以及它们之间的过渡关系，这与两阶段构造演化模型的预测一致。该模型还能够解释川东褶皱带中一些复杂的构造现象，如褶皱的叠加和变形、断裂的形成和演化等。在两阶段构造演化过程中，不同阶段形成的褶皱和断裂相互作用，导致了构造的复杂性。通过对两阶段构造演化模型的分析，可以更好地理解川东褶皱带的形成机制和演化历史，为进一步研究该地区的地质构造提供了重要的理论框架。六、结论与展望6.1研究主要结论通过对中生代川东褶皱带的数值模拟研究，本研究取得了以下主要结论：构造演化过程与特征：清晰揭示了川东褶皱带在中生代的构造演化过程，明确其经历了两阶段的构造变形。第一阶段，隔槽褶皱区在区域构造应力场及海相盖层与变质基底之间滑脱层的共同作用下首先形成，呈现出背斜宽、向斜窄的典型特征；第二阶段，隔档褶皱区在新的应力场和志留系滑脱层的控制下逐渐形成，表现为向斜宽、背斜窄的构造样式。模拟结果精确展示了褶皱形态从初始的微小弯曲逐渐发展为复杂褶皱构造的过程，与实际地质构造中的褶皱形态演化高度契合。在模拟过程中，早期地层在水平挤压应力作用下出现轻微弯曲，随着应力持续作用，隔槽褶皱区和隔档褶皱区的褶皱逐渐发育成熟，其形态、规模和轴向变化都与实际地质观测结果一致。应力应变与褶皱形成机制：深入分析了应力应变分布与褶皱形成的内在联系。在构造演化过程中，应力应变分布呈现出明显的规律性，且与褶皱构造的发育紧密相关。在褶皱形成初期，应力主要集中在边界和岩石力学性质相对薄弱的部位，为褶皱的起始提供了必要的动力条件。随着褶皱的发展，应力和应变逐渐向褶皱区域集中，促使褶皱不断发育和演化。在褶皱形成的后期，应力和应变在褶皱和断裂部位达到相对稳定状态，但仍保持较高值，这反映了褶皱和断裂构造的稳定性和持续性。通过对应力应变分布云图的详细分析，清晰地观察到了应力集中区域的迁移和变化，以及应变在不同阶段的扩展和分布特征，进一步验证了应力应变与褶皱形成的密切关系。关键影响因素：系统探讨了岩石力学性质、构造应力和滑脱层等因素对褶皱形成的重要影响。不同岩石的力学参数差异决定了其在构造应力作用下的变形行为，弹性模量较高的岩石易发生脆性变形，弹性模量较低的岩石易发生塑性变形。构造应力的大小、方向和作用方式对褶皱的形成和演化起着决定性作用，水平挤压应力是导致褶皱形成的主要动力，其大小和方向的变化直接影响褶皱的幅度、波长、轴向和轴面倾向。滑脱层能够有效降低上覆地层与下伏地层之间的摩擦力，使上覆地层更容易发生滑动和变形，从而促进褶皱的形成。在川东褶皱带中，不同深度的滑脱层对褶皱的形态和分布产生了显著影响，西部隔档式褶皱主要受志留系滑脱层控制，东部隔槽式褶皱主要受海相盖层与变质基底之间的滑脱层控制。模拟结果验证：将数值模拟结果与实际地质构造进行了全面对比验证，结果表明模拟结果在褶皱形态和断裂构造等方面与实际地质情况具有高度的一致性。模拟得到的隔档式褶皱和隔槽式褶皱的形态、轴面倾向以及断裂的走向、形态和组合特征等，都与实际地质观测结果相符，充分验证了本次数值模拟所采用的模型、参数设置以及模拟方法的合理性和准确性，为进一步研究川东褶皱带的形成机制和演化历史提供了坚实