燕山地区早侏罗世岩浆活动热供给的数值模拟 Numeric Simulation of Early Jurassic Magma-heating Supply in the Yanshan Area.pdfVIP

第33卷第4期 地 震 V01．33，No．4 2013年10月 EARTHQUAKE ()ct．。2013 燕山地区早侏罗世岩浆活动 热供给的数值模拟+ 刘 翠1’2，石耀霖2，乔彦超2，邓晋福1，李宁1，段培新1 100049) (1．中国地质大学，北京100083；2．中国科学院计算地球动力学重点实验室。北京 摘要：燕山地区燕山期发育大量的岩浆岩，如此大规模的岩浆活动(热)的来源一直是一个谜。 本文利用有限元的热传导模拟反演可能的热供给。模型假设温度1250℃、面积905．86km2半 圆面形玄武岩浆底侵到36．6～50 km深的陆壳底部，陆壳的地温梯度随着克拉通的破坏而逐 渐升高。数值模拟的主要结论：①早侏罗世时期的地温曲线表明，要使围岩发生熔融，所需要 的底侵玄武岩岩浆量是非常大的，基于燕山地区实际产出的早侏罗世时期的酸性岩，反演产生 81 053 km2的花岗岩需要的最少底侵玄武岩岩浆量为1km!；②早侏罗世(J。)的模型温度场表 明，随着时间的演化，在40～50 km处水平方向的围岩温度升高明显，形成很明显的一条高温 带或熔融带，该高温带或者熔融带可能成为后期的构造薄弱面。加上榴辉岩相变造成的密度增 加引起下坠，达到临界值时，就会出现垮落，因此本文从数值模拟的角度，支持拆沉作用的发 牛。 关键词：岩浆活动；早侏罗世；拆沉；数值模拟；燕山地区 中图分类号：P315 文献标识码：A 文章编号：1000一3274(2013)04一0257一12 引言 燕山地区是燕山运动最为典型的地区，是华北地区燕山期克拉通活化最为理想的研究 区之一。“华北克拉通破坏”是目前我国乃至全球研究的热点，有关学者对此进行了深入、 详细的研究和探讨，包括地质、地球化学、地球物理等‘1”]。而数值模拟作为弥补无法直观 探测的困难和海量数据的需要，在探讨华北克拉通可能的破坏机制以及相关参数的获得方 面正起到越来越重要的作用[0~1“。。 目前对于华北克拉通破坏的若干问题仍然存在很大的争议，其中之一就是燕山期华北 克拉通破坏的机制到底是热侵蚀作用还是拆沉作用。1“’1”16j。对于该区燕山期的地质演 化，拆沉作用能否发生，具不具备相关的数学地质条件；另外，燕山地区在J—K时期，发 7~1 育有大量的岩浆岩Ll 9j，邓晋福等[13将早侏罗世(J，)到早白垩世(K，)共划分了4幕造山 *收稿日期：2013一02—18；修改回日期：201305—28 基金项目： 质调查局项目(1212011121075；12】2010911028)及国家留学基金委资助 作者简介： 刘翠(1973一)，女，山东邹城人，博士，主要从事岩浆岩及岩浆数值模拟等研究。 万方数据 地 震 33卷 过程(不含早白垩世晚期(K})的伸展阶段)，其中早、中侏罗世(J。和J。时期)的岩浆岩量 比较少，主要为中基性的玄武岩、安山岩类，而J。和K}发育大量的酸性岩。因此，为何该 区在燕山期随着时间演化，岩浆岩的量越来越多，并且越来越偏酸性。上述这些都是需要 回答的问题。 本文结合燕山地区出露的岩浆岩，对早侏罗世(J。)的岩浆活动开展二维数值模拟，以 更加切合实际的限定条件来探讨燕山期华北克拉通破坏的热机制和来源。 1原理、方法及参数 1．1 数值模拟的原理及方法 地球内部热的传播方式包括热传导、热对流、热辐射等。热传导是介质内无宏观运动 时的传热现象，其在固体、液体和气体中均可发生。但严格而言，只有在固体中才是纯粹 的热传导，而流体即使处于静止状态，其中也会由于温度梯度所造成