华北克拉通燕山期岩浆活动的数值模拟

摘要

为了定量给出华北克拉通燕山期岩浆活动的规模,在合理构建华北燕山期地质-物理模型的基础上,利用有限元方法,在忽略玄武岩浆内对流只考虑热传导的假定下,反演了该区燕山期形成花岗岩所需要的侵入玄武岩浆的量.通过计算表明:（1）在压力1GPa、地幔的玄武岩浆底侵作用于下地壳英云闪长岩围岩的情况下,英云闪长岩围岩的纯熔融量与底侵玄武岩纯结晶总量的比值约为0.161;假设仅有25%的熔浆可分离出来,则分离出的熔浆与玄武岩结晶总量比约为0.04025.（2）数值模拟得到最后发生相变的整体区域为-50.285～-36.867km.与前人研究得到的现今下地壳下部和壳幔过渡带33～40m和40～50km接近,表明此数值模拟方法具有可行性.