青藏高原东缘构造应力场转换的动力学机制及其数值模拟研究

青藏高原东缘处于印度 - 欧亚板块碰撞的关键构造区域。该区域构造应力场演化规律和动力学机制一直是地球动力学领域研究热点。这里是高原物质向东侧逸散主要通道，也是强震活动集中发生区域。构造应力场时空转换特征直接影响地壳变形模式和地震危险性评估。

构造应力场是地壳内部因构造活动形成的应力状态，其转换过程体现板块运动边界条件、深部动力学过程和浅部响应之间复杂关联。从动力学机制看，青藏高原东缘应力场转换受三方面因素影响：印度板块持续向北挤压、高原隆升引发的重力势能差以及东亚块体阻挡效应。这些动力因素相互作用，让区域应力场呈现明显分区性和阶段性特点。

开展数值模拟研究是探究这一转换机制有效手段。关键是建立符合实际地质结构的力学模型，运用有限元或离散元方法模拟不同边界条件下应力场时空变化过程。开展这类研究有系统操作步骤，要依据地质和地球物理资料构建三维地质模型，确定合理岩石力学参数和边界条件，通过数值计算模拟应力场动态变化，结合震源机制解、地应力测量等野外观测数据验证和校正模拟结果。

这项研究实际应用价值显著。能加深对陆 - 陆碰撞带构造变形机制理论认知，为区域地震危险性预测和资源勘探提供动力学依据。在工程抗震设计和地质灾害防治领域，准确掌握构造应力场转换规律对评估场地稳定性、制定防灾策略有直接指导作用。

开展青藏高原东缘构造应力场转换的数值模拟研究，能完善大陆动力学理论体系，服务国家重大战略需求，有重要科学价值和现实意义。

青藏高原东缘是印度 - 欧亚板块碰撞效应朝东传递的关键区域，该区域构造应力场特征及其转换现象是认识区域动力学演化的重要窗口。研究区域涵盖松潘 - 甘孜地块、扬子地块西缘、川滇菱形块体，还有龙门山、鲜水河断裂带等核心构造单元。自晚新生代开始，这里历经了多期较为明显的构造应力场转换。

整合野外断层擦痕与节理产状统计、地震机制解、GPS位移场、古地磁以及同位素测年等多源数据，能够系统地还原不同时期的应力状态。松潘 - 甘孜地块早期应力场主要是近南北向主压应力，呈现出逆冲 - 走滑应力机制，到了晚期应力方向渐渐向东偏转，形成北东向挤压与南东向伸展的共轭系统。扬子地块西缘的应力场从早期的近东西向挤压，慢慢转变为晚期的北西 - 南东向拉张。川滇菱形块体长时间受到南东向逃逸构造的控制，主压应力方向从早期的北东向，逐渐变成晚期的近东西向。

应力场转换的时间节点可以通过地层不整合面、古应力遗迹测年、GPS速率突变等证据来确定，大致是以中更新世作为分界。早期（大约3.5 - 0.8Ma）主要受到印度板块直接挤压的控制，而晚期（大约0.8Ma到现在）则是由高原物质东逸和深部地幔流动共同驱动。在空间分布方面，龙门山断裂带呈现出应力方向的梯度变化特征，鲜水河断裂带记录了走滑分量增强的转换特点。

此次转换的核心特征是主压应力轴从近南北向朝着北东东向进行系统性的偏转，这一现象反映出印度板块碰撞从前期的刚性挤入到后期的侧向扩展的演化过程，为区域地震危险性评价和构造建模提供了关键的动力学边界条件。

构造应力场转换的动力学模型构建是数值模拟研究的核心部分。这个过程本质上就是建立一个符合地质实际情况的数学物理框架，从而对区域构造变形的驱动机制以及响应过程进行定量描述。

在构建模型的时候，第一步要做的是依据研究区的地质构造背景来明确几何边界条件。就拿青藏高原东缘的模型范围来说，需要把川滇块体、巴颜喀拉块体等主要活动块体都涵盖进去，而边界的确定要结合GPS观测速率场以及主要断裂带的分布情况来进行。断层的三维形态需要参考深部地球物理探测结果，要具体刻画它的倾角、走向还有延伸深度。就像鲜水河断裂带的展布特征，需要进行精确的建模，这样才能够反映出它对区域应力传递所起到的控制作用。

