论文摘要
复杂结构地壳非均匀介质的非弹性构造变形是近年构造物理研究的一个热点。地震和火山喷发等地质灾害除了浅部地壳活动外，也涉及深部岩石的粘弹性变形在时间过程中应力松弛并传递给浅部地壳弹性层的结果，除了考虑岩石圈纵向的不均匀还必须考虑横向的不均匀。火山区地壳(包括岩浆囊、通道、围岩等)介质是一种最典型的复杂结构的非均匀地球介质，岩浆囊是地球内部活跃的低粘滞性体，其围岩是粘弹性转换带，在应力传递过程中扮演重要角色，岩浆通道和裂隙的形成及火山区的地震是应力传递和释放的结果。长白山火山2002年至2004年出现了一个扰动期，火山区震群频发，GPS结果也显示出上拱膨胀的特点。火山的活动不仅与深部的活动有关，还受到浅部的构造变形的影响，因此除了理解岩浆活动的作用，还需要理解区域应力场和局部应力场共同作用下的火山区浅部地壳的活动特性。本论文以长白山天池火山地区为例，研究水平方向和垂直方向结构与介质不均匀的局部应力场特征。从概念模型出发，用低粘滞系数体代替岩浆囊，把火山区浅部地壳看作为非均匀粘弹性组合，分析含有岩浆囊的浅部地壳的粘弹性结构的应力分布、传递及其物理机制；同时模拟断层与低粘滞性体共同存在的条件下应力的调整；用粘弹性三维有限元模型，在区域性构造应力和局部岩浆的共同作用下，模拟火山区应力场特征，探讨粘滞性差异导致的局部应力异常以及与火山地震发生机制的关系，说明岩浆粘滞性在火山活动和区域应力调整中的重要影响。1、首先分析在没有岩浆活动情况下，横向粘弹性结构对应力场的影响。含有低粘滞性体的分层粘弹性组合的模型分析表明，在常速度条件下，应力传递不仅在垂向不均匀，而且在水平向也不均匀，它是由低粘滞性体与周围相对刚性的介质间的应变差异产生的，这种内部的不均匀性最终会在地表弹性层体现出来。由于介质粘弹性分层，弹性层中应力可以更快地集中，而低粘滞系数体的存在除了影响了局部应力分布外，还加快了局部弹性层中的应力积累速度。模拟结果说明在深度一定的情况下，低粘滞性体横纵比越大，对局部形变及应力场的影响就越小，反之，这种影响就越大。距地表越近的低粘滞性体引起区域应力场显著变化的可能性就越大，在浅部地壳，低粘滞性体赋存深度对局部区域的影响较之低粘滞性体形状大。低粘滞性体上方刚性壳体厚度影响着应力在其中的积累，厚度越小应力积累效果越明显，积累速度也越快。同时发现低粘滞性体达到一定深度以下时在上部壳体中的应力积累就基本没有差别，即对局部应力场和地表变形的影响趋于稳定。2、分析断层对含有低粘滞性体的地壳应力场的影响。模拟了具有不同的断层和低粘滞性体组合方式(断层分别位于低粘滞性体上方，低粘滞性体单侧和两侧)的模型的应力分布，发现断层对构造背景应力场的传递具有一定的阻隔作用，多数情况下断层的存在会减弱低粘滞性体对周围应力场的影响；当低粘滞性体两侧存在断层时，在距离低粘滞性体一定范围之外时，这种削弱作用较大；随着断层与低粘滞性体距离的缩短，这种削弱作用反而逐渐变小；另外，与没有断层存在的情况相比，断层的存在不仅会使局部应力场弱化，而且这种弱化作用会随着断层数量的增多而增强。3、利用火山区浅部地壳模型研究了岩浆活动时(岩浆压力作用下)的应力分布与变形特征。分析了2002年以来长白山固定台站以及流动台监测到的Ms≥1的地震数据，根据震级与时间的关系将长白山火山地震震群划分为了四个震群类型，同时发现长白山火山震群有周期发生的特点，大约每隔半年左右都会有一次大的震群。然而地震活动并没有持续活跃，2005至今长白山地震活动又趋于平静。因此推测2002年至2004年的地震活动可能是由于2002年6月29日汪清M7.2深震触发作用引发的一个较短时间的岩浆压力扰动过程，这种岩浆压力扰动过程与地壳和地幔介质的粘滞性应力传递活动有关，因为这种应力传递过程会通过耦合变形引起了上地壳应力场的调整，加上岩浆囊周围的差异应力，导致浅部地壳中形成网状裂隙，导致震群的发生。用三维粘弹性耦合模型模拟了在岩浆压力增大条件下，周围地壳应力场的变化特征，得到地表最大剪应力集中在岩浆囊上部以及岩浆囊端部对应上方表壳，这也解释了长白山地区火山地震发生在岩浆囊上方以及火山侧翼有很多喷泉(该区域最大剪应力大，易形成裂隙、通道)。4、考虑长白山火山地震是岩浆囊膨胀和区域(挤压)应力场的共同作用的结果，设计两个对比模型：只考虑区域应力场和考虑区域应力场与岩浆囊膨胀共同作用。模拟结果发现岩浆囊膨胀过程中考虑区域应力场的最大剪应力比不考虑岩浆囊膨胀应力场的大，这是由于岩浆囊膨胀过程中岩浆囊上方的区域东西向挤压应力场受到了抑制，而与垂直方向的挤压力却增大很多，即岩浆囊膨胀导致的局部应力与区域应力的共同作用会导致更大的剪切应力，这是岩浆压力扰动阶段火山地震多数为剪切型地震的一种解释。5、讨论了地壳与岩浆粘滞性比的影响及其在火山区变形过程中的作用。火山喷发的形式以及地表的隆起量都受粘滞性比的影响，地壳与岩浆粘滞性比值越小，岩浆上涌过程中地表的隆升就越显著，达到一定临界条件时大规模喷发的可能性就越大。通过粘弹性模型计算出了火山区变形的特点：上拱膨胀。同时用有限元模型与Mogi模型进行了对比：Mogi点源模型无法体现不同形状岩浆囊的结果，无法进行应力场分析，强调等效体积的作用弱化了其他参数的作用，而有限元模型可以完善这些不足。通过模拟，认为岩浆的活动与上地壳的形变关系可以分为三个阶段：①岩浆上涌在岩浆囊积累(压力积累，但此时囊内压力不能使上部地壳发生形变)→②压力大于外部压力时地表开始隆起→③囊内压力等于或小于外部压力时停止隆升，①～③阶段可视为一个周期。本论文的主要创新点在于：利用三维粘弹性有限元模型研究了岩浆压力作用下长白山火山区的非均匀地壳应力分布，并用以解释长白山天池火山区地震的时空分布及其机制特征，探讨了岩浆粘滞性、岩浆囊、岩浆活动在火山区域构造中的作用。研究了局部低粘滞性体对区域构造的重要影响，在分层粘弹性模型的基础上加入了横向非均匀结构的影响，通过模拟低粘滞性体和一些构造元素(断层以及不同边界条件)的组合，说明组合体内的应力分布不均匀以及差应力的突变可能是引发板内局部区域地震的动力因素。