青藏高原下地壳物质流动充分与必要条件数值模拟研究.pdfVIP

·270· 中国地球物理2011 (8)地球内部结构及其动力学 青藏高原下地壳物质流动的充分及必要条件的数值模拟研究 onthe and conditions Numericalsimulation sufficient necessary crustalchannelflowunder Plateau aboutthelower Tibetan 杨辉+滕吉文王谦身 皮娇龙 Hui Jiwen Pi WangQianshen Yang Teng Jiaolong 中国科学院地质与地球物理研究所 100029北京 flow)的概念以后，近年来诸多的有关 自20世纪80年代Morgan(1987)等人提出通道流(Channel 青藏高原下地壳通道流的文章不断发表，在他们的模型中主要是设定上地壳和上地幔的物质结构近似于刚 性．而之间的中地壳或下地壳为强度明显低的塑性流变层，且可以流动。这个模型在解释一些地区的构造 应力及变形时得到了很好的应用。但是，在一些工作中，根据地壳流的空间分布和深部探测资料(如电阻 and 率，Rippe 界处，特别是Moho界面处的流动速率最大。但是青藏高原地壳的底界反射波能量很强，属于高速的物性， 因而这种现象很难从介质属性和层圈分层机理上给予合适的解释。 et et 已开展的相关数值模拟工作(Beaumontal，2004；Roydenal，1997)中，均将下地壳设置为黏度系 数很低的物质组构，则下地壳便可以产生通道流的效应。同样，通道流的存在必须是低黏度物质被高黏度 的物质所包围，这样，在通道两侧，即上、下边界处就必须存在剪切力系的作用。 然而，下地壳这种特异黏滞系数的设定通常并非是针对具体地域的真实地壳介质属性测定而得，这便 牵涉到下地壳是软弱物质的真实依据所在。分层结构的地球内部各圈层介质力学属性的黏滞系数可以通过 理论计算，也可以通过高温高压实验、重力学、大地电磁、地震学和大地测量等多种途径获得。可是各种 方法针对同一地区所给出的结果差异的误差范围在通常情况下可以达到l～2个数量级，甚至达到更高的 数量级。在误差范围内的取值便能造就数值模拟中通道流现象的产生，由此可见，这样的数值模拟结果可 以作为讨论问题时的参考，但若作为定量的模型则缺乏足够的说服依据。而青藏高原的下地壳和地幔盖层 均为地震波速度随深度增加而递增的高速固态物质，这种边界限定了下地壳不易任意流动。若下地壳要相 对于上地壳和上地幔盖层有明显的运动，则不但需要有足够大的边界力施加于下地壳介质上，而且在其物 质运移的过程中上、下界面会产生一系列的物理一力学响应，如摩擦产生的热效应、介质同源同深的破裂 效应及一系列小地震的发生。但是这类现象并未在周边介质属性和地震波场的研究中得以验证。 为了探讨青藏高原深部物质流动的充要条件，我们基于多年来在高原进行的高精度人工源深部地震探 测和对地壳与上地幔结构的研究所给出的高原地区壳、幔介质属性和分层结构的平均模型(滕吉文，2011 年青藏高原学术会议报告)构建了数值模拟的模型。在地震波速度结构中，青藏高原上、中地壳间存在一 抬升。这便构成高原内部壳、幔物质运移的边界条件。 假设青藏高原岩石圈的基础由上地壳、中地壳、F地壳与上地幔盖层四层介质组成，二维数值模拟工 作构建了三组模型，分别考虑：①模型仅由这四层介质组成：②中地壳与下地壳之间加一软弱薄层：③下 地壳的上、下界面，即中地壳与下地壳之间，下地壳与上地幔盖层之间，均施加一软弱薄层。 研究结果表明，模型l中，运动的驱动力必须足够强大，且下地壳的