准确识别并且量化驱动因素是保证模型合理的关键。研究区应力场的主要动力来源是印度板块对欧亚板块持续的远程推挤，这种作用力的大小和方向能够通过板块运动学模型反演得到。与此同时深部地壳流的上涌效应、块体之间的相互作用，比如龙门山断裂带东西两侧的力学差异，也都需要考虑进去。这些因素一起构成了模型的边界荷载和内部体力条件。在概念模型当中，这些驱动因素会通过应力传递路径和断层活动性关联起来。例如远程推挤应力在遇到高强度块体的时候会发生偏转，进而使得局部应力场方向发生转换。

当把概念模型转化成数值力学模型时，需要选择合适的本构关系和破坏准则。岩石的弹塑性本构能够比较好地描述地壳材料在长期构造应力作用下的变形行为，它的屈服准则通常采用Mohr - Coulomb模型，其表达式为：$\tau = c+\sigma$n\tan\phi ，这里面的$\tau$代表的是剪应力，$c$指的是黏聚力，$\sigma$n是正应力，$\phi$表示的是内摩擦角。断层的力学行为是通过摩擦定律来表征的，比如速率 - 状态摩擦定律就可以模拟断层滑动过程中的弱化现象。关键参数的确定需要综合岩石力学实验结果、地质和地球物理观测数据，像地壳介质密度与波速的对应关系、断层闭锁段的摩擦系数等。模型是不是合理最终要通过和实际观测进行对比来验证，这里面包括地表位移场、地震活动分布以及应力方向测量结果，通过这样的验证来确保数值模拟能够真实地再现研究区构造应力场转换的动力学过程。

数值模拟是探究构造应力场转换动力学机制的重要技术手段。数值模拟的核心在于搭建数学物理模型，把复杂地质过程转化为能够计算的数值问题。本次研究选用FLAC3D软件进行模拟，这款软件是基于拉格朗日有限差分法开发的，它非常适合模拟大变形的地质过程，能够很好地处理材料和几何的非线性问题，与青藏高原东缘构造变形强烈的特点相匹配。

模拟方法的物理依据是岩石力学的本构关系，其核心控制方程是应力平衡方程：

这里面，\(\sigma_{ij}\)代表应力张量，\(\rho\)指的是岩石密度，\(g_i\)为重力加速度分量。在进行计算的时候采用显式差分格式，通过虚拟网格节点的运动方程一步步去求解系统反应。
几何模型严格按照实际地质构造的特点来构建。该模型覆盖了青藏高原东缘的关键构造单元，其水平方向的尺寸是800公里乘600公里，深度一直延伸到岩石圈底部，也就是120公里的位置。网格划分采用非均匀加密的方式，在断层带附近，网格大小设定为1到2公里，而在离构造区比较远的地方，网格放宽到5到8公里，这样做既能够保证计算精度，又可以提升运算效率。像龙门山断裂、鲜水河断裂这些主要的断裂带，都按照实际产状被放进模型里，断层的倾角是根据野外地质测量以及地震勘探所得到的数据来确定的。
力学参数的选择结合了室内岩石力学实验和区域地质研究的结果。上地壳的花岗岩密度设定为2700千克每立方米，弹性模量为60吉帕，泊松比是0.25；下地壳的镁铁质岩密度调整到2900千克每立方米，弹性模量增加到80吉帕。断层带参数是通过敏感性分析来确定的，黏聚力处于1到5兆帕之间，内摩擦角在25度到35度之间，这些数值参考了汶川地震科学钻探测到的数据。模型采用了莫尔 - 库仑破坏准则，该准则能够比较好地模拟岩石的脆性破裂情况。
边界条件的设置是为了还原印度板块挤压欧亚板块所产生的效果。在模型底部施加了完全固定的约束，以此来模拟岩石圈底部的刚性支撑。在东西两侧边界施加了水平方向的位移约束，在南部边界施加了向北的匀速位移载荷，其速率为5毫米每年，这和GPS测量到的现在地壳运动速率是一样的。北部边界被设定为自由边界，这样可以让模型在受到挤压的时候能够自然变形。

初始应力场的设定是根据区域地震矩张量解和水压致裂测量的数据。通过施加和实测值一致的差应力，使得模型在模拟开始之前就处于接近真实地质情况的应力环境。为了验证参数的可靠性，研究开展了系统的敏感性分析，也就是逐个改变断层摩擦系数（±20%）、边界加载速率（±30%）、岩石流变参数（±25%），查看这些变化对应力场分布产生多大的影响，最后确定最佳的参数组合。这种参数校验的方法能够保证模拟结果符合地质实际情况，为后续的分析奠定可靠的数值基础。

数值模拟结果清楚地呈现出青藏高原东缘构造应力场在时间和空间方面的演变特点。通过数值模拟的办法，得到了不同地质历史时期应力分布情况、应变积累状态、断层位移量以及块体运动速度场的数据。

