幔源岩浆在地壳中分层侵位的控制因素:二维热-力学模拟.pdf
《幔源岩浆在地壳中分层侵位的控制因素:二维热-力学模拟.pdf》由会员分享,可在线阅读,更多相关《幔源岩浆在地壳中分层侵位的控制因素:二维热-力学模拟.pdf(15页珍藏版)》请在咨信网上搜索。
1、书书书 ():岩石学报 :崔晓娜,陈林 幔源岩浆在地壳中分层侵位的控制因素:二维热 力学模拟 岩石学报,():,:幔源岩浆在地壳中分层侵位的控制因素:二维热 力学模拟崔晓娜,陈林 ,中国科学院地质与地球物理研究所,岩石圈演化国家重点实验室,北京 中国科学院大学地球与行星科学学院,北京 ,收稿,改回 :,():,:,(,):()(),(),(),()(),;摘要幔源岩浆在地壳内的上升和聚集样式不仅依赖于岩浆自身的性质,还取决于围岩的强度和热状态。已有数值和物理模型大多关注岩浆自身物性对其上升过程的影响,而对围岩流变强度或热状态如何影响岩浆的上升和聚集过程的研究相对薄弱,尤其是周期性的幔源岩浆在壳
2、内分层侵位的受控因素仍然不清楚。本文利用二维热 力学数值模拟方法,通过发展本文受国家自然科学基金项目()和国家重点研发计划项目()联合资助第一作者简介:崔晓娜,女,年生,博士生,地球动力学专业,:通讯作者:陈林,男,年生,研究员,主要从事地球动力学数值模拟研究,:多期岩浆脉冲和岩墙生成算法,研究了岩浆从深部地幔上升至地壳内部侵位的动力学过程,系统测试了地壳(围岩)地温梯度(上、下地壳的地温梯度分别以 和 表示)和地壳强度对岩浆上升过程和聚集样式的影响。模拟结果表明:()地壳地温梯度对岩浆的侵位深度有重要影响,岩浆侵入冷地壳(),岩浆主体在岩石圈深度聚集,地表的相对高差小于 ;岩浆侵入温地壳()
3、在下地壳底部聚集形成岩浆房,上升至上、下地壳界面,岩浆房上方的地表地形呈现中心拗陷两翼隆起的形态,地表最大高程可达 ;岩浆侵入热地壳()仅在下地壳底部聚集形成岩浆房,地形演化特征与温地壳背景的情况类似,但最大高程小于 。()在同等地壳地温梯度条件下,上、下地壳的相对强度决定了岩浆的聚集样式:下地壳的强度越弱,岩浆更易在下地壳聚集形成岩浆房;上地壳的强度越弱,岩浆更易在地表喷发。进一步分析表明:岩浆的聚集样式受地壳地温梯度与地壳流变分层性的共同控制,地壳越热且流变强度分层性越显著则越有利于岩浆在地壳中的多层级侵位;每一期岩浆脉冲的供给均会导致岩浆房内部的超压值骤增。我们的模拟结果对理解火山喷发前
4、壳内岩浆的赋存状态及岩浆活动区的地形演化具有启示意义。关键词岩浆;地壳;地温梯度;粘度;数值模拟中图法分类号 上地幔部分熔融形成的低粘度、低密度的基性岩浆在地壳中聚集和侵位的过程通常会伴随岩浆房、岩墙或岩席的形成(,;,)。岩浆的侵位和喷发往往会导致地震、海啸等自然灾害的发生(,;,),也会形成铜矿、金矿等金属矿床(,;杨林等,)。因此,理解幔源岩浆从深部上升至地表的动力学过程具有重要的科学意义。高分辨率地球物理成像结果揭示,我国东北长白山火山下方的地壳内部存在多层级的岩浆聚集区(,;,)。岩石地球化学数据表明,地表出露的玄武岩呈现不同期次的岩浆补给(郭文峰等,;,)、岩浆房多次聚集的特征(,
5、;,),且在地表喷发的火山数量远小于地下赋存岩浆的体积(,)。由此可见,幔源岩浆在地壳中的侵位具有多期次岩浆供给和分层聚集的特点。然而,幔源岩浆侵位过程中产生的分层聚集样式受何种因素控制仍不清楚。针对幔源岩浆的上涌过程,前人已开展了大量数值和物理模拟研究。早期数值模拟研究主要关注岩浆侵位对围岩带来的热影响(,;,),或是探讨单期幔源岩浆在均质地壳的聚集过程(,;,)。近期的数值模拟研究主要关注岩浆的传输通道(,;,)。围岩的流变强度及其分层性决定了岩浆的传输方式,包括熔融底辟、减压通道或张性裂隙(,),而岩墙和底辟可以共同促进岩浆向上传播(,)。伴随岩浆侵位,侵入体周围发育的裂隙由初期呈径向展
6、布的剪裂隙转换为向上扩展的主干张性断裂(卢靖雯等,),同时岩浆侵位会对围岩施加附加变形,导致近场与远场的小尺度变形(常成和罗纲,)。()利用三维热 力学模拟方法探究地壳流变和温度结构对基性岩侵位形态的影响,最近 ()考虑了岩浆熔融抽取过程中的岩浆成分演化,模拟壳内多期次玄武岩席聚集形成岩浆房的机制。物理模拟研究早期主要探究岩墙向岩席的转换机制,通过将流体注入固体明胶的模拟实验,验证了当固体明胶存在刚度对比层、层间弱接触面或外部应力时,可以导致垂向传播的裂隙运动轨迹转为水平向(,;,)。最近,()利用物理模拟方法探究了岩浆沿岩席长距离水平传播的形成机制,发现岩浆的流变学性质具有重要影响。然而,物
7、理模拟受相似比的限制,难以直接运用于理解实际的地质过程(,)。综上,前人的相关研究要么没有考虑岩浆的多期补给,要么没有考虑岩墙作为不同层级岩浆房之间的连接通道,尤其是针对多期幔源岩浆在地壳中分层侵位过程及其控制因素的研究仍然很薄弱。在本研究中,我们聚焦幔源岩浆从深部到地壳的聚集运移过程,通过在二维热 力学数值模拟程序中引入多期岩浆脉冲和岩墙生成算法,模拟幔源岩浆从深部地幔上升到地壳的动力学过程,系统测试地壳地温梯度和流变强度的影响,探究幔源岩浆在地壳内部分层侵位的控制因素。我们将首先介绍数值模拟方法与模型设计;然后描述地壳地温梯度和地壳强度对壳内岩浆赋存状态的影响;最后,讨论幔源岩浆在地壳内聚
8、集样式的控制因素。数值模拟方法我们使用开源的二维热 力学 程序(,)模拟幔源岩浆多期次侵位的动力学过程。本程序在求解动量守恒方程、质量守恒方程与能量守恒方程时,联合交错网格有限差分方法和粒子标记(,)技术。在节点上解方程,在 上传递物性参数,可以有效避免利用有限差分方法求解连续性方程时可能引起的数值发散问题。同时,该方法考虑了岩石的粘 弹 塑性流变本构关系和岩石的部分熔融,可以在时间与空间尺度上精细模拟岩浆侵位的动力学过程。控制方程()质量守恒 连续性方程 岩石学报 ,():()式中,和 分别代表速度矢量的水平和垂直分量。()动量守恒 斯托克斯方程 ()(,)()式中,和 分别代表水平和垂直坐
9、标;是密度,为压力 、温度 、岩石成分 和部分熔融百分数 的函数;是重力加速度;、分别为偏应力张量,其本构关系为:()()()式中,、为应变率张量;代表有效粘度。()能量守恒 热传导方程 (,)(,)()()()式中,是等压热容;代表温度对时间的物质导数;和分别代表水平和垂直方向的热流;是热导率,与压力 、温度 和岩石成分 有关(,);是热膨胀系数;代表放射性生热率,依赖于岩石成分;代表剪切生热,表示形变过程中由机械能转化的热能;代表绝热生热,表示随着压力变化而产生的热量;代表潜热,表示岩石发生相变过程中释放或消耗的热。岩石流变本程序考虑了岩石的粘 弹 塑性()流变,体偏应变率 包括三个分量:
10、()()()()()(),()式中,为剪切模量;为偏应力的目标共旋转时间导数;为岩石的屈服应力;为第二偏应力不变量 ();为塑性势;为塑性因子。以上方程中的偏应力分量 是基于粘 弹 塑性本构关系以显式的一阶时间有限差分形式表示为粘 弹性应力的目标时间导数:(),()式中,为弹性时间步长;是上一个时间步长由 技术求解并经过平流和旋转后的偏应力张量(,;,);为粘 弹因子;为有效的粘塑粘度。()()其中,式中,是气体常数;、和 是岩石物理实验约束的流变参数,分别代表物质常数、活化能、活化体积和应力指数;为 粘 度 因 子,用 于 线 性 改 变 物 质 的 有 效 粘 度(,)。由于实验室条件下物
11、质成分的变化,所测得的流变参数也具有很大的不确定性。为了将这些结果应用到实际地质过程,我们使用粘度因子对物质在实验室条件下测得的有效粘度进行线性调节(刘丹红和陈林,)。在这里,我们通过改变粘度因子调控地壳的强度。熔融岩石(,为岩石的部分熔融百分数)的粘度为 。一方面,我们假定岩浆上升是一个十分快速的地质过程,在该过程中没有明显的热损失,因此粘度基本没有变化;另一方面,将部分熔融岩石的粘度取为常数,可以简化一些复杂因素的影响,便于我们重点考察围岩的粘 弹 塑性流变对幔源岩浆上升过程带来的影响。此外,在模拟过程中,我们设定粘度的上限为 ,下限为 ,采用的流变参数见表 。物质的塑性变形遵从 屈服准则
12、(,):()()其中 是岩石的内聚力;是有效内摩擦角。崔晓娜等:幔源岩浆在地壳中分层侵位的控制因素:二维热 力学模拟表 数值模型中用到的主要物理参数 物质流变)()()()()()()()()()()()空气 上地壳湿石英 ()下地壳斜长石 ()岩石圈干橄榄石 ()软流圈干橄榄石 ()熔体 ()基性岩浆 ()参考文献),注:其他通用物质参数:热膨胀系数 ,压缩系数 ;)()(),式中 代表岩石在 和 条件下的密度;)();();();)当 ,()();当 ,;当 ,()();当 ,;当 ,;当 ,();当 ,();当 ,;)();();();();();();();();();()部分熔融根据
13、岩石物理实验数据的约束,程序考虑了岩石的部分熔融行为(,;,)。假设部分熔融百分数 与温度为简单的线性关系:,()式中,为岩石的固相线,为岩石的液相线,各种岩石的固相线和液相线见表 。部分熔融岩石的有效密度 取决于部分熔融百分数:()()式中,和 分别代表岩石的固相密度和熔融密度,它们都是关于温度和压力的函数,计算表达式见表 。其中,我们对幔源岩浆密度的选取参照了 ()的结果,并且与该作者汇编的结果进行了对比,在 的压力范围内,我们的结果中岩浆的密度与前人的实验结果和理论计算结果均一致(电子版附图 ),说明岩浆密度的选取具有合理性。部分熔融对有效热容量()和热扩散系数()的影响由下式给出(,)
14、:()()()()()多期岩浆脉冲的实现火山喷发需要体积足够大、低粘度、低密度、部分熔融且持续供应的岩浆热源,以提供浮力、维持高温(,)。我们通过预设岩石圈之下的岩浆房来模拟幔源岩浆上升的动力学过程,在已有的二维粘 弹 塑性数值模拟程序(,)的基础上加入多期岩浆脉冲算法,在岩石圈以下多期次、定量的稳定提供岩浆热源,尽量保证岩浆通道不变形。算法实现细节为:()我们添加新一期岩浆脉冲的判断条件为:深度 以下,部分岩浆通道和岩浆源区中熔体的总含量小于 ;()当熔体含量满足上述条件时,我们用性质相同但标号不同的 补充新一期次的幔源岩浆;()在模拟过程中,我们总共提供六期岩浆脉冲,为岩浆的持续上涌提供动
15、力。岩墙的生成算法岩墙作为拉伸作用形成的薄弱带,它的形成与围岩的应力场有关(,;,)。在本文中,拉伸应力标记为正、挤压应力标记为负。岩墙作为岩浆房之间的连接通道(,;,),岩浆沿岩墙传播(,;,),其宽度随深度变化(,)。目前对于岩墙的热 力学模拟主要探究岩墙对岩浆上升过程的影响,例如,()认为花岗岩岩浆房会形成辐射状岩墙,为了探究岩浆上升过程中岩墙和底辟的相互作用,因此设定临界的应力和第二应力不变量作为垂直状岩墙形成的条件。()通过岩石热 力学模拟方法,考虑岩墙 岩石学报 ,():图 岩墙形成的示意图 作为玄武岩熔融提取的通道。我们通过探究幔源岩浆侵位对地壳应力状态带来的影响,认为垂直状岩墙
16、是岩浆上升的通道,辐射状岩墙是岩浆房体积扩张时围岩的同期破裂,表征岩浆房的发育。因此,本文中考虑岩墙是在岩浆侵位对上覆围岩造成拉伸环境中形成的(,),作为岩浆向浅部运移的通道。我们在程序中加入岩墙生成算法,岩墙的形成和扩展由岩浆房超压值和地壳应力共同控制。岩浆房超压值为岩浆房内部的总压强与静岩压力的差值(岩浆房内部的超压值相差无几),岩浆房内部的超压值会随着每一期次的岩浆补给而增加、随着岩浆房的体积增大而减小,超压值演化表征了岩浆房的发育过程;当岩浆底侵时,岩浆聚集使得岩浆房之上产生应力集中的区域,即岩浆房的发育对上覆的上地壳与下地壳形成了水平的拉伸应力,这种拉伸应力即为地壳应力,地壳的应力场
17、决定了岩墙垂直向发育。以下地壳岩墙的形成过程为例,具体算法如下:()首先,由于多期次的岩浆脉冲供应,壳内岩浆房发育至一定体积,将岩浆房内部超压值小于 作为岩浆房发育结束的时刻与判断是否形成岩墙的初始时刻。()随后,当岩浆房之上的地壳应力大于 时,发育长度为岩浆房顶部至上、下地壳界面、宽度为 的岩墙(见图)。经过多次测试,宽度的岩墙是保证岩浆上升的最窄宽度,过窄的岩墙会导致通道快速冷却、无法输运岩浆,因此我们在模拟中设定岩墙的最小宽度为 。模型设计为了模拟幔源岩浆上升的动力学过程,我们设计了如图所示的初始模型。模型在水平方向上的宽度为 ,垂向深度为 ,包含 上地壳、下地壳和 岩石圈地幔,代表了一
18、个典型的大陆岩石圈,具体物性参数见表 。我们初始在岩石圈底部预设一个半径为 、温度为 的圆形幔源岩浆区域,在顶部设置了一个宽度为 的岩浆通道,连接到地壳底部。整个模型空间被 个节点组成的规则网格离散化,空间分辨率为 ,总共 万个粒子()随机分布在网格空间中,粒子分辨率高达 。模型的所有力学边界均为自由滑边界。在模型的上部,我们设计一个厚 的伪空气层,其密度和粘度分别为 、,用于模拟岩浆上升过程中产生的地表起伏(,;,)。模型的热边界条件:顶边界固定为 ,侧边界为绝热边界,底边界为 固 定 为 (,;,),整个空气层温度固定为 ,岩石圈底边界以下为绝热地温梯度 。为了区分地壳的不同热状态,我们采
19、用分段函数赋予了地壳两种初始温度结构:第一种温度结构初始固定莫霍面的温度,设定地表和莫霍底界面的初始温度分别为 和 ,通过改变上、下地壳界面的初始温度(、)来调节上、下地壳的地温梯度(、),让温度随深度线性增加;第二种温度结构初始固定岩石圈底部的温度,设定地表和岩石圈底界面的初始温度分别为 和 ,通过改变莫霍面的初始温度(、)来调节地壳和岩石圈的地温梯度,同样让温度随深度线性增加。模拟结果为了探究来自地幔深部的多期岩浆在地壳中分层聚集的动力学过程,我们基于多期岩浆脉冲和岩墙算法的 程序,开展了 组数值模拟实验(表 )。在模拟中,系统测试了地壳地温梯度和地壳强度对岩浆上涌过程的影响。在这里,我们
20、选择上地壳的粘度因子 、下地壳的粘度因子 、的模型作为参考模型(即表 中的 )。参考模型图 为参考模型(,)的演化结果。来自岩石圈底部的幔源岩浆穿过岩石圈到达莫霍面,初始阶段(;图 ),岩浆在莫霍面之上聚集形成椭圆状岩浆房,伴随同时的垂向增厚和横向扩展,粘度结果显示岩浆房周围出现局部的低粘度区域(图 ),表征岩浆房的扩张趋势。随后,由于岩浆侵位对上覆地壳造成的拉伸力,下地壳形成岩浆通道使得岩浆由此上升(图 ),水平拉伸应力瞬时消散(图 )。时,岩浆上涌至上、下地壳界面,多期次岩浆脉冲带来的高超压导致岩浆通道不断变宽,在岩浆房内部,二期岩浆及多期岩浆以韧性底辟形式对流,促进多期岩浆混合及岩浆崔晓
21、娜等:幔源岩浆在地壳中分层侵位的控制因素:二维热 力学模拟图 模型设计()为模型全空间;()为()中黑色方框的放大区域岩性色标见标注,不同的紫色表示六期基性岩浆,不同的绿色表示六期熔体,蓝色为冷却后的基性岩浆,红色为冷却后的熔体 ();(),表 数值模拟实验 模型 熔体含量()()()()相关图件 图 、附图 附图 附图 附图 图 、附图 附图 附图 附图 附图 附图 附图 岩石学报 ,():续表 模型 熔体含量()()()()相关图件 附图 附图 图 、附图 附图 图 、附图 附图 附图 附图 附图 图 图 图 房上升(图 ),粘度结果显示岩浆房之上局部的低粘度区域(图 ),表征岩浆的上升趋
22、势。最终 时,上、下地壳界面之上岩浆聚集形成扁平状岩浆房,但岩浆房顶部的应力累积不足(图 ),达到稳态至形成双层聚集的岩浆房,此时下地壳岩浆房体积扩张,并且内部卷入极少量的地壳和岩墙等物质(图 ),上地壳与下地壳减薄程度不同的原因在于两者对于岩浆侵位带来的空间问题响应不同。在 的时间尺度内,岩浆房内部保持活跃,进行持续的岩浆对流。图 中 为参考模型的水平偏应力随时间的变化,下地壳与上地壳岩浆侵位均对岩浆房之上的局部区域形成水平拉伸应力。上地壳累积的拉伸应力在脆韧性转换深度处最大,向地表逐渐降低,向深部为韧性区应力累积很弱;下地壳累积的拉伸应力在上、下地壳界面处最大,向深部减小,同样韧性区应力累
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 岩浆 地壳 中分 层侵位 控制 因素 二维 力学 模拟
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。