基于S波接收函数的四川地区岩石圈埋深解析与动力学机制探究

基于S波接收函数的四川地区岩石圈埋深解析与动力学机制探究一、引言1.1研究背景与意义四川地区地处中国大陆的西南部，是中国地质构造最为复杂的区域之一。它不仅是扬子克拉通西缘的关键组成部分，还受到印度板块与欧亚板块强烈碰撞的远程效应影响，经历了多期次的构造运动和演化，造就了现今复杂多样的地质构造格局。在板块碰撞的持续作用下，印度板块向北推挤，使得青藏高原不断隆升，并向东挤出物质，四川地区首当其冲，成为了构造应力的强烈作用区域。这种强烈的构造活动导致了四川地区岩石圈结构的显著变化，岩石圈的厚度、物质组成和力学性质在不同区域存在明显差异。而岩石圈作为地球最外层的刚性固体圈层，它的结构和演化对地球的动力学过程有着至关重要的影响，其不仅控制着板块的运动和相互作用，还与地震、火山等地质灾害的发生密切相关。因此，深入研究四川地区岩石圈的结构和演化，对于理解中国大陆的构造演化、地球动力学过程以及地质灾害的发生机制都具有十分重要的意义。四川地区的岩石圈结构研究在地球科学领域占据着关键地位。它处于多个板块和构造单元的交汇地带，扬子克拉通西缘的古老基底与青藏高原东缘的新生代构造变形带在此相互作用。这种独特的地质位置使得四川地区成为研究岩石圈演化和板块相互作用的天然实验室。通过对该地区岩石圈结构的研究，我们可以深入了解印度板块与欧亚板块碰撞远程效应下，岩石圈的变形机制和响应过程，这对于揭示中国大陆的构造演化历史、重建古板块运动轨迹以及理解地球内部动力学过程都提供了重要的线索和依据。此外，四川地区频发的地震灾害也使得岩石圈结构研究变得更为紧迫。例如，2008年的汶川大地震和2013年的芦山地震，都给当地人民的生命财产造成了巨大损失。这些地震的发生与岩石圈的结构和深部动力学过程密切相关，深入研究岩石圈结构有助于我们更好地理解地震的孕育和发生机制，提高地震预测和灾害防治的能力。在岩石圈结构的研究方法中，S波接收函数方法凭借其独特的优势脱颖而出。当地震波在地球内部传播时，会在不同介质的界面上发生反射、折射和转换等现象，而S波接收函数正是利用这些特性，通过对地震台站记录的远震S波波形进行处理和分析，有效提取出地球内部不同深度界面的信息。相较于其他方法，S波接收函数方法具有较高的分辨率，能够清晰地分辨出岩石圈内部的精细结构，如莫霍面、岩石圈-软流圈边界等重要界面的深度和起伏变化，这为研究岩石圈的厚度和结构提供了精确的数据支持。同时，该方法还能对岩石圈内部的速度结构变化进行有效探测，通过分析转换波的到时和振幅等信息，获取岩石圈内部不同深度的速度变化情况，从而推断岩石圈的物质组成和力学性质。而且，S波接收函数方法不受地形和地质条件的限制，无论是在山区、平原还是海洋等不同地形条件下，都能够有效地获取地震数据并进行分析，具有广泛的适用性和应用前景。在四川地区复杂的地质条件下，S波接收函数方法能够充分发挥其优势，为研究该地区岩石圈结构提供了强有力的技术手段，有助于我们深入了解岩石圈的深部结构和动力学特征。1.2国内外研究现状在岩石圈结构的研究领域，S波接收函数方法凭借其独特优势，已成为国内外学者广泛应用的重要手段之一。国外学者在早期就开展了相关研究，例如，[国外学者名字1]利用S波接收函数对北美地区的岩石圈进行了深入研究，通过分析大量的地震数据，成功绘制出该地区岩石圈底部界面的起伏图像，揭示了岩石圈厚度在不同区域的变化特征，发现了岩石圈与软流圈之间存在复杂的相互作用关系。[国外学者名字2]则对欧洲部分地区进行研究，运用先进的数据处理技术和成像方法，从S波接收函数中提取出了岩石圈内部更为精细的结构信息，指出岩石圈内部存在多个速度变化界面，这些界面的存在对地震波的传播路径和速度产生了显著影响。国内学者也在该领域取得了一系列丰硕成果。张耀阳等人基于跨越华夏块体至四川盆地西部的130个线性流动地震台站及其附近90个固定台网台站的观测资料，采用S波接收函数波动方程叠后偏移方法，开展了华南大陆岩石圈结构研究。成像结果显示，研究区岩石圈结构复杂，不同构造单元之间差异显著，构造边界带附近小尺度变化强烈。150km以上的厚岩石圈主要位于四川盆地，不足100km的薄岩石圈主要分布于川东褶皱带和华夏块体，雪峰山下方岩石圈厚度显著增加，且以雪峰山为界岩石圈结构和性质存在着显著的东西差异。通过结合其它地球物理观测得到的地壳-上地幔结构信息，提出四川盆地还保留着厚而冷的克拉通岩石圈根，且岩石圈地幔具有结构分层特征；雪峰山可能是扬子克拉通与华夏块体在西南部的边界；雪峰山以东区域可能经历了岩石圈的减薄和改造，且华南岩石圈的减薄与华北相似，都主体发生在东部地区，造成现今南北重力梯度带两侧强烈的结构差异，该研究结果为认识华南大陆的构造演化及其深部动力学提供了地震学约束。在青藏高原及周边地区的研究中，南科大学者郭震和陈永顺课题组利用一条跨越青藏高原东北缘地区宽频带地震观测剖面所记录的远震体波波形资料，联合P波和S波接收函数方法，构建了青藏高原东北缘高分辨率的岩石圈界面结构图像。研究结果显示，岩石圈（LAB）界面以南祁连缝合带(SQS)为界，呈现出差异增厚特征，在已经发育完全的松潘-甘孜(SG)和中亚造山带(CAOB)下方岩石圈较薄，而正在生长的西秦岭(QKW)和祁连造山带(QL)和阿拉善(Alxa)地块下方岩石圈较厚，在南祁连缝合带（SQS）下方，存在明显的地壳和岩石圈界面错断（~10km），并认为南祁连缝合带是目前青藏岩石圈和亚洲岩石圈的汇聚边界，青藏岩石圈地幔与亚洲岩石圈地幔在南祁连缝合带的汇聚与南北两侧的岩石圈差异增厚主导着青藏高原东北缘新生代以来的隆升、扩展和构造变形，该研究首次明确了岩石圈变形对于高原隆升和扩展的重要性，为研究青藏高原东北缘的隆升及其构造变形机制提供了新的地震学约束。针对S波接收函数传统处理流程中常因人为设定参数而引入诸多不可控误差的问题，广州地球化学研究所的张周和邓阳凡研究员提出了一种通用的估计策略——GC_SRF。该策略通过网格搜索和相关性分析，可以精准确定SRF计算时的旋转角度和时窗长度两个主要参数，同时有效避免人为干扰的引入。通过对比分析，选用30%的均方根误差(RMSE)比例对SRF进行有效质量控制。理论测试表明，该方法能从全波场理论地震图中恢复转换于LAB的微弱地震信号，进一步结合偏移成像等方法能够获得清晰的岩石圈几何形态。将该策略应用于青藏高原东南缘，揭示了该区岩石圈几何结构强烈不均一性，结果显示四川盆地岩石圈厚度约为160km且在100km深度附近存在清晰的MLD，腾冲火山岩石圈厚度约为90km，而川滇地区下方的LAB却无法从成像结果中清晰分辨，研究认为川滇地区下方无法清晰识别的LAB指示了地幔柱作用下岩石圈改造后的重新生长特征，也正是这一复杂过程造成该地区复杂的地震波速结构、低衰减和高电导率等地球物理异常。尽管国内外学者在利用S波接收函数研究岩石圈结构方面取得了众多成果，但针对四川地区的研究仍存在一定的不足。四川地区地质构造复杂，受到多种构造运动的叠加影响，然而目前的研究在台站分布上还不够密集，导致对一些局部区域的岩石圈结构分辨率较低，无法准确刻画岩石圈内部的精细结构和深部动力学过程。在数据处理和解释方面，虽然已有多种方法和技术，但对于四川地区复杂地质条件下的数据处理，仍缺乏针对性和有效性，不同方法得到的结果之间存在一定的差异，使得对岩石圈结构的认识存在不确定性。此外，对于岩石圈结构与地质灾害之间的关系研究还不够深入，如何将岩石圈结构的研究成果应用于地震灾害的预测和防治，仍是亟待解决的问题。1.3研究内容与方法本研究旨在利用S波接收函数深入剖析四川地区岩石圈埋深及动力学机制，具体研究内容涵盖以下几个关键方面：首先，通过对四川地区多个地震台站记录的远震S波数据进行系统收集和整理，获取高质量的地震波形资料，为后续的分析提供坚实的数据基础。这些数据将成为我们窥探岩石圈内部结构的关键窗口，其准确性和完整性直接影响研究结果的可靠性。其次，运用先进的S波接收函数计算方法，对收集到的地震波形进行精确处理，提取出岩石圈内部不同深度界面的转换波信息。在这个过程中，我们将采用如脉冲反褶积、旋转坐标等一系列技术手段，以确保能够从复杂的地震波信号中有效分离出反映岩石圈结构的关键震相，为后续的成像和分析提供清晰、准确的信息。在获取转换波信息后，我们将进一步利用共转换点叠加（CCP）和波动方程偏移成像等技术，对岩石圈内部的结构进行详细成像。通过CCP叠加，可以将来自不同方位和距离的地震波转换点聚集在一起，增强信号的信噪比，从而更清晰地揭示岩石圈内部界面的起伏和变化。而波动方程偏移成像则能够根据波动理论，将地震波的传播路径和反射、折射信息进行精确模拟，构建出高分辨率的岩石圈结构图像，直观展示岩石圈的厚度变化、内部界面的形态以及可能存在的异常结构。通过这些成像技术，我们将全面了解四川地区岩石圈的三维结构特征，包括莫霍面、岩石圈-软流圈边界等重要界面的深度和形态变化，以及岩石圈内部可能存在的低速层、高速体等异常结构。这些信息将为深入研究岩石圈的动力学机制提供关键的约束条件，帮助我们揭示岩石圈的演化历史和现今的动力学状态。为了深入探讨四川地区岩石圈的动力学机制，我们将结合岩石圈结构的成像结果，综合考虑区域地质构造背景、地球物理场特征以及板块运动等因素。通过分析岩石圈结构与这些因素之间的内在联系，我们将尝试揭示岩石圈变形的驱动力来源、变形机制以及演化过程。例如，考虑到四川地区位于印度板块与欧亚板块碰撞的远程效应区域，我们将研究板块碰撞如何通过应力传递和物质迁移影响岩石圈的结构和动力学过程；同时，结合区域内的地震活动、地热流分布等地球物理场特征，探讨岩石圈内部的热-动力耦合作用对岩石圈演化的影响。此外，我们还将与周边地区的岩石圈结构和动力学特征进行对比分析，进一步明确四川地区在区域构造演化中的位置和作用，揭示其独特的动力学机制。在研究方法上，本研究将以S波接收函数方法为核心，并结合其他相关的地球物理方法和地质资料进行综合分析。S波接收函数方法作为研究岩石圈结构的重要手段，具有高分辨率、对深部结构敏感等优势，能够有效提取岩石圈内部的精细结构信息。在数据处理过程中，我们将严格遵循数据处理流程和质量控制标准，确保数据的准确性和可靠性。同时，为了提高研究结果的精度和可靠性，我们将采用多种方法进行交叉验证和对比分析。例如，将S波接收函数成像结果与P波接收函数成像结果、面波层析成像结果以及重力、磁力等地球物理资料进行对比，综合分析不同方法得到的结果之间的一致性和差异性，从而更准确地解释岩石圈结构和动力学特征。此外，我们还将结合地质构造演化历史、岩石学和地球化学等方面的研究成果，从多个角度深入探讨岩石圈的动力学机制，为全面理解四川地区的地质构造演化提供有力的支持。二、S波接收函数的基本原理与方法2.1S波接收函数的理论基础S波接收函数作为研究地球内部结构的重要工具，其理论基础根植于地震波传播理论和地球介质的弹性性质。当地震波在地球内部传播时，由于地球介质并非均匀连续，而是存在着不同密度和弹性参数的界面，地震波在这些界面上会发生复杂的传播行为，包括反射、折射和转换等现象，这些现象蕴含着丰富的地球内部结构信息，而S波接收函数正是基于对这些现象的深入研究和巧妙利用而发展起来的。在地震波传播过程中，当S波入射到地球内部不同介质的界面时，会产生一系列复杂的波场响应。假设地球内部存在一个简单的两层模型，上层介质的密度为\rho_1，剪切波速度为v_{s1}，下层介质的密度为\rho_2，剪切波速度为v_{s2}。当S波以入射角\theta从上层介质入射到界面时，根据弹性波理论，会在界面上产生反射S波和透射S波，同时还会产生转换P波，这些波的传播方向和振幅大小都与界面两侧介质的性质以及入射角密切相关。根据广义反射透射系数矩阵方法，可以精确计算出各个波的反射系数和透射系数，进而得到它们的振幅和相位信息。这些计算结果表明，不同震相的波到达观测台站的时间和振幅特征，能够反映出界面的深度、介质速度的变化以及界面的起伏情况等重要信息。S波接收函数的提取过程，本质上是从复杂的地震波记录中分离出与地下界面转换波相关的信息。在实际地震观测中，地震台站记录到的地震波是多种波的叠加，包括直达波、反射波、转换波以及各种多次波等，这些波相互干涉，使得原始地震记录十分复杂。为了提取出S波接收函数，需要采用一系列的数据处理技术和方法。首先，利用等效震源的概念，假设观测台站下方介质对地震波垂直分量的响应是一个\delta函数，这一假设基于深源远震地表位移的垂直分量表现为尖脉冲的时间函数与仪器响应的褶积，且后续续至震相较小的观测事实，以及理论上陡角度入射P波所产生的转换波及地壳内部鸣震震相的垂直分量较小的计算结果。在此假设下，将谱振幅比值修改为频率比值，并通过反变换将其转换到时间域，从而计算出观测台站下方介质对等效震源在水平分量上的响应，即得到S波接收函数。在实际操作中，通常采用三分量地震仪记录地震波，通过对垂直分量和水平分量的地震记录进行处理，具体而言，用垂直分量的地震记录对水平分量的地震记录进行反褶积处理，在频率域中，将垂直分量地震记录的频谱作为分母，水平分量地震记录的频谱作为分子，得到频率域的响应函数，再将其反变换回时间域，就可以得到介质结构响应的径向分量和切向分量，也就是我们所需要的径向接收函数和切向接收函数。通过这样的处理，最大程度地消除了震源时间函数及传播路径对远震体波波形的影响，保留了相位信息，使得时间域的接收函数能够更直观地反映地下界面的特征。S波接收函数对岩石圈内部结构信息的反映是多方面且深入的。岩石圈作为地球最外层的刚性固体圈层，其内部存在着多个重要的界面，如莫霍面（地壳与上地幔的分界面）、岩石圈-软流圈边界等，这些界面两侧的介质性质存在明显差异，当S波传播到这些界面时，会产生强烈的转换波信号，这些信号被记录在S波接收函数中，成为我们探测岩石圈内部结构的关键线索。例如，莫霍面作为地壳和上地幔的分界面，其深度的变化反映了地壳厚度的变化，而S波接收函数中来自莫霍面的转换波震相的到时，可以精确计算出莫霍面的深度。通过对不同台站S波接收函数的分析，可以绘制出莫霍面的起伏图像，揭示地壳厚度在不同区域的变化特征。对于岩石圈-软流圈边界，由于岩石圈和软流圈的物质组成和力学性质存在显著差异，S波在该边界上也会产生明显的转换波信号，通过分析S波接收函数中与该边界相关的转换波特征，如振幅、到时等，可以推断岩石圈-软流圈边界的深度和形态变化，了解岩石圈的厚度分布情况。此外，S波接收函数还能够对岩石圈内部的速度结构变化进行有效探测，通过分析转换波的到时和振幅等信息，可以获取岩石圈内部不同深度的速度变化情况，进而推断岩石圈的物质组成和力学性质。例如，在岩石圈内部，如果存在低速层或高速体，这些异常结构会导致S波接收函数中转换波的到时和振幅发生异常变化，通过对这些异常变化的分析，可以识别出岩石圈内部的异常结构，并探讨其形成机制和地质意义。2.2数据处理流程本研究的数据收集主要来源于四川地区多个地震台站的观测记录，这些台站分布广泛，涵盖了四川盆地、川西高原、川滇边界等不同地质构造区域，能够全面地记录该地区的地震活动。台站所使用的地震仪器均为高精度的宽频带地震仪，能够精确地捕捉到地震波的各种信息，包括地震波的振幅、频率、相位等，为后续的研究提供了高质量的数据基础。数据收集的时间跨度较长，以确保能够获取足够数量的远震事件，从而提高研究结果的可靠性和代表性。在数据收集过程中，严格按照相关的观测规范和标准进行操作，确保数据的准确性和完整性。对数据进行实时监测和质量控制，及时发现和处理可能出现的数据异常情况，如仪器故障、噪声干扰等，以保证收集到的数据能够真实地反映地震波的传播特征。数据预处理是S波接收函数计算的关键环节，其目的在于去除数据中的噪声和干扰，提高数据的质量和信噪比，为后续的计算和分析提供可靠的数据基础。在数据预处理过程中，首先进行仪器响应校正，由于地震仪在记录地震波时，会对地震波的波形和频率产生一定的影响，导致记录到的数据与实际的地震波存在差异，因此需要通过仪器响应校正来消除这种影响，恢复地震波的原始特征。根据地震仪的型号和参数，获取相应的仪器响应文件，利用专业的数据处理软件，对地震记录进行仪器响应校正，将记录到的地震数据从仪器响应的影响中解放出来，得到真实的地震波信号。去除噪声和干扰是数据预处理的重要步骤，地震数据中常常包含各种噪声和干扰，如环境噪声、电磁干扰、仪器噪声等，这些噪声和干扰会严重影响S波接收函数的计算结果，降低数据的信噪比。为了有效地去除这些噪声和干扰，采用了多种滤波技术，如带通滤波、高通滤波、低通滤波等。根据地震波的频率特征和噪声的频率分布，选择合适的滤波参数，设计滤波器对地震数据进行滤波处理。通过带通滤波，可以保留地震波的有效频率成分，去除高频和低频噪声；高通滤波可以去除低频噪声和干扰，突出地震波的高频信号；低通滤波则可以去除高频噪声和干扰，保留地震波的低频成分。在实际应用中，通常会根据数据的特点和噪声的情况，综合运用多种滤波技术，以达到最佳的去噪效果。除了滤波技术，还采用了一些其他的去噪方法，如自适应噪声对消、小波变换去噪等，进一步提高数据的质量。在完成仪器响应校正和去噪处理后，还需要对地震数据进行重采样和时窗选取。由于不同台站的地震数据采样率可能不同，为了便于后续的数据处理和分析，需要将所有数据重采样到相同的采样率。根据研究的需要和数据的特点，选择合适的重采样率，利用插值算法对数据进行重采样，使所有数据具有相同的时间分辨率。时窗选取是为了截取地震波信号中与S波接收函数计算相关的部分，去除其他无关的信号，提高计算效率和准确性。根据地震波的传播时间和S波接收函数的特征，确定合适的时窗长度和位置，在地震记录中选取包含S波震相的时窗，将时窗内的数据用于后续的S波接收函数计算。在时窗选取过程中，需要注意避免选取到噪声和干扰较大的部分，同时要确保时窗内包含足够的S波信号信息，以保证S波接收函数计算的准确性。S波接收函数的计算是本研究的核心步骤，其计算方法主要基于地震波传播理论和信号处理技术。在计算过程中，首先利用三分量地震记录，将垂直分量作为参考，对水平分量进行反褶积处理，以提取出S波接收函数。具体而言，将垂直分量的地震记录作为等效震源，假设台站下方介质对地震波垂直分量的响应为一个\delta函数，这是基于深源远震地表位移垂直分量表现为尖脉冲与仪器响应褶积且续至震相较小，以及理论上陡角度入射P波产生的转换波及地壳内部鸣震震相垂直分量较小的事实。在此假设下，将水平分量的地震记录与垂直分量进行反褶积运算，在频率域中，将垂直分量地震记录的频谱作为分母，水平分量地震记录的频谱作为分子，得到频率域的响应函数，再通过傅里叶反变换将其转换回时间域，从而得到介质结构响应的径向分量和切向分量，即径向接收函数和切向接收函数。在反褶积过程中，为了提高计算结果的稳定性和准确性，通常会采用一些反褶积算法，如脉冲反褶积、最小熵反褶积等，并根据数据的特点和计算要求选择合适的参数。例如，在脉冲反褶积中，需要设置合适的脉冲宽度和白噪水平等参数，以确保能够有效地提取出S波接收函数。为了提高S波接收函数的质量和可靠性，还可以采用叠加技术对多个接收函数进行叠加处理。由于不同地震事件的接收函数可能存在一定的噪声和误差，通过叠加多个接收函数，可以有效地降低噪声的影响，提高信号的信噪比，使接收函数中的有效信号更加突出。在叠加过程中，通常会根据地震事件的震中距、方位角等参数对接收函数进行分类和加权，对于震中距和方位角相近的接收函数给予较大的权重，以提高叠加效果。同时，还可以采用一些统计方法对叠加结果进行质量评估，如计算均方根误差、相关系数等，以确保叠加后的接收函数具有较高的质量和可靠性。质量控制是确保S波接收函数数据可靠性和有效性的重要环节，通过一系列严格的质量控制措施，可以有效筛选出高质量的接收函数，为后续的研究提供可靠的数据支持。在质量控制过程中，首先对S波接收函数的信噪比进行评估，信噪比是衡量接收函数质量的重要指标，它反映了信号与噪声的相对强度。通过计算接收函数中信号部分的能量与噪声部分的能量之比，得到信噪比。设定一个合适的信噪比阈值，对于信噪比低于阈值的接收函数，认为其质量较差，可能受到噪声的严重干扰，将其剔除。在实际应用中，通常会根据数据的特点和研究的要求，确定合适的信噪比阈值，一般来说，信噪比阈值可以设置在3-5之间，具体数值需要根据实际情况进行调整。除了信噪比评估，还对接收函数的波形特征进行检查，一个高质量的接收函数应该具有清晰的波形特征，能够准确反映地下界面的信息。检查接收函数的波形是否存在异常，如波形的突变、振荡等，这些异常可能是由于噪声干扰、数据处理错误或地下结构异常等原因引起的。对于波形特征异常的接收函数，进行仔细分析和排查，找出原因并进行相应的处理。如果是噪声干扰引起的，可以尝试采用更有效的去噪方法进行处理；如果是数据处理错误导致的，需要重新检查数据处理过程，修正错误；如果是地下结构异常导致的，则需要结合其他地球物理资料进行深入分析，以确定其地质意义。还通过对比不同台站的接收函数，来检查接收函数的一致性和可靠性。由于不同台站的地质条件和地震波传播路径可能存在差异，因此不同台站的接收函数应该具有一定的相似性，但也会存在一些差异。如果发现某些台站的接收函数与其他台站的接收函数存在明显的不一致，如震相的到时、振幅、波形等方面存在较大差异，需要对这些台站的数据进行重新检查和分析，找出原因并进行处理。可能是由于台站的仪器故障、数据处理错误或地质条件的特殊性等原因导致的，通过进一步的调查和分析，可以确定问题的根源，并采取相应的措施进行解决，以确保不同台站的接收函数具有较高的一致性和可靠性。2.3成像方法与技术成像原理是基于S波接收函数中转换波的到时信息，通过合理假设地球内部介质的速度模型，利用地震波传播理论，将时间域的接收函数转换为深度域的图像，从而清晰展现岩石圈内部不同深度界面的位置和形态。在实际操作中，通常假设地球介质为层状结构，根据不同层间的速度差异以及S波在界面上的转换关系，运用射线理论或波动方程理论来计算地震波的传播路径和到时。射线理论假设地震波沿直线传播，通过计算射线在不同介质层中的传播时间和折射、反射角度，来确定转换波的到时和传播路径，进而实现对岩石圈内部结构的成像；波动方程理论则更全面地考虑了地震波的波动特性，通过求解波动方程，能够更精确地模拟地震波在复杂介质中的传播过程，包括波的干涉、绕射等现象，从而获得更准确的成像结果。在岩石圈结构成像中，常用的成像方法包括共转换点叠加（CCP）和波动方程偏移成像等，它们各自具有独特的特点和适用范围。共转换点叠加（CCP）是一种基于射线理论的成像方法，其原理是将来自不同方位和距离的地震波转换点聚集在同一深度位置上进行叠加。在实际地震观测中，由于地震波的传播路径和入射角不同，同一深度界面的转换波会在不同的时间和位置被记录到。CCP方法通过对地震波的射线路径进行计算和分析，将具有相同转换点深度的接收函数进行叠加，这样可以有效增强信号的信噪比，突出岩石圈内部界面的反射信息，使得原本微弱的信号得到增强，从而更清晰地揭示岩石圈内部界面的起伏和变化。该方法在岩石圈结构相对简单、横向变化较小的区域具有较好的成像效果，能够快速有效地获取岩石圈内部主要界面的大致位置和形态信息。例如，在一些稳定的克拉通地区，岩石圈结构相对均一，CCP方法能够清晰地成像出莫霍面和岩石圈-软流圈边界等主要界面，为研究该地区的岩石圈厚度和结构提供了重要的数据支持。波动方程偏移成像则是基于波动方程理论的成像方法，它充分考虑了地震波的波动特性，能够更精确地处理地震波在复杂介质中的传播过程。该方法通过对地震波场进行数值模拟，将观测到的地震数据反向传播回地下，根据波动方程的解来确定地下界面的位置和反射系数。在波动方程偏移成像中，常用的算法包括克希霍夫积分偏移、有限差分偏移和有限元偏移等。克希霍夫积分偏移是一种基于射线理论的积分方法，它通过对地震波的传播路径进行积分，来计算地下界面的反射响应；有限差分偏移则是将波动方程离散化，通过差分计算来求解地震波场；有限元偏移则是利用有限元方法对地下介质进行离散化，求解波动方程得到地震波场。这些算法各有优缺点，在实际应用中需要根据具体情况选择合适的算法。波动方程偏移成像能够适应复杂的地质构造，对于岩石圈内部存在复杂的速度变化、断层、褶皱等地质构造的区域，能够提供更准确的成像结果，清晰地展现岩石圈内部的精细结构和复杂构造特征。例如，在构造活动强烈的造山带地区，岩石圈结构复杂，存在大量的断层和褶皱，波动方程偏移成像能够准确地成像出这些复杂构造，为研究该地区的岩石圈变形机制和动力学过程提供了关键的信息。在四川地区岩石圈研究中，不同成像方法的适用性需要根据该地区复杂的地质构造特点来综合考虑。四川地区地质构造复杂多样，既有稳定的四川盆地，又有构造活动强烈的川西高原和川滇边界地区。对于四川盆地，岩石圈结构相对较为稳定，横向变化较小，CCP方法能够较好地成像出岩石圈内部的主要界面，如莫霍面和岩石圈-软流圈边界，为研究盆地的岩石圈厚度和结构提供了有效的手段。然而，在川西高原和川滇边界地区，由于受到印度板块与欧亚板块碰撞的强烈影响，岩石圈结构复杂，存在大量的断层、褶皱以及速度异常区域，此时CCP方法的局限性就会凸显出来，难以准确地反映出岩石圈内部的精细结构和复杂构造。而波动方程偏移成像方法则能够充分发挥其优势，通过考虑地震波的波动特性和复杂介质的影响，准确地成像出这些复杂构造区域的岩石圈结构，为研究该地区的岩石圈变形机制和动力学过程提供了更精确的信息。因此，在四川地区岩石圈研究中，为了获得更全面、准确的岩石圈结构信息，通常会结合使用CCP和波动方程偏移成像等多种成像方法，相互补充和验证，以提高研究结果的可靠性和精度。三、四川地区地质构造背景3.1区域构造特征四川地区位于中国大陆西南部，地处扬子板块西缘，是青藏高原东缘的重要组成部分。其大地构造位置独特，处于特提斯-喜马拉雅构造域与滨太平洋构造域的结合部位，这使得该地区受到了多种构造运动的强烈影响，形成了现今复杂多样的地质构造格局。从板块构造的角度来看，印度板块持续向北推挤欧亚板块，这种强烈的碰撞作用不仅导致了青藏高原的隆升，还使得高原物质向东挤出，四川地区首当其冲，成为了构造应力的集中释放区域。在这种远程效应的影响下，四川地区的岩石圈发生了显著的变形和调整，形成了一系列规模宏大的断裂带和褶皱构造。四川地区主要的构造单元包括四川盆地、川西高原和川滇边界构造带。四川盆地作为扬子板块的重要组成部分，具有稳定的刚性基底，其基底由前震旦纪变质岩系组成，经历了多期构造运动的改造，形成了相对稳定的构造格局。在漫长的地质历史时期中，四川盆地经历了海相沉积、陆相沉积等多个阶段，沉积了巨厚的地层，这些地层记录了盆地的演化历史和构造运动的影响。盆地内部构造相对简单，地层平缓，主要发育一些小型的褶皱和断裂构造，但在盆地边缘，由于受到周边构造单元的挤压作用，构造变形较为强烈，形成了一系列的逆冲断层和褶皱带，如龙门山断裂带、华蓥山断裂带等。川西高原位于青藏高原东缘，是青藏高原向四川盆地过渡的地带。该地区构造活动强烈，岩石圈变形显著，主要受到印度板块与欧亚板块碰撞的直接影响。川西高原内部分布着多条大型断裂带，如鲜水河断裂带、安宁河断裂带、龙门山断裂带等，这些断裂带控制了该地区的构造格局和地震活动。鲜水河断裂带是一条左旋走滑断裂，历史上曾发生多次强烈地震，如1973年的炉霍7.6级地震，其地震活动频繁，对当地的地质构造和地貌演化产生了重要影响。安宁河断裂带则是一条兼具逆冲和走滑性质的断裂，它控制了安宁河谷的形成和演化，也是该地区地震活动的主要发震构造之一。龙门山断裂带更是备受关注，2008年的汶川8.0级特大地震就发生在该断裂带上，此次地震释放出巨大的能量，造成了极其严重的人员伤亡和财产损失，也使得龙门山断裂带成为全球研究的焦点。川西高原的地层以浅变质岩系和火山岩系为主，这些岩石经历了强烈的构造变形和变质作用，形成了复杂的褶皱和断裂构造，其岩石变形特征和构造样式反映了该地区强烈的构造活动和复杂的动力学过程。川滇边界构造带是四川地区与云南地区的交界地带，也是多个构造单元的交汇区域。该构造带受到印度板块与欧亚板块碰撞以及太平洋板块向欧亚板块俯冲的共同影响，构造运动十分复杂。川滇边界构造带内分布着一系列南北向的断裂带，如小江断裂带、则木河断裂带等，这些断裂带构成了川滇菱形块体的边界，控制了块体的运动和变形。小江断裂带是一条活动性很强的左旋走滑断裂，历史上多次发生强烈地震，如1833年的嵩明8.0级地震，其地震活动对周边地区的地质构造和生态环境造成了严重影响。则木河断裂带同样具有强烈的活动性，它与小江断裂带相互作用，共同影响着川滇边界地区的构造演化和地震活动。在川滇边界构造带，地层岩性复杂多样，既有古老的变质岩系，也有新生代的沉积岩和火山岩，这些不同时代和类型的岩石在强烈的构造作用下相互叠置，形成了独特的地质构造景观。3.2地层与岩石分布四川地区地层分布广泛，经历了漫长而复杂的地质演化过程，形成了丰富多样的地层序列。从太古宙到新生代，各个地质时期的地层在四川地区均有不同程度的出露，这些地层记录了该地区地质历史时期的构造运动、沉积环境变迁以及生物演化等重要信息，为研究岩石圈结构提供了重要的地质依据。太古宙和元古宙地层主要出露于四川盆地边缘和川西高原部分地区。在康定、泸定一带，出露有康定杂岩，其主要由片麻岩、混合岩等组成，是太古宙时期形成的结晶基底，经历了强烈的变质作用和构造变形，岩石结构复杂，矿物定向排列明显。在会理、会东等地，广泛出露中元古代的会理群，该群地层下部为优地槽建造，以火山岩、碎屑岩为主，反映了当时强烈的构造活动和火山喷发环境；上部为冒地槽建造，主要由浅变质的碎屑岩和碳酸盐岩组成，表明沉积环境逐渐趋于稳定。这些古老地层的存在，为研究四川地区早期岩石圈的形成和演化提供了关键线索，它们的岩石特征和构造变形历史，反映了岩石圈在漫长地质时期的演化过程。古生代地层在四川地区分布较为广泛，岩性丰富多样。寒武系地层在四川盆地边缘和川西南地区有出露，主要为海相沉积的碎屑岩和碳酸盐岩，含有丰富的三叶虫化石，这些化石不仅是地层划分和对比的重要依据，还反映了当时温暖浅海的沉积环境。奥陶系地层同样以海相沉积为主，岩性包括石灰岩、页岩等，在川南地区可见到厚层的石灰岩，其中发育有丰富的生物礁，表明当时的海洋生态环境适宜生物生长和礁体发育。志留系地层则以碎屑岩为主，常夹有页岩和泥岩，在川东北地区，志留系地层中含有较多的笔石化石，这些笔石的分布和演化特征，为研究该地区的古地理环境和生物演化提供了重要信息。泥盆系和石炭系地层在四川地区也有一定分布，主要为海相沉积的碳酸盐岩和碎屑岩，它们的沉积特征反映了当时的海平面变化和构造运动对沉积环境的影响。古生代地层的分布和岩性特征，反映了四川地区在古生代时期处于海洋环境，受到全球海平面变化和区域构造运动的共同影响，岩石圈在这种环境下经历了长期的沉积和改造过程。中生代地层在四川地区占据重要地位，是研究岩石圈演化的关键时期。三叠系地层在四川盆地广泛分布，下三叠统主要为海相沉积的石灰岩和碎屑岩，中三叠统逐渐过渡为海陆交互相沉积，含有大量的双壳类、腕足类等化石，反映了当时海陆环境的频繁交替。上三叠统则以陆相沉积为主，主要为砂岩、页岩和煤层，在川西地区，上三叠统的须家河组是重要的含煤地层，煤层的形成与当时温暖湿润的气候和茂密的植被有关，同时也表明该地区在晚三叠世时期已经逐渐由海洋环境转变为陆地环境。侏罗系地层主要为陆相沉积的砂岩、泥岩和页岩，在四川盆地内广泛分布，形成了独特的红层地貌，如川中地区的红色丘陵。这些红层的形成与当时炎热干燥的气候条件有关，红色的砂岩和泥岩中常含有石膏、岩盐等蒸发岩矿物，反映了干旱的沉积环境。白垩系地层同样以陆相沉积为主，岩性与侏罗系相似，但在沉积厚度和岩相变化上有所不同，白垩系地层在四川盆地东部地区沉积较厚，反映了当时该地区的构造沉降作用较强。中生代地层的演化过程，反映了四川地区岩石圈在这一时期经历了从海洋到陆地的转变，以及构造运动对沉积环境的深刻影响，岩石圈的结构和性质在这个过程中也发生了显著变化。新生代地层在四川地区主要为第四系地层，广泛分布于盆地内部和河谷地带。第四系地层主要由松散的沉积物组成，包括冲积物、洪积物、坡积物、残积物等，其厚度和岩性在不同地区差异较大。在成都平原，第四系地层主要为冲积物和洪积物，厚度可达数十米，由砂、砾石、黏土等组成，这些沉积物是在河流和洪水作用下形成的，反映了该地区在新生代时期的水动力条件和地貌演化。在山区，第四系地层则以坡积物和残积物为主，厚度相对较薄，主要由岩石风化产物和少量的坡积物组成，它们的分布和特征与山区的地形、气候和岩石性质密切相关。新生代地层的形成与现今的地质环境密切相关，它们记录了近期岩石圈的表面过程和环境变化，对于研究现代地质作用和岩石圈的表层结构具有重要意义。四川地区的岩石类型丰富多样，不同类型的岩石分布与地层和构造密切相关，对岩石圈结构产生着重要影响。岩浆岩在四川地区有广泛分布，主要包括侵入岩和喷出岩。在川西地区，分布着大量的花岗岩侵入体，如康定杂岩中的花岗岩，这些花岗岩是在深部岩浆侵入地壳的过程中形成的，其形成与板块碰撞和俯冲作用有关。花岗岩具有较高的硬度和强度，对岩石圈的力学性质产生重要影响，使得该地区的岩石圈在局部区域表现出较强的刚性。在晚二叠世时期，四川地区发生了大规模的火山喷发活动，形成了广泛分布的峨眉山玄武岩，这些玄武岩覆盖了四川盆地西部和南部的大片区域。玄武岩的喷发与地幔柱活动有关，大量玄武岩的喷出改变了岩石圈的物质组成和结构，使得该地区岩石圈的厚度和密度发生变化，对岩石圈的热结构和动力学过程产生了深远影响。沉积岩是四川地区分布最广泛的岩石类型之一，其分布与地层和沉积环境密切相关。在四川盆地，沉积岩主要包括砂岩、页岩、石灰岩等。砂岩主要分布于河流相和三角洲相沉积环境中，如侏罗系和白垩系地层中的砂岩，它们具有较高的孔隙度和渗透率，对地下水的储存和运移起着重要作用。页岩则多分布于湖泊相和海相沉积环境中，如三叠系和侏罗系地层中的页岩，页岩的富含有机质，是重要的烃源岩，同时页岩的力学性质相对较弱，在构造应力作用下容易发生变形和破裂。石灰岩主要分布于海相沉积环境中，如古生代地层中的石灰岩，石灰岩具有较高的抗压强度，但在地下水的溶蚀作用下，容易形成喀斯特地貌，改变岩石圈的表面形态和结构。沉积岩的分布和性质对岩石圈的结构和动力学过程产生着多方面的影响，它们的沉积过程记录了岩石圈的演化历史，而其力学性质和物质组成则影响着岩石圈的变形和稳定性。变质岩在四川地区主要分布于川西高原和盆地边缘的构造活动带。在龙门山断裂带和鲜水河断裂带等区域，出露有大量的变质岩，如片岩、板岩、千枚岩等。这些变质岩是在强烈的构造运动和高温高压作用下形成的，岩石中的矿物发生了重结晶和定向排列，形成了明显的片理构造。变质岩的形成与板块碰撞、俯冲和地壳深部的构造变形密切相关，它们的存在反映了岩石圈在构造活动带经历了强烈的改造和变形过程。变质岩的力学性质和结构特征与原岩和变质程度有关，一般来说，变质程度较高的变质岩具有较高的强度和稳定性，而变质程度较低的变质岩则相对较软，容易发生变形。变质岩的分布和特征对岩石圈的结构和动力学过程产生着重要影响，它们是岩石圈构造变形的产物，同时也影响着岩石圈在构造活动带的变形和演化。3.3地震活动特点四川地区地震活动频繁，其地震活动在空间分布、震级特征等方面表现出明显的特点，这些特点与该地区复杂的地质构造和岩石圈结构密切相关。从空间分布来看，四川地区的地震主要集中在几个特定的区域，这些区域往往与大型断裂带的分布高度吻合。龙门山断裂带是四川地区地震活动最为强烈的区域之一，2008年汶川8.0级特大地震和2013年芦山7.0级地震就发生在该断裂带上。龙门山断裂带位于青藏高原东缘与四川盆地的交界处，由于印度板块与欧亚板块的强烈碰撞，青藏高原物质向东挤出，在龙门山地区受到四川盆地刚性地块的阻挡，导致应力高度集中，从而引发了强烈的地震活动。鲜水河断裂带同样是地震频发的区域，该断裂带是一条左旋走滑断裂，历史上曾多次发生强烈地震，如1973年炉霍7.6级地震。鲜水河断裂带的地震活动与青藏高原内部的构造变形和物质运移密切相关，其走滑运动导致了断裂带两侧岩石的强烈变形和破裂，为地震的发生提供了条件。安宁河断裂带、小江断裂带等区域也时有地震发生，这些断裂带控制了川滇地区的构造格局和地震活动，它们的活动特征和力学性质对该地区的地震活动产生了重要影响。在这些主要断裂带之间，也存在一些中小规模的地震活动，它们的分布与区域内的次级断裂和构造变形有关，虽然震级相对较小，但同样对区域的地震活动格局产生影响。在震级特征方面，四川地区地震震级跨度较大，从微小地震到特大地震均有发生。其中，6级以上的中强地震时有发生，这些地震往往会造成较大的人员伤亡和财产损失。如前文所述的汶川8.0级特大地震，是新中国成立以来破坏性最强、波及范围最广的一次地震，其释放的能量巨大，对四川地区的地质构造和生态环境造成了深远影响。芦山7.0级地震、九寨沟7.0级地震等也都给当地带来了严重的灾害。除了这些强震，四川地区还存在大量的中小地震，这些地震震级一般在3-5级之间，虽然单个地震的影响相对较小，但它们的频繁发生反映了该地区地壳的活跃程度。这些中小地震往往是地壳内部应力调整的表现，它们的活动特征和分布规律对于研究区域地壳的稳定性和地震危险性具有重要意义。通过对地震震级-频度关系的分析，可以发现四川地区的地震活动符合古登堡-里克特关系，即震级越高，发生的频率越低；震级越低，发生的频率越高。这一关系为评估该地区的地震危险性和地震灾害风险提供了重要的理论依据。地震活动与岩石圈结构之间存在着紧密的内在联系。岩石圈结构的不均一性和地质构造的复杂性是导致地震活动的重要原因。在四川地区，不同构造单元的岩石圈厚度、物质组成和力学性质存在显著差异，这种差异使得岩石圈在受力时容易产生应力集中和变形不协调，从而引发地震。例如，在青藏高原东缘与四川盆地的交界处，由于岩石圈厚度的急剧变化和物质性质的差异，当受到板块碰撞的应力作用时，应力容易在边界区域集中，导致岩石的破裂和地震的发生。断裂带作为岩石圈中的薄弱部位，其活动性对地震活动起着关键作用。龙门山断裂带、鲜水河断裂带等大型断裂带，由于长期的构造活动，断裂带内岩石破碎，力学强度降低，容易在应力作用下发生错动和滑动，从而引发地震。岩石圈内部的低速层、高速体等异常结构也会影响地震波的传播路径和速度，进而影响地震的发生和传播特征。低速层的存在可能导致地震波能量的衰减和聚焦，增加地震的破坏力；高速体则可能改变地震波的传播方向，使得地震的影响范围和破坏程度发生变化。通过对地震活动与岩石圈结构关系的研究，可以更好地理解地震的孕育和发生机制，为地震预测和灾害防治提供科学依据。四、基于S波接收函数的四川地区岩石圈埋深研究4.1观测数据与台站分布本研究的数据主要来源于中国地震科学探测台阵、中国地震局地震台网以及部分区域地震台网的观测记录，这些数据涵盖了四川地区多个地震台站在不同时间段内记录的远震事件，时间跨度从[开始时间]至[结束时间]，共计获取了[X]个远震事件的地震波形数据。所使用的地震台站均配备了高精度的宽频带地震仪，能够准确地记录地震波的三分量信息，包括垂直分量（Z分量）和两个水平分量（径向分量R和切向分量T），仪器的频率响应范围为[具体频率范围]，能够有效捕捉到不同频率成分的地震波信号，为后续的S波接收函数计算和岩石圈结构研究提供了高质量的数据基础。四川地区的地震台站分布广泛，覆盖了四川盆地、川西高原、川滇边界等不同地质构造区域，共计[X]个台站。在四川盆地，台站分布相对较为密集，这是因为四川盆地是人口密集和经济发展的重要区域，对地震监测的需求更为迫切。这些台站均匀地分布在盆地内部及其周边地区，能够全面地监测盆地内的地震活动和岩石圈结构特征。例如，在成都平原及其周边地区，布置了多个台站，如成都台、德阳台、绵阳台等，这些台站可以有效地监测四川盆地内的地震波传播特征，为研究盆地内岩石圈的厚度变化和内部结构提供了丰富的数据。在川西高原，台站分布相对稀疏，但仍覆盖了主要的构造单元和断裂带，如鲜水河断裂带、龙门山断裂带等区域都有台站分布。这些台站对于监测川西高原地区的地震活动和岩石圈结构变化具有重要意义，能够捕捉到该地区由于板块碰撞和构造活动导致的岩石圈变形信息。在川滇边界地区，台站分布也较为关键，因为该地区是构造活动强烈的区域，地震频发。在小江断裂带、则木河断裂带等附近布置的台站，如西昌台、攀枝花台等，能够及时记录到该地区的地震事件，为研究川滇边界地区岩石圈的复杂结构和动力学过程提供了重要的数据支持。台站分布对研究结果的影响是多方面的。台站的空间分布密度直接影响研究结果的分辨率。在台站分布密集的区域，能够获取更多的地震波信息，对岩石圈结构的成像分辨率更高，可以更精确地揭示岩石圈内部界面的起伏和变化。在四川盆地内部，由于台站分布密集，通过对这些台站记录的地震波数据进行分析，可以清晰地成像出莫霍面和岩石圈-软流圈边界等重要界面的深度和形态变化，对岩石圈厚度的计算精度也更高。相反，在台站分布稀疏的区域，如川西高原的部分地区，由于地震波信息有限，成像分辨率相对较低，可能会导致对岩石圈结构的认识存在一定的不确定性。台站的分布位置与地质构造的关系也对研究结果产生影响。如果台站位于重要的构造单元或断裂带附近，能够更准确地捕捉到由于构造活动导致的岩石圈结构变化信息，为研究构造活动对岩石圈的影响提供关键数据。位于龙门山断裂带附近的台站，能够记录到地震波在断裂带附近的传播特征变化，从而推断出断裂带对岩石圈结构的影响，如岩石圈的变形、破裂以及物质运移等情况。台站分布的均匀性也会影响研究结果的可靠性。如果台站分布不均匀，可能会导致某些区域的地震波信息缺失或不足，从而影响对整个研究区域岩石圈结构的全面认识。因此，在研究过程中，需要综合考虑台站分布的各种因素，通过合理的数据处理和分析方法，尽可能地提高研究结果的准确性和可靠性。4.2岩石圈埋深的计算结果通过对四川地区多个地震台站的S波接收函数进行处理和分析，运用共转换点叠加（CCP）和波动方程偏移成像等技术，得到了该地区岩石圈埋深的计算结果。结果显示，四川地区岩石圈埋深呈现出显著的空间变化特征，不同构造单元的岩石圈厚度存在明显差异。在四川盆地，岩石圈厚度相对较大，平均厚度约为150-160km。从空间分布来看，盆地内部岩石圈厚度较为均匀，变化较小，但在盆地边缘地区，由于受到周边构造单元的影响，岩石圈厚度出现了一定程度的变化。在龙门山断裂带附近，岩石圈厚度有所增加，这可能是由于青藏高原物质向东挤出，在龙门山地区受到四川盆地刚性地块的阻挡，导致岩石圈发生变形和增厚。通过对该区域多个台站的S波接收函数成像分析，发现龙门山断裂带下方岩石圈内部存在明显的速度变化界面，这些界面的存在可能与岩石圈的变形和构造活动有关，进一步证实了岩石圈厚度增加的结论。在盆地南部，靠近川滇边界地区，岩石圈厚度则略有减小，这可能与川滇地区强烈的构造活动和岩石圈的改造有关。川西高原地区岩石圈厚度相对较薄，平均厚度约为80-100km。该地区岩石圈厚度的变化与构造活动密切相关，在构造活动强烈的区域，如鲜水河断裂带、龙门山断裂带等附近，岩石圈厚度变化较大，呈现出明显的不均匀性。在鲜水河断裂带下方，岩石圈厚度在短距离内出现了较大幅度的变化，这可能是由于断裂带的走滑运动导致岩石圈发生强烈的变形和破裂，使得岩石圈结构变得复杂，厚度出现异常变化。通过对该区域地震台站的S波接收函数分析，发现岩石圈内部存在多个速度异常区域，这些异常区域与断裂带的分布和活动特征密切相关，进一步表明了构造活动对岩石圈厚度的影响。在川西高原的一些相对稳定区域，岩石圈厚度相对较为稳定，但仍明显小于四川盆地的岩石圈厚度。川滇边界构造带岩石圈结构复杂，岩石圈厚度变化较大，呈现出不规则的分布特征。在小江断裂带、则木河断裂带等附近，岩石圈厚度变化尤为显著，部分区域岩石圈厚度可能小于80km，而在一些局部区域，由于构造作用的影响，岩石圈厚度可能会有所增加。小江断裂带作为川滇边界构造带的主要断裂之一，其活动强烈，对岩石圈结构产生了深远影响。通过对小江断裂带附近台站的S波接收函数研究，发现岩石圈内部存在复杂的速度结构和界面起伏，这表明岩石圈在断裂带的作用下发生了强烈的变形和改造，导致岩石圈厚度和结构的复杂性增加。川滇边界构造带还受到印度板块与欧亚板块碰撞以及太平洋板块向欧亚板块俯冲的共同影响，多种构造应力的叠加使得该地区岩石圈结构更加复杂，岩石圈厚度的变化也更加难以预测。为了更直观地展示四川地区岩石圈埋深的空间变化特征，制作了岩石圈厚度分布图（图1）。从图中可以清晰地看出，四川盆地岩石圈厚度较大，呈现出相对稳定的块状分布；川西高原岩石圈厚度较薄，且在断裂带附近存在明显的厚度变化梯度；川滇边界构造带岩石圈厚度变化复杂，呈现出不规则的分布形态。不同构造单元之间岩石圈厚度的差异，反映了各构造单元在地质历史时期经历的构造运动和演化过程的不同。四川盆地作为扬子板块的稳定组成部分，经历了长期的构造稳定期，岩石圈得以逐渐增厚和稳定；而川西高原和川滇边界构造带则受到强烈的构造活动影响，岩石圈发生了变形、破裂和改造，导致岩石圈厚度变薄且变化复杂。[此处插入岩石圈厚度分布图]四川地区岩石圈埋深的计算结果表明，该地区岩石圈结构复杂，不同构造单元的岩石圈厚度存在显著差异，这种差异与区域地质构造背景和构造活动密切相关。通过对岩石圈埋深的研究，为深入理解四川地区的岩石圈动力学机制提供了重要的基础数据和约束条件。4.3结果分析与讨论四川地区岩石圈埋深呈现出显著的空间变化特征，不同构造单元岩石圈厚度的差异与该地区复杂的地质背景密切相关。四川盆地作为扬子板块西缘的稳定地块，其岩石圈厚度较大且相对稳定，这主要归因于其古老的刚性基底和长期的构造稳定。盆地基底由前震旦纪变质岩系组成，经历了多期构造运动的改造，形成了相对稳定的构造格局，这种稳定的基底结构为岩石圈的增厚提供了基础。在漫长的地质历史时期中，四川盆地经历了海相沉积、陆相沉积等多个阶段，沉积了巨厚的地层，这些地层的堆积进一步增加了岩石圈的厚度。此外，四川盆地位于板块内部相对稳定的区域，受到的构造应力相对较小，岩石圈在长期的地质演化过程中得以保持相对稳定的状态，没有受到强烈的构造变形和破坏，从而使得岩石圈厚度能够维持在较大的水平。川西高原岩石圈厚度较薄，这与该地区强烈的构造活动密切相关。川西高原处于青藏高原东缘，是印度板块与欧亚板块碰撞的前沿地带，受到板块碰撞的强烈挤压作用，岩石圈发生了强烈的变形和改造。印度板块持续向北推挤欧亚板块，导致青藏高原物质向东挤出，川西高原首当其冲，岩石圈在这种强烈的构造应力作用下发生破裂、变形和物质迁移，使得岩石圈厚度减薄。该地区广泛分布的断裂带，如鲜水河断裂带、龙门山断裂带等，这些断裂带的活动导致岩石圈内部的物质重新分布和结构调整，进一步加剧了岩石圈的减薄。断裂带的走滑运动和逆冲运动使得岩石圈发生破裂和错动，部分岩石圈物质被挤出或发生变形，从而导致岩石圈厚度的减小。川滇边界构造带岩石圈结构复杂，厚度变化较大，这是多种构造应力叠加的结果。该构造带不仅受到印度板块与欧亚板块碰撞的影响，还受到太平洋板块向欧亚板块俯冲的远程效应影响，多种构造应力在该区域相互作用，导致岩石圈发生复杂的变形和调整。小江断裂带、则木河断裂带等主要断裂的活动，控制了该地区岩石圈的结构和厚度变化。这些断裂带的活动性质多样，既有走滑运动，又有逆冲运动，它们的相互作用使得岩石圈内部的应力分布极为复杂，岩石圈在不同部位受到不同方向和大小的应力作用，从而导致岩石圈厚度的不规则变化。川滇边界构造带还存在强烈的岩浆活动和热液活动，这些活动对岩石圈的物质组成和结构产生了重要影响，进一步增加了岩石圈结构的复杂性和厚度变化的多样性。与其他相关研究成果相比，本研究得到的四川地区岩石圈埋深结果在整体趋势上具有一定的一致性，但在局部细节上仍存在一些差异。张耀阳等人利用S波接收函数波动方程叠后偏移方法研究华南大陆岩石圈结构时发现，150km以上的厚岩石圈主要位于四川盆地，不足100km的薄岩石圈主要分布于川东褶皱带和华夏块体，这与本研究中四川盆地岩石圈厚度较大、川西高原和川滇边界部分地区岩石圈厚度较薄的结果相符。在一些局部区域，不同研究的结果存在差异。部分研究中四川盆地岩石圈厚度的具体数值与本研究略有不同，这可能是由于研究方法、数据处理过程以及台站分布的差异所导致。不同的研究方法在数据处理和成像过程中可能会引入不同的误差和不确定性，从而影响对岩石圈厚度的计算结果。台站分布的差异也会对研究结果产生影响，如果台站分布不均匀或在某些关键区域缺失，可能会导致对岩石圈结构的采样不全面，从而影响对岩石圈厚度的准确估计。这些差异的产生原因主要包括以下几个方面。不同研究采用的数据处理方法和成像技术存在差异，每种方法都有其自身的优缺点和适用范围，对数据的处理和解释方式也不尽相同，这可能导致在相同的地质条件下得到不同的岩石圈厚度结果。在S波接收函数的计算过程中，不同的反褶积算法、时窗选取方法以及叠加技术等都可能对结果产生影响；在成像过程中，不同的成像算法对地震波传播路径的模拟和对界面反射信息的提取方式也存在差异，这些差异都可能导致成像结果的不同。研究中所使用的数据质量和台站分布情况也会对结果产生重要影响。如果数据存在噪声干扰、仪器故障或数据缺失等问题，可能会影响S波接收函数的计算精度和成像质量，从而导致对岩石圈厚度的估计出现偏差。台站分布的不均匀性会导致对某些区域的岩石圈结构采样不足，使得研究结果不能准确反映该区域的真实情况。如果在川西高原等构造复杂区域台站分布稀疏，可能会遗漏一些岩石圈结构的细节信息，从而导致对该区域岩石圈厚度的估计不准确。地质构造的复杂性和岩石圈结构的横向变化也是导致研究结果差异的重要原因。四川地区地质构造复杂多样，岩石圈结构在横向和纵向上都存在显著的变化，不同研究在数据采集和分析过程中可能对这些变化的敏感性不同，从而导致结果的差异。在一些构造活动强烈的区域，岩石圈内部可能存在复杂的速度变化和界面起伏，不同研究对这些复杂结构的成像能力和解释方式不同，也会导致对岩石圈厚度的认识存在差异。五、四川地区岩石圈动力学机制探讨5.1岩石圈变形机制四川地区岩石圈变形特征显著，呈现出复杂多样的形态。在不同构造单元，岩石圈变形表现出明显的差异。在四川盆地，由于其处于相对稳定的构造环境，岩石圈变形以较为和缓的褶皱变形为主，地层呈现出平缓的褶皱形态，褶皱的幅度和波长相对较大，反映了岩石圈在长期的地质演化过程中受到的构造应力相对较弱且较为均匀。在川西高原和川滇边界构造带，岩石圈变形则表现得极为强烈，以断裂变形和强烈的褶皱变形为主要特征。这些地区分布着众多大型断裂带，如龙门山断裂带、鲜水河断裂带、小江断裂带等，断裂带两侧的岩石发生了显著的错动和位移，导致岩石圈结构的强烈破坏和变形。在龙门山断裂带，2008年汶川地震使得断裂带两侧的山体发生了巨大的错动，垂直位移可达数米，水平位移也十分显著，这种强烈的断裂变形对岩石圈的结构和稳定性产生了深远影响。这些地区的褶皱变形也十分复杂，褶皱的轴面往往发生倾斜和倒转，褶皱的紧闭程度较高，反映了岩石圈在强烈的构造应力作用下发生了复杂的变形过程。岩石圈变形机制主要包括弹性变形和塑性变形两种基本类型，它们在四川地区岩石圈变形过程中都起着重要作用，且相互影响、相互制约。弹性变形是指岩石在受力时发生变形，当外力消失后，岩石能够恢复到原来的形状和大小的变形过程。在四川地区，当岩石圈受到较小的构造应力作用时，会发生弹性变形。在地震孕育过程中，岩石圈在构造应力的缓慢积累下，会发生弹性变形，储存应变能。随着应力的不断增加，岩石的弹性变形逐渐增大，当应力达到岩石的弹性极限时，岩石将发生塑性变形或破裂。弹性变形具有可逆性，它在岩石圈变形的早期阶段起到了重要的作用，能够调节岩石圈内部的应力分布，延缓岩石圈的破裂和变形。塑性变形是指岩石在受力超过其屈服强度后，发生不可逆的永久变形的过程。在四川地区，由于受到印度板块与欧亚板块碰撞的强烈影响，岩石圈受到巨大的构造应力作用，塑性变形广泛存在。在川西高原和川滇边界构造带，岩石在高温高压和强烈的构造应力作用下，矿物晶体发生位错滑移、晶界滑动等塑性变形机制，导致岩石的形状和结构发生永久性改变。在龙门山断裂带附近的岩石中，通过显微镜观察可以发现矿物晶体的位错滑移现象，这是塑性变形的重要证据。塑性变形使得岩石圈的结构变得更加复杂，岩石的力学性质也发生了改变，从而影响了岩石圈的变形和演化过程。塑性变形还与岩石的流变性质密切相关，不同类型的岩石具有不同的流变特性，这决定了它们在塑性变形过程中的行为和方式。在高温高压条件下，岩石的流变性质会发生变化，使得塑性变形更加容易发生。在实际的岩石圈变形过程中，弹性变形和塑性变形并非孤立存在，而是相互作用、相互转化的。在构造应力作用下，岩石首先发生弹性变形，随着应力的逐渐增加，当达到岩石的屈服强度时，岩石开始发生塑性变形。在塑性变形过程中，岩石内部的应力得到一定程度的释放和调整，可能会导致部分塑性变形转化为弹性变形。在地震发生时，岩石圈中的应力突然释放，原本处于塑性变形状态的岩石可能会发生弹性回跳，产生地震波，这就是地震的弹性回跳理论的基本原理。岩石圈变形还受到多种因素的影响，如岩石的类型、温度、压力、应变速率等。不同类型的岩石具有不同的力学性质和变形行为，花岗岩等刚性岩石在变形过程中往往表现出较强的脆性，容易发生破裂和断裂；而页岩等柔性岩石则具有较好的塑性，更容易发生塑性变形。温度和压力对岩石的变形机制也有重要影响，在高温高压条件下，岩石的塑性变形能力增强，弹性变形相对减弱；而在低温低压条件下，岩石则更容易表现出弹性变形特征。应变速率的大小也会影响岩石的变形机制，当应变速率较快时，岩石更容易发生脆性破裂；当应变速率较慢时，岩石则有更多的时间发生塑性变形。5.2深部热结构与动力学四川地区深部热结构呈现出复杂的特征，这与该地区独特的地质构造背景和岩石圈演化历史密切相关。深部热结构主要受到地幔热流、岩石圈热导率以及地壳内部放射性元素衰变等因素的综合影响。地幔热流作为地球内部热量的主要来源之一，对深部热结构起着基础性的控制作用；岩石圈热导率决定了热量在岩石圈中的传输效率和分布特征；地壳内部放射性元素衰变产生的热量也在一定程度上影响着深部热结构的局部变化。通过对四川地区地温梯度的研究发现，不同构造单元的地温梯度存在显著差异。在四川盆地，地温梯度相对较低，平均地温梯度约为25-30℃/km，这主要是由于四川盆地岩石圈厚度较大，热流向上传输受到一定阻碍，导致地温梯度相对较小。盆地内沉积层较厚，且岩石的热导率相对较低，也进一步减缓了热量的传输，使得地温梯度维持在较低水平。在川西高原，地温梯度相对较高，平均地温梯度可达35-45℃/km，这与该地区岩石圈较薄、构造活动强烈密切相关。川西高原处于青藏高原东缘，受到印度板块与欧亚板块碰撞的强烈影响，岩石圈发生强烈变形和减薄，地幔热流更容易向上传输，导致地温梯度升高。该地区广泛分布的断裂带为热液活动提供了通道，使得深部热量能够更有效地向上传导，进一步增加了地温梯度。川滇边界构造带地温梯度变化较为复杂，在一些构造活动强烈的区域，地温梯度较高，而在相对稳定的区域，地温梯度则相对较低。这是因为川滇边界构造带受到多种构造应力的叠加影响，岩石圈结构复杂，深部热流的传输受到多种因素的制约，导致地温梯度呈现出复杂的变化特征。深部热结构对岩石圈动力学过程有着重要影响，热对流和热传导是深部热结构影响岩石圈动力学的主要方式。热对流是地幔物质在温度梯度作用下发生的大规模循环运动，它在岩石圈动力学中扮演着关键角色。在四川地区，热对流对岩石圈的运动和变形产生了重要影响。在川西高原和川滇边界构造带，由于地幔热对流的作用，地幔物质上升，对岩石圈底部产生向上的作用力，导致岩石圈发生隆升和伸展变形。这种热对流驱动的岩石圈变形与区域内的构造运动相互作用，进一步加剧了岩石圈的复杂性。热对流还会导致岩石圈内部物质的重新分布和结构调整，对岩石圈的物质组成和力学性质产生影响。热传导是热量通过物质分子的热运动从高温区域向低温区域传递的过程，它在岩石圈的热量传输和温度分布中起着重要作用。在四川地区，热传导对岩石圈的稳定性和变形也产生了影响。在四川盆地，由于岩石圈较厚，热传导相对缓慢，使得岩石圈内部温度分布相对均匀，岩石圈的稳定性较高。而在川西高原和川滇边界构造带，由于岩石圈较薄且构造活动强烈，热传导速度相对较快，岩石圈内部温度变化较大，导致岩石圈的力学性质发生改变，容易发生变形和破裂。热传导还会影响岩石圈内部的化学反应和矿物相变，进而影响岩石圈的物质组成和结构。在高温高压条件下，岩石中的矿物可能会发生相变，导致岩石的物理性质发生变化，这对岩石圈的动力学过程产生了重要影响。深部热结构与岩石圈动力学之间存在着密切的耦合关系，它们相互作用、相互影响，共同塑造了四川地区现今的地质构造格局。深部热结构的变化会导致岩石圈动力学过程的改变，而岩石圈动力学过程的变化也会反过来影响深部热结构。在印度板块与欧亚板块碰撞的过程中，深部热结构发生了显著变化，地幔热流重新分布，这导致了岩石圈的变形和运动，形成了川西高原和川滇边界构造带的复杂构造格局。而岩石圈的变形和运动又会改变深部热流的传输路径和分布特征，进一步影响深部热结构。这种耦合关系在地震活动、岩浆活动等地质现象中也有明显体现。地震活动往往发生在岩石圈的薄弱部位，而这些薄弱部位的形成与深部热结构和岩石圈动力学的相互作用密切相关。岩浆活动则是深部热结构和岩石圈动力学共同作用的结果，地幔热对流和岩石圈的变形为岩浆的产生和上升提供了条件。5.3板块相互作用与岩石圈演化四川地区位于印度板块与欧亚板块碰撞带的东侧，其岩石圈演化受到板块相互作用的深刻影响。印度板块持续向北推挤欧亚板块，这种强烈的碰撞作用是四川地区岩石圈演化的主要驱动力之一。在碰撞过程中，印度板块的物质向北楔入欧亚板块之下，导致青藏高原不断隆升，并向东挤出物质，四川地区首当其冲，成为了构造应力的强烈作用区域。这种板块间的相互作用不仅改变了岩石圈的物质组成和结构，还对岩石圈的厚度和变形特征产生了深远影响。印度板块与欧亚板块的碰撞对四川