可视化的结果显示早期阶段最大主应力方向是北东东向，应力集中区主要出现在鲜水河断裂带和龙门山断裂带交汇的位置，此地应变能密度较高且应变张量场呈现出明显的剪切特征。随着模拟时间不断推进，应力方向逐渐顺时针转动，最后变成南东东向，应力集中区也朝着东南方向移动，这情况和GPS观测到的现今地壳运动速率场分布情况十分相符。计算得出的断层位移数据显示鲜水河断裂带主要是左旋走滑运动，而龙门山断裂带呈现的是逆冲分量和右旋走滑的复合运动特征。块体运动速度矢量也显示出川滇块体相对于华南块体有明显往东南方向逃逸的趋势。

对模拟结果和实际地质观测数据进行对比验证，发现两者既存在高度的一致性，又有一些内在的差异。模拟得到的震源机制解，特别是P轴方位角，和区域内历史地震震源机制解统计结果在整体分布格局上大致是一样的，都能够反映出区域构造应力场的主压应力方向。模拟的块体运动速度场与GPS连续站观测的地表位移速率在方向和大小方面匹配得比较好，这说明数值模型可以有效地捕捉宏观地壳变形的行为。不过在局部区域像龙门山断裂带中南段，模拟得到的逆冲分量比地质剖面显示的长期构造缩短量稍微小一些。这种差异有可能是因为模型把深部岩石圈的流变学参数进行了简化，或者是没有充分考虑到古构造残余应力所产生的影响。

对转换过程的动力学机制进行分析可以发现，区域构造应力场的转换是多种因素一起作用而产生的结果。印度板块持续不断地向北推挤欧亚板块，这是区域变形的根本动力来源。当推挤方向稍微往东边偏一点的时候，会直接使得东缘地区主压应力方向出现系统性的偏转。下地壳和岩石圈地幔的物质流动也就是“地壳流”的侧向扩展和调整，给上地壳的块体运动提供了重要的解耦边界条件，让应力能够高效率地传递到远离碰撞带的前陆地区。块体边界的几何形态以及断层本身的力学性质是重要的约束条件。例如龙门山断裂带产状陡立，对川滇块体向东的运移起到了阻碍作用，使得应力在这个地方高度集中并且转化为垂向隆起。

分析应力传递路径可知，应力主要沿着大型断裂带高效地传导。来自青藏高原内部的应力经由昆仑断裂和鲜水河断裂传递到川滇块体，进而让龙门山断裂带被激活。不同区域的转换表现存在差异，根本原因在于岩石圈介质不均匀。就像四川盆地坚硬的克拉通基底形成了强大的阻挡力量，而西侧的松潘 - 甘孜褶皱带相对比较软弱，更容易发生变形。这个动力学模型明确关联上了区域构造演化过程，也就是从早期以碰撞挤压为主的情况，转变到后期以块体侧向挤出和深部地壳流调整为主的过程，这为解释青藏高原东缘复杂构造地貌的形成提供了定量化的力学依据。其动力学平衡方程是可以表达成这样的形式：

其中$\sigma$代表的是应力张量，$\rho$表示的是密度，$g$指的是重力加速度，这个方程是模拟的基础条件。

这项研究对青藏高原东缘构造应力场转换的动力学机制和数值模拟进行了系统分析，然后得出了下面这些结论。

青藏高原东缘处于印度板块与欧亚板块碰撞挤压的关键位置，该区域构造应力场转换一方面体现出区域地壳变形复杂过程，另一方面也体现出深部动力作用和浅部地质响应存在相互关联。研究发现，这片区域应力场转换主要受板块边界力、深部地幔对流以及岩石圈强度空间差异影响，并且东向挤出和南东向旋转叠加效应在此表现得十分明显。

在数值模拟方面，搭建三维有限元模型，把实际地质和地球物理数据结合起来，成功将应力场从近东西向向南东向的动态转换过程进行了还原，同时也对边界条件和介质属性对应力分布的调控作用进行了验证。模拟结果显示应力场转换不是简单线性变化，而是呈现出阶段性特点，这种现象和区域断裂系统活动性存在紧密联系。

在实际应用当中，研究成果能为认识青藏高原东缘地震活动孕育机制提供重要依据，尤其是在强震危险区判定上能起到指导作用。在对比不同模型参数敏感性之后，明确了岩石圈流变结构和断裂带摩擦系数对应力场转换的影响程度，这为区域动力学模型优化提供了量化约束。

总体来说，这项研究加深了对青藏高原东缘构造演化过程的理解，还为地质灾害防治以及资源勘探提供了理论支撑，体现出构造应力场研究在地球动力学领域具有实际应用价值。