利用3-D+FEM%2bFDM方法模拟地壳中岩浆传播%3a以长白山火山为例.pdfVIP

·276· 中国地球物理201I (8)地球内部结构及其动力学 FEM+FDM方法模拟地壳中岩浆传播：以长白山火山为例 利用3．D of incrust 3D Numericalsimulation FEM+FDM magmapropagationusing —acase on volcano studyChangbaishan 程旭‘陈祖安 ChenXuChenZuan 中科院地质与地球物理研究所 北京市朝阳区北土城西路19号 100029 由地幔到地表的岩浆的运移和定位是一个重要并且困难的问题。本文通过数值模拟的方法，研究火山 喷发之前，岩浆从岩浆房开始在地壳中运移的动力学过程；研究岩浆房的超压力、岩浆性质(包括温度、 密度、粘度)和围岩性质(包括断裂韧性、粘弹性、密度、构造应力)对岩浆运移机制的影响：研究不同 条件下岩浆以多种不同类型速度与方向传播，会给地面带来什么样的不同信息(如地表变形和地震活动性 和各类火山地震事件等)，为预报火山喷发工作提供科学依据。 以往的运移机制包括整体的流体驱动破裂和局部的岩脉传播机制。来自熔体的浮力和压力的应力差和 局部的构造应力导致裂纹的张开与岩浆的传播。前人的工作中提出了各种流体驱动破裂的模型，本研究在 此基础上，考虑了黏性流体的影响，加入了浮力的影响，并用一个过压的岩浆源代替了恒稳的流量。本文 建立了地壳中岩浆传播的3维模型，利用有限元和边界元的方法计算出岩脉传播的形态和速度。边界元方 法用来分析裂隙中岩浆的流动，有限元方法用来计算岩脉传播过程中裂隙的打开。本文亦讨论了影响岩脉 传播的主要因素，以及相关的地质学意义。 建立模型时考虑从岩浆方中上升进入围岩的岩脉，看作是在有限空间中被岩浆所充满的3维边界裂纹。 裂纹的宽度看作常数，张开位移为W，裂纹的长度为a(O，随时间，变化。 假设岩浆为黏性流体，更具润滑理论，在流体充满的空间中，黏性流体的流动性行为有雷诺方程给出： Ow 1 a／，a， 、＼ 一Ot。面言I∥面【p+肪舻)J 其中，∥为岩浆的动力黏度，肪为岩浆的密度，胁lgz为岩浆受到的浮力，P为流体受到的差应力： t P Pe-(p,gz+ao) 仇为围岩弹性形变产生的应力，Pe=陆乜)，肛为围岩密度，A矿为垂直于岩脉平面的构造应力。 考虑压力P。(z)与张开位移啡)之间的弹性关系。对于非连续变形分析中的接触问题，当pe(z)已知， “z)可根据有限元方法求出。假设地壳中存在一断层(裂纹)，当裂纹尖端的压力P。0)大于构造虑力，初 始闭合的断层(裂纹)会张开，岩浆便会侵入进去。如此的过程会在一定条件下持续下去。包含闭合断层 (裂纹)的弹性问题转化为了一个接触问题。 根据已有研究成果，长白山火山并没有熄灭，也有可能在一定条件下再次活动，甚至喷发。本文以长白 山火山为应用对象，探讨它可能再次喷发的条件，并研究岩浆运移过程传到地面所得到的信息的特征。为监 测人员分析数据进行预报提供理论基础。因此，将理论和数值模拟研究应用到该地区，寻找火山喷发前的岩 浆运移过程所产生的前兆信号及其变化规律，对于预测火山喷发，减少灾害损失，具有十分重要的现实意义。 本研究由国家自然科学基金资助。 参考文献 andarrest simulationsofthe offluid·filledfracturesinthe 41： [1】Dahm,T．Numericalpropagation 623—五38．