论文部分内容阅读
关于华北克拉通岩石圈减薄方式,国内外学者通过各种研究定性的提出了不同的机制。本篇论文根据最近几年国内外发表的有关资料,通过数值模拟方法,主要对地质学家提出的华北克拉通岩石圈热对流侵蚀和拆沉减薄机制分别进行了初步模拟。模拟的目的主要在于探讨不同机制发生的条件、特点、可能造成的地质过程有什么差别,从而为探讨华北克拉通减薄问题提供参考。 论文主要围绕华北克拉通岩石圈减薄机制进行以下研究:①针对燕山期华北克拉通的岩浆活动进行了数值模拟工作,考虑在热传导的机制上,玄武岩浆底侵围岩,使其部分熔融,部分熔融的围岩最终变为花岗岩,现今华北克拉通有大规模的花岗岩区,我们通过统计花岗岩的量可以反演底侵的玄武岩浆的量,而这部分玄武岩在岩石圈加厚后,可以相变为密度加大的榴辉岩,通过这种方式估计出的榴辉岩的量可以约束我们计算华北克拉通岩石圈拆沉时榴辉岩的量,在计算过程中,我们考虑了相变对这个过程的影响,更符合实际情况。②根据地质学家关于华北克拉通岩石圈减薄机制的理论,我们数值计算了岩石圈热侵蚀减薄机制的形成条件和特点,主要是采用了上地幔小尺度方腔地幔对流的程序,考虑了粘滞系数随温度变化,模型中热侵蚀在上地幔底部有热扰动,使其温度增加时,会造成上地幔系统活化,对流加剧而造成岩石圈的减薄,我们讨论了两种等效粘滞系数情况下,不同上地幔底边界温度对计算结果的影响,并得到了每种情况的岩石圈减薄曲线。③数值计算了岩石圈拆沉减薄机制的产生条件和特点,拆沉机制认为岩石圈先挤压加厚下地壳玄武岩相变为榴辉岩,而榴辉岩的密度大于地幔的平均密度,这样下地壳的重力不稳定造成上地幔系统的不稳定对流,使岩石圈发生减薄破坏,通过我们的计算讨论了不同规模的榴辉岩拆沉,会有不同的模式,我们也计算得到了不同规模榴辉岩造成的岩石圈拆沉量。 通过以上对华北克拉通岩石圈减薄机制的数值模拟研究,主要得到以下结论:①通过计算表明在1GPa条件下,地幔的玄武岩浆底侵作用于下地壳英云闪长岩围岩的情况下,英云闪长岩围岩纯的熔融量与底侵玄武岩纯的结晶总量的比值约为0.161;假设仅有25%的熔浆可分离出来,那么,可分离出的熔浆与玄武岩结晶总量比为约0.04025。即如果生成总量为1km厚的花岗岩的岩基需要24.85km的底侵玄武岩岩浆的结晶来提供热量。而24.85km的底侵玄武岩岩浆需要约248.5km的地幔物质参与(假设二辉橄榄岩10%的熔融程度,并且全部局部熔浆均可分离出来)。由于玄武岩浆内热对流过程传热更有效率,可以从深部提供大量的热量,所以实际的玄武岩厚度可能远低于24.85km。数值模拟结论得到最后时刻发生相变的整体区域为-36.867到-50.285km。相比于樊祺诚等人对汉诺坝捕虏体研究反演得到的现今下部下地壳和壳幔过渡带33-40和40-50km,也比较合理,因此得到了岩石学验证。②通过我们的计算证实了原来稳定存在的克拉通,如果发生在上地幔底部温度扰动升高的情况,由于浮力驱动的小尺度地幔对流加剧能够使岩石圈发生大规模减薄,减薄速率在mm/y的量级,减薄的时间尺度基本都在十几个百万年;计算过程中我们讨论了初始参考等效粘滞系数分别为η0=1.0E22Pa s和η0=1.0E23Pa s的两种情形,在这两种条件下,我们分别计算了底边界温度从原来1673K增加为1773K,1873K,1973K,2073K的4种情况。通过不同端员的计算我们知道岩石圈可以从200km减薄到100km-126km左右,这符合现今地球物理观测的结果。初始参考等效粘滞系数和底边界温度是影响减薄速率的重要因素。③通过本文的研究,证实中生代华北克拉通岩石圈受挤压,如果有足够多的玄武岩相变为密度较大的榴辉岩,则榴辉岩重力失稳能够引起下地壳的拆沉,造成岩石圈的大规模减薄。通过我们对比不同规模的榴辉岩减薄方式发现,当榴辉岩较大规模时更可能发生双管道拆沉,而规模小时发生单管道拆沉。我们对比不同底边温度的结果发现,虽然热减薄和拆沉减薄机制,从力学机制上都能造成岩石圈减薄,但数值模拟提示它们可能对应着不同的地表地质特征。热对流减薄,中心区域变热变薄,中心区域是一个伸展区域,发展一个大的变质核杂岩区,岩浆的发展时序是从中心向两侧对称分布,并且向两侧减弱;拆沉减薄可能会形成中心区域受到挤压加厚,反而变冷。在中心区域的两侧发生两个减薄伸展区域,可能对应两个变质核杂岩区域,岩浆的发展时序是从两侧向中心减弱发展。这一结果对地质上研究华北克拉通减薄机制具有参考价值。