青藏高原及其中国地区大陆地壳岩石层变形演化动力学数值模拟.pdf

摘 要 青藏高原作为全球构造活动最剧烈的地区之一，其独特的构造和地球物理环 境使其成了大陆动力学研究的天然实验室。其地壳岩石层变形运动、应力一应变 场、演化和隆升过程以及动力学机制一直是地球科学家们研究的热点。本文利用 固体力学理论和有限元方法分别进行了GPs观测速度约束层边界F的青藏高原 及邻近区域的二维稳态运动学，和印度板块挤压作用下近20万年中国大陆地壳 岩石层变形的三维动力学演化数值模拟。 在_维稳态运动学数值模拟中，本文以GPs观测结果为边界约束，利用有 限元分析方法模拟了青藏高原及邻区的地壳运动速度场、应力．应变场。为了对 比不同的地壳物质性质和断层接触对青臧高原变形、应力一应变场分布的影响， 模型按照青藏高原及临近区域的材料性质和接触方式分成了四种不同的情况： 元的网格划分做了改进和创新。对比四种模拟结果之问以及它们和观测数据的异 同显示：】．在GPs观测的时间尺度内，DP塑性材料和弹性材料性质对模拟青 藏高原变形的影响没有大的差异；2．可滑移的断层的存在对高原地壳岩石层形 变运动和应力．应变场影响很大，特别是增加了青藏高原的东西向拉张趋势和侧 向滑移的速度；3．在短时间尺度内，弹性非连续接触模型(M．ET)预测结果最 接近观测情况；4．模拟结果与GPS观测吻合程度较高，计算结果较为可靠。本 文在M．ET模型基础上对阿尔金和羌塘附近局部区域利用三个GPs观测点数据 对模型进一步约束，从而得到接近青藏高原真实情况的形变运动场和构造应力 场。数值模拟的结果表明，青藏高原现今地表应力场和深部应力场在很大程度上 是同源的：青藏高原区域中部和南部拉张应力和挤压应力同时存在，且拉张应力 4土1) 作用受断层和其他作用力影响很人：阿尔金断裂带当前滑移速率较低(7 mm／a，现今处于活动性较差的状态；现今的青藏高原变形中，起主导作用的是 地壳的分散式缩短和增厚。 为了研究青减高原及巾国大陆地壳岩石层近期的动力学演化过程，本文以印 度板块向北挤压作为大陆变形运动的主要驱动力。考虑岩石层(地壳和上地幔顶 部)的情况，将中国大陆地壳按照地质构造性质划分为10个不同的构造块体， 利用ANsYs有限元软件建立了青藏高原及中国大陆地壳岩石层变形的三维动力 学有限元模型。模型假定各块体之间运动连续，并将青藏高原及周边区域视为幂 指数蠕变性材料。 IIIM 论文对四种不同边界条件模型(MI，MII，M IV)进行了数值模拟。 模拟结果显示，东部边界和南部边界的红河断裂带区域为自由边界的模型MI 结果比较合理，同时蠕变幂指数(岩石力学参数)n二3时结果与实际情况比较吻 合。在此基础上，计算了青藏高原及邻区高原的隆升速率、中国大陆地壳岩石层 形变速率以及应力场，并将计算结果与实际观测资料(如GPS)和其他模拟工作 进行对比。结果表明，本文模型预测的结果与观测数据有较好的吻合。 三维动力学数值模拟显示，1．印度板块挤压是青藏高原及中国大陆地壳岩 石层变形演化的最主要的驱动力，决定了青藏高原的变形运动格局和演化过程。 同时也不同忽视菲律宾板块和太平洋板块的影响。2．在印度板块的挤压作用下， 现今青藏高原地壳变形以南北向缩短和增厚为主，红河断裂带东西向自由滑移运 动对青藏高原地壳水平运动影响很大。3．近20万年内青藏高原东西方向上的隆 升速度变化较小，但在南北方向上隆升速度存在明显的差异。4．中国大陆地壳 岩石层垂直变形速度和应力场分布有明显的区域性特征。应力和垂直形变速率分 布为西高东低。5．近20万年中国大陆岩石层形变运动和应力场格局较为稳定。 另外，模拟结果与观测结果之间的存在着一’定的差异，暗示着重力势、断层作用、 上地幔小尺度对流以及剥蚀等因素在青藏高原及中国大陆地壳岩石层的变形和 演化过程中可能起着重要的作用。 关键词： 青藏高原、GPs、速度场、应力场、隆升、演化 Abstract Asits tectonicstructureand P1ateau s口ecial geophysicalenviro砌ent，Tibet 口